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CHAPTER  I 
INTRODUCTION 

The  moisture  content  of  air  has  a  major  influence  on  many  industrial 
operations  conducted  today.  Because  of  this,  it  is  often  necessary  to 
limit  the  moisture  content  of  air  on  a  year-round  basis. 

Dehumidification,  the  removal  of  water  from  air,  is  essential  in  the 
handling  of  a  variety  of  products.  With  the  storage  of  materials  and 
products,  maintaining  a  low  air-moisture  content  level  retards  the 
formation  of  corrosion  and  may  eliminate  the  need  for  surface-protecting 
materials.  In  the  food  industry,  a  reduced  humidity  level  slows  the  growth 
of  molds  and  increases  product  storage  life.  In  manufacturing,  the  use  of 
dehumidified  air  during  critical  processes  often  results  in  a  higher 
quality  final  product.  Also,  worker  productivity  is  related  to  personal 
comfort,  which  is  influenced  by  the  moisture  content  of  the  air. 

The  actual  dehumidification  process  depends  on  the  temperature  and 
pressure  relationships  which  exist  between  the  air,  its  moisture  content, 
and  surrounding  objects.  Air,  in  general,  is  composed  of  a  mixture  of 
different  gases  with  each  exerting  an  individual  partial  pressure  which, 
when  summed,  equals  the  total  overall  pressure  of  the  mixture.  Moisture 
contained  in  air  is  in  the  form  of  a  vapor  and  therefore  also  maintains  a 
partial  pressure  as  part  of  the  overall  mixture  pressure.  This  partial 
pressure  cannot  exceed  the  saturation  pressure  of  water  for  the  temperature 
of  the  mixture.  Because  saturation  pressure  varies  directly  with 
temperature,   it  is  possible  to  limit  the  moisture  content  of  air  by 


lowering  the  temperature  of  the  mixture,  forcing  some  of  the  water  to 
condense  out.  Also,  differences  in  vapor-pressure  will  drive  moisture  from 
an  area  of  high  vapor-pressure  to  an  area  of  low  vapor-pressure. 
Therefore,  dehumidification  may  also  be  accomplished  through  the  use  of 
sorbent  materials,  which  utilize  a  low  vapor-pressure  to  extract  excess 
moisture  from  the  air. 

Commercial  dehumidification  of  air  in  most  instances  is  accomplished 
by  the  use  of  one  or  both  of  the  following  types  of  systems:  desiccant 
systems  and  refrigeration  systems.  Both  methods  have  achieved  widespread 
popularity  in  industry  and  the  selection  of  the  method  to  be  utilized  often 
depends  on  the  industrial  operation  to  be  performed.  Desiccant  systems, 
which  utilize  sorbent  materials,  are  best  suited  for  operations  which 
require  high  air  temperature-low  moisture  content  levels.  Refrigeration 
systems,  which  dehumidify  by  cooling  the  air,  are  often  utilized  in 
situations  were  air  of  high  moisture  content  must  be  processed  and 
relatively  low  final  air  temperatures  are  required.  Combinations  of  the 
two  systems  are  frequently  used  to  handle  air  processing  situations  which 
fall  between  these  two  extremes. 

Both  of  the  methods  presented  above  are  capable  of  handling  large 
quantities  of  air  for  dehumidification.  The  selection  of  the  actual  system 
to  be  utilized  is  often  based  on  a  required  exit  air  condition,  initial 
cost,  and  the  energy  requirements  of  the  system  at  design  conditions.  This 
procedure,  however,  may  not  select  the  most  energy  efficient  system 
available  because,  on  an  annual  basis,  the  majority  of  the  system 
operational  time  is  spent  at  conditions  far  less  severe  than  those  of 
design.  As  a  result,  a  method  for  predicting  system  energy  requirements  on 
an  annual  basis  should  be  included  in  the  selection  process.   To  accomplish 


this,  ASHRAE   commissioned  Research  Project  298  RP  "Analysis  of  Open 
Sorption  and  Refrigeration  Humidity  Control  Systems"  . 

The  objective  of  ASHRAE  Research  Project  298  RP  was  to  develop  an 
analytical  method  for  estimating  the  annual  energy  consumption  and  cost- 
effectiveness  of  the  different  dehumidification  systems  available, 
including  refrigeration,  open-sorption  (desiccant)  dehumidification,  and 
combinations  of  both  refrigeration  and  open-sorption  dehumidification. 
From  this,  the  optimum  dehumidification  system  for  the  desired  application 
may  be  determined. 

1.1  OBJECTIVE 


The  objective  of  this  study  was  to  develop  mathematical  models  for 
predicting   the  energy  requirements  of  refrigeration   systems   when 
subjected  to  variations  in  both  dehumidification  load  and  ambient  air 
conditions  which  would  influence  system  operation  on  a  year-round  basis. 
The   study   focused   primary   on   reciprocating   vapor-compression 
refrigeration  systems   and  was  performed  in  support  of  ASHRAE  Research 
Project  298  RP,  "Analysis  of  Open-Sorption  and  Refrigeration  Humidity 
Control  Systems"  . 
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Treatment  of  open-sorption  dehumidification  systems  may  be  found 
in  Reference  [1 ]. 


CHAPTER  II 
REFRIGERATION 

Dehumidification  of  air  by  refrigeration  requires  the  removal  of 
sensible/latent  energy  from  the  air  stream  being  processed.  For  this 
process,  the  air  must  be  exposed  to  a  low  temperature  medium  for  a 
significant  period  of  time,  allowing  the  transfer  of  enough  energy  to  take 
place  between  the  air  stream  and  the  medium  for  the  air  to  be  cooled  to  a 
predetermined  level.  At  this  level,  the  air  is  capable  of  supporting  less 
moisture  than  before  and  any  excess  moisture  is  forced  to  condense  out  of 
the  mixture  to  insure  that  the  saturation  pressure  of  water  is  not 
exceeded.  A  typical  refrigeration  process  for  the  cool ing/dehumidification 
of  air  is  illustrated  on  a  psychrometric  chart  in  Figure  2.1. 

Refrigeration  systems  may  be  subdivided  into  two  general  categories: 
Absorption  and  Vapor-Compression.  Both  methods  utilize  a  low  pressure- 
temperature  fluid  to  absorb  heat  from  the  desired  medium  (air)  and  a  high 
pressure-temperature  fluid  to  later  reject  the  heat  to  another  medium.  The 
difference  between  the  two  systems  is  in  the  type  of  fluids  utilized  and 
the  way  in  which  the  high  and  low  pressure  sides  of  the  systems  are 
generated. 

2.1  ABSORPTION  REFRIGERATION 


In  the  absorption  system  (Figure  2.2),  refrigerant  vapor  at  a  low 
pressure  is  absorbed  in  a  solution  (the  absorbent)  and  the  resulting 
solution  is  pumped  to  a  high  pressure  container  (the  generator).   At  the 
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Figure  2.2  Schematic  for  Absorption  Refrigeration  systems. 


generator,  heat  is  added  to  the  solution,  forcing  the  absorbent  to  reject 
the  refrigerant.  The  liquid  absorbent  is  then  returned  to  the  absorber 
while  the  refrigerant  vapor  is  sent  to  the  condenser.  At  the  condenser, 
the  refrigerant  rejects  energy  to  the  surroundings  and  condenses.  The 
liquid  refrigerant  then  flows  from  the  condenser,  through  a  pressure 
reduction  valve,  and  into  the  evaporator.  Here,  heat  is  absorbed  from  the 
desired  medium,  allowing  the  refrigerant  to  revaporize.  Finally,  the 
refrigerant  returns  to  the  absorber  to  repeat  the  process. 

For  the  above  mentioned  system  to  operate  efficiently,  a  pair  of 
substances  must  be  utilized  which  have  such  a  strong  attraction  for  one 
another  that  on  contact  one  will  immediately  absorb  the  other.  Two  such 
refrigerant-absorption  pairs  have  achieved  widespread  commercial  use: 
Ammonia-Water  and  Water-Lithium  Bromide.  Also,  efficient  system  operation 
requires  that  a  relatively  high  temperature  source  must  be  available  at  the 
generator  while  a  method  of  external  cooling  is  utilized  at  both  the 
absorber  and  the  condenser. 

Capacity  control  for  an  absorption  refrigeration  system  may  be 
accomplished  by  varying  the  energy  source  input  at  the  generator  or  by 
restricting  the  flow  rate  of  the  solution  between  the  absorber  and  the 
generator.  By  either  method  the  amount  of  refrigerant  boiled  off  in  the 
generator  is  reduced,  resulting  in  a  stronger  solution  in  the  absorber. 
This  solution  absorbs  refrigerant  at  a  slower  rate  than  before  and  an 
increase  in  the  pressure-temperature  level  of  the  evaporator  results. 
Finally,  the  higher  temperature  level  of  the  evaporator  slows  the  rate  of 
energy  transfer  from  the  desired  medium  to  the  level  needed  to  match  the 
desired  load. 


2.2  Vapor -Compression  Refrigeration 

The  primary  function  of  all  refrigeration  systems  is  to  transfer  heat 
from  a  relatively  low  temperature  medium  to  a  medium  at  a  higher 
temperature.  To  accomplish  this,  the  vapor-compression  refrigeration 
system  operates  on  a  cycle  based  on  the  energy  absorbed  or  rejected  by  a 
single  refrigerant  as  it  shifts  between  its  liquid  and  vapor  phases.  By 
controlling  the  pressures  (and  therefore  temperatures)  at  which  the  phase 
changes  occur,  it  becomes  possible  to  transfer  the  energy  between  the  two 
mediums,  with  the  requirement  that  some  additional  energy  in  the  form  of 
work  also  be  supplied  to  the  system. 

A  schematic  of  a  typical  vapor-compression  system  is  presented  in 
Figure  2.3.  This  system  is  composed  of  four  primary  components: 
compressor,  condenser,  expansion  value,  and  evaporator,  arranged  to  provide 
a  high  pressure-temperature  side  (condenser)  where  energy  is  rejected  from 
the  system  and  a  low  pressure-temperature  side  (evaporator)  where  energy  is 
absorbed.  The  remaining  components  (compressor  and  expansion  valve)  are 
utilized  to  provide  for  the  transition  of  the  refrigerant  between  the  two 
sides  of  the  system  and  also  to  regulate  the  flow  rate  of  the  refrigerant 
between  the  other  components.  A  description  of  the  operation  cycle  of  a 
vapor-compression  refrigeration  system  follows. 

2 . 2a  Vapor-Compression  Refrigeration  Cycle 

In  the  vapor-compression  cycle,  refrigerant  vapor  initially  at  a  low 
pressure-temperature  level  is  compressed  to  a  higher  level  and  flows  to  the 
condenser.  At  the  condenser,  the  high  temperature  level  of  the  refrigerant 
allows  energy  to  be  transferred  out,  resulting  in  the  condensation  of  the 
refrigerant  at  the  elevated  pressure.  From  here,  the  liquid  refrigerant 
passes  through  a  pressure  reduction  (expansion)  valve  and  flows  into  the 
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Figure  2.3   Schematic  for  Vapor-Compression  Refrigeration  systems. 
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Figure  2.4   Ideal  Vapor-Compression  Refrigeration  Cycle. 
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ENTHALPY 
Figure  2.5  a 
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Figure  2.5  b 


STEPS  .         PROCESS 

1  ...  2  Pressure  drop-Temperature  rise  in  suction  line. 

2-3  Compression  of  the  refrigerant  vapor. 
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5-6  Condensation  of  the  refrigerant. 

6-7  Subcooling  of  the  refrigerant  liquid. 
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8-9  Evaporation  of  the  refrigerant. 

9  -  1  Superheating  of  the  refrigerant  vapor. 


Figure  2.5    Actual  Vapor-Compression  Refrigeration  Cycle. 
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evaporator.  The  relatively  low  temperature  level  of  the  refrigerant  caused 
by  the  pressure  reduction  allows  the  refrigerant  to  absorb  heat  from  the 
desired  medium,  resulting  in  its  revaporization.  Finally,  the  refrigerant 
vapor  flows  back  to  the  compressor  to  repeat  the  cycle  again. 

The  above  mentioned  cycle  is  illustrated  in  its  idealized  form  on 
Pressure-Enthalpy  (Figure  2.4a)  and  Temperature-Entropy  diagrams  (Figure 
2.4b). 

The  actual  refrigeration  cycle  deviates  from  the  ideal  cycle  primarily 
because  of  additional  heat  transfer  to  or  from  the  surroundings  and 
pressure  drops  associated  with  fluid  flow.  Also,  vapor  entering  the 
compressor  will  likely  be  in  a  superheated  condition  to  eliminate  the 
possibility  of  damaging  the  machine  by  attempting  to  compress  a  liquid. 
Refrigerant  leaving  the  condenser  is  usually  subcooled  below  its  saturation 
point  to  prevent  the  flashing  of  some  of  the  refrigerant  to  its  vapor  phase 
before  reaching  the  expansion  valve,  which  would  severely  limit  the  flow 
capacity  of  the  valve  due  to  the  high  specific  volume  of  vapor  relative  to 
liquid.  Thus,  the  actual  refrigeration  cycle  would  probably  approach  those 
illustrated  in  Figures  2.5a  and  2.5b. 

2.2b  Vapor-Compression  Refriaeration/Dehumidification  Systems 

Commercial  vapor-compression  refrigeration  systems  utilized  for  the 
dehumidification  of  air  generally  fall  into  one  of  two  categories:  DX-coil 
systems,  and  chilled-water  systems.  Chilled-water  systems  are  generally 
favored  over  DX-coil  systems  for  large,  complex  multi-coil  operations  while 
DX-coil  systems  are  used  extensively  for  the  smaller,  single  coil 
operations.  For  either  type  of  system,  the  basic  schematic  illustrated  in 
Figure  2.3  is  appropriate. 
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DX-Coil  Systems 

With  DX-coil  systems,  the  air  stream  to  be  processed  is  cooled 
directly  by  contact  with  the  evaporator  (Figure  2.6),  which  is  maintained 
at  a  temperature  lower  than  the  final  desired  temperature  of  the  air. 
Thus,  the  energy  removed  from  the  air  stream  is  utilized  directly  to 
evaporate  the  refrigerant. 

Chilled-Water  Systems 

With  chilled-water  systems,  the  low  temperature  of  the  evaporator  is 
utilized  to  cool  a  supply  of  water.  This  water  is  then  circulated  to 
another  coil(s)  to  absorb  the  energy  necessary  for  the 
cooling/dehumidification  of  the  air  stream  (Figure  2.7).  Thus,  the  energy 
removed  from  the  air  stream  is  first  utilized  to  heat  the  water  supply 
which  is  then  used  to  cause  the  vaporization  of  the  refrigerant  in  the 
evaporator. 
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CHAPTER  III 
VAPOR-COMPRESSION  SYSTEM  MODELING  -  FULL  LOAD  OPERATION 

In  determining  the  energy  requirements  of  a  refrigeration  system,  the 
type  of  model  utilized  is  dependent  on  how  the  system  loading  is  varied. 
If  the  loading  is  dynamic  in  nature,  then  the  equations  which  describe  the 
interaction  of  the  various  components  will  generally  include  differential 
equations.  If,  however,  the  loading  is  under  relatively  steady-state 
conditions,  then  the  equations  describing  system  interaction  will  be  more 
algebraic  in  nature. 

The  objective  of  this  study  is  to  develop  a  model  which  may  be 
utilized  to  predict  the  overall  energy  requirements  of  reciprocating  vapor- 
compression  refrigeration/dehumidification  systems  when  subjected  to 
variable  loading  conditions.  In  general,  the  operational  dynamic  time 
constants  of  a  dehumidification  system  are  quite  short  compared  to  those  of 
typical  variable  dehumidification  loads.  Thus,  the  assumption  of  steady- 
state  loading  is  valid  provided  the  variable  load  is  broken  down  into 
sections  of  time  spent  at  each  loading  level  for  the  model  to  consider 
individually.  This  was  the  approach  utilized  for  the  modeling  of  the 
refrigeration/dehumidification  systems  in  this  study. 

3.1  System  Modeling  -  Full  Load  Operation 

As  discussed  earlier,  the  function  of  the  refrigeration  system  is  to 
transfer  energy  from  a  relatively  low  temperature  medium  to  another, 
usually  higher  temperature  medium.  The  basic  system  to  be  modeled,  as 
illustrated  in  Figure  3.1,   is  composed  of  the  following  components: 
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Figure  3.1       Basic  Vapor-Compression  Refrigeration  system. 
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reciprocating  compressor,  condenser,  evaporator,  and  expansion  valve.  A 
refrigerant,  such  as  R-22,  is  circulated  between  the  various  components, 
providing  the  means  by  which  heat  transfer  may  occur. 

With  steady-state  system  loading,  two  different  modeling  approaches 
may  be  utilized  to  predict  overall  performance.  The  first,  as  presented 
in  ASHRAE's  "Procedures  for  Simulating  the  Performance  of  Components  and 
Systems  for  Energy  Calculations"  [2],  involves  the  use  of  energy  balance 
equations  for  each  system  component.  These  equations  are  coupled  together 
by  common  variables,  such  as  refrigerant  temperatures  and  mass  flow  rates. 
This  coupling,  in  turn,  allows  the  equations  to  be  manipulated  (e.g.,  by 
use  of  Successive  Substitution  or  Newton-Raphson  iteration)  to  determine 
system  performance  for  a  given  set  of  conditions,  as  shown  by  Stoecker  [3], 
[4],  The  use  of  this  approach  in  system  modeling  provides  excellent 
agreement  with  manufacturers  equipment  catalog  data  (usually  within  +1- 
1/2%)  but  requires  the  recalculation  of  the  constants  in  the  energy  balance 
equations  for  each  change  in  equipment  size  (e.g.  constants  for  a  15  ton 
system  will  not  be  appropriate  for  a  30  ton  system). 

The  second  approach,  as  developed  by  Hittle  [5],  Leverenz  and  Bergan 
[6],  involves  the  use  of  an  input/output  model  to  directly  relate 
performance  to  the  input  conditions  for  a  subsystem  (condensing  unit, 
chiller  unit,  etc.)  of  the  overall  system.  The  equations  used  in  the  model 
are  based  on  the  overall  performance  data  of  the  subsystem,  eliminating  the 
need  to  solve  the  energy  balance  equations  for  each  of  the  subsystems 
individual  components.  Thus,  the  number  of  variables  in  the  overall  model 
is  reduced,  simplifying  the  iteration  process.  Also,  the  generalized 
nature  of  the  model  allows  for  the  use  of  the  same  equation  constants  for  a 
wide  range  of  equipment  sizes  while  still  providing  performance  predictions 
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within  +  3-1/2%  of  manufacturers  catalog  data.  However,  some  of  the 
components  of  a  system,  such  as  Direct-Expansion  evaporator  coils,  can  not 
be  categorized  in  a  generalized  subsystem  model  because  of  the  wide 
variation  in  model  size  which  may  be  placed  in  a  particular  system  by  the 
designer.  Therefore,  an  overall  system  model  based  on  generalized 
subsystems  input/output  models  may  still,  depending  upon  the  system,  also 
require  one  or  more  of  its  components  to  be  modeled  individually.  For  this 
reason,  both  component  and  the  subsystem  modeling  procedures  will  be 
presented  in  the  following  sections. 

A  summary  of  the  equations  for  the  component,  subsystem,  and  Leverenz- 
Bergan  models  is  presented  in  Appendix  A. 

3.2  Component  Modeling  -  Fuji  Load  Operation 

Reciprocating  Compressor 

The  primary  function  of  the  compressor  is  to  increase  the  pressure  and 
temperature  of  the  refrigerant  gas  from  the  level  at  which  it  vaporized  to 
a  level  which  would  give  the  gas  the  ability  to  condense  back  to  a  liquid 
(Figure  2.4  —  Step  2-3).  To  accomplish  this,  a  reciprocating  compressor 
utilizes  one  or  more  pistons  moving  in  cylinders  equipped  with  suction 
(inlet)  and  discharge  (outlet)  valves  which  allows  the  refrigerant  gas  to 
be  pumped  to  the  required  pressure  and  temperature  levels.  The  amount  of 
refrigerant  gas  which  the  compressor  handles  under  these  varying  operating 
conditions  plays  a  fundamental  role  in  both  the  overall  load  capacity  of 
the  system  and  the  energy  input  requirements  of  the  system. 

The  primary  concerns  in  modeling  a  compressor  involve  the  describing 
of  the  mass  flow  rate  which  the  compressor  can  supply  and  also  the  power 
input  for  the  compressor  in  providing  that  flow  rate.  The  compressor  mass 
flow  rate  '  is  a  function  of  the  specific  volumes  of  the  gas  at  both  the 
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suction  and  discharge  sides  of  the  compressor,  ^lso,  as  the  pressure 
difference  between  the  two  sides  increases,  the  mass  flow  rate  decreases. 
This  is  due  to  the  residual  gas  left  in  the  cylinder  after  discharge  which 
must  first  be  expanded  back  to  the  suction  pressure  level  before  more  gas 
can  be  drawn  into  the  cylinder. 

To  make  the  modeling  process  more  convenient,  several  simplifications 
are  available.  Instead  of  describing  compressor  operation  in  terms  of 
specific  volumes  and  pressures,  the  operational  parameters  are  often 
described  in  terms  of  saturation  temperature  (for  which  a  unique  saturation 
pressure  and  specific  volume  would  exist  for  the  saturated  refrigerant 
vapor).  Also,  the  compressor  mass  flow  rate  may  be  combined  with  the 
enthalpy  rise  of  the  refrigerant  in  the  evaporator  into  a  term  called  the 
compressor  capacity.  Thus,  compressor  operation  may  be  described  in  terms 
of  compressor  capacity  (Q  )  and  power  requirements  (P),  using  saturation 
temperature  as  the  variables  (T  ,T  ). 

Qe  -  f(Te,Tc) 
P-f(Tt.Te) 

where   Q  »  Compressor  Capacity 
P  ■  Compressor  Power 

T  =  Saturated  Suction  (evaporator)  Temperature 
T  =  Saturated  Discharge  (condenser)  Temperature 

The  following  equations  are  typical  of  those  used  in  refrigerant  compressor 

modeling  [2]. 

Qe  *  al  +  a2Te  +  a3Te2  +  Vc  +  H1  c     +  VeTc  +  VA 

+  VeTc2  +  VeV  t3'1' 
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P  =  bl  +  b2Te  +  b3Te2  +  b4Tc  +  b5TC2  +  Wc 

♦  b7Te2Tc  +  b8TeTc2  +  bgTe2Tc2  (3.2) 

Condenser 

The  primary  function  of  the  condenser  is  to  transfer  energy  from  the 
refrigerant  gas  to  another  medium,  thereby  allowing  the  refrigerant  to 
undergo  a  phase  change  from  a  vapor  to  a  liquid  state.  The  energy  rejected 
includes  the  compressor  capacity  (which  is  also  the  energy  absorbed  to  the 
evaporator)  and  the  mechanical  work  of  compression,  as  given  below. 

Q  =  Q  +  P  (3.3) 

c    e 

where  Q  =  heat  rejected  at  the  condenser 
Q  =  compressor  capacity 
P  =  compressor  power 

The  heat  rejection  in  the  condenser,  as  illustrated  in  Figure  2.5 
(points  4  through  7),  approaches  that  of  a  constant  pressure  process. 
Also,  the  assumption  of  a  single  temperature  for  the  refrigerant  through 
most  of  the  condenser  may  be  justified  by  the  following  two  reasons:  1) 
the  majority  of  the  energy  transfer  occurs  because  of  the  phase  change  of 
the  refrigerant,  and  2)  the  increased  temperature  difference  due  to  the 
initial  desuperheating  of  the  refrigerant  is  offset  by  lower  heat  transfer 
coefficients  [4].  Because  the  overall  heat  transfer  coefficients  do  not 
change  appreciably,  a  constant  heat-exchanger  effectiveness  may  be  assumed 
for  the  condenser.  Thus,  the  heat  transfer  at  the  condenser  (Qc)  may  be 
described  in  terms  of  the  condensing  temperature  (T  )  and  the  conditions  of 
the  medium  (air  or  water)  to  which  the  energy  is  being  transferred,  as 
shown  below. 
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Q  =  f(T  ,  T  ,  m  ) 

where  Q  =  Condenser  Load 
c 

T  =  Condensing  Temperatures 

T  =  Entering  medium  temperature 
m 

m  =  Medium  mass  flow  rate 
m 

In  general  form,   the  rate  of  heat  transfer  at  the  condenser  may  be  modeled 
by  the  following  equation, 

Q  .  F,  +  F,  (T  -  T)  (3.4) 

xc    1    2  v  c    nr 

where  F.   and  F_  are  constants  denoting  the  physical  capacity  of  the 

condenser  for  a  preset  air  or  water  mass  flow  rate. 

Another  condenser  characteristic  to  be  included  in  the  model  is  the 

amount  of  refrigerant  subcooling  which  the  condenser  is  capable  of 

supplying.   Each  degree  of  refrigerant  subcooling  increases  compressor 

capacity  by  approximately  1/2%  by  reducing  the  percentage  of  the  liquid 

flashed  during  expansion,   thereby  providing  refrigerant  to  the  evaporator 

in  a  condition  somewhat  closer  to  that  of  a  saturated  liquid  [7]   (Figure 

2.4  —  points  3  to  4  and  Figure  2.5  —  points  6  through  8).  The  percentage 

increase  in  compressor  capacity  is  generally  presented  in  manufacturers' 

catalog  data  as  a  function  of  the  temperature  difference  between  the 

refrigerant  and  the  entering  air  or  water  temperature,   and  is  represented 

by  the  equation  below. 

subcooling  =  F3  +  F4  (T.  -  TJ  +  F5(T.  -  TJ2  +  ...    (3.5) 

where  subcooling  =  %  increase  in  compressor  capacity  due  to  subcooling 

F_,  F, ,  F.  =  Constants 
3   4'   5 

For  some  applications,   it  may  also  be  necessary  to  know  the 
temperature  of  the  medium  (air  or  water)  leaving  the  condenser.   This  may 
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be  found  from  the  energy  balance  equation  for  the  medium,  as  given  below. 

Qc  "  K   Cp(Tm,out  "  Tm,in>  <3'6> 

where    C  =  specific  heat 

T_   =  leaving  medium  temperature 

T  .  =  entering  medium  temperature 

Evaporators 

The  primary  function  of  the  evaporator  is  to  transfer  energy  from  the 
desired  medium  (e.g.  air  or  water)  into  the  refrigerant,  thereby  allowing 
the  refrigerant  to  undergo  a  phase  change  from  a  liquid  to  a  vapor  state. 

As  with  the  condenser,  the  energy  transfer  in  the  evaporator 
approaches  that  of  a  constant  refrigerant  pressure-temperature  process. 
However,  unlike  the  condenser,  the  assumption  of  a  constant  heat-exchanger 
effectiveness  is  not  valid  due  to  variations  in  the  overall  heat  transfer 
coefficient  (which  are  caused  by  changes  in  the  boiling  heat  transfer 
coefficient  with  variations  in  evaporator  load)  [4].  Instead,  a  higher 
order  polynomial  must  be  used  for  the  evaporator  as  compared  to  the 
condenser.  Thus,  heat  transfer  at  the  evaporator  may  be  described  in  terms 
of  saturated  evaporating  temperature  (T  )  and  the  conditions  of  the  medium 
(air  or  water)  from  which  the  energy  is  being  transferred,  as  shown  below. 

Q  =  f(T  ,  m  ,  h  ) 
ve     e   in   m' 

where      Q        ■  Evaporator  load 

T         =  Saturated  Evaporating  Temperature 

m        ■  Medium  mass  flow  rate 
m 

h         ■  Entering  medium  enthatlpy 
For      modeling     purposes,       the     evaporators     are     divided     into       two 
categories:     DX-coil     and       Chilled  water. 
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DX-coil.  Direct  Expansion  Coils,  as  discussed  in  Chapter  2,  utilize 
direct  contact  between  the  air  stream  and  the  evaporator  to  provide  the 
cooling/dehumidification  desired  for  the  air.  Because  the  energy  transfer 
from  the  air  stream  involves  both  latent  and  sensible  energies,  the  heat 
transfer  rate  should  be  described  in  terms  of  entering  air  enthalpy  and  not 
dry-bulb  temperature.  For  modeling  convenience,  entering  air  wet-bulb 
temperature  is  generally  substituted  for  the  entering  air  enthalpy  in  the 
energy  balance  equations.  Thus,  DX-coil  evaporator  operation  may  be 
described  by  the  equations  [2] 

Qe  =  cmnne,  ♦  c2ie  +  c3xe2  +  c4iwb  ♦  c5iwb2  +  c5Teiwb  *  c7Te2Twb 

C8TeTwb2  +  C9Te2Twb2l  (3'7) 

CF  =  C1Q  +  Cu  fpm  +  C12  fpm2  (3.8) 

where  Q  =  Evaporator  load 

T  =  Saturated  evaporating  temperature 
T  ,  =  Entering  air  wet-bulb  temperature 
AREA  =  Face  area  of  coil 

CF  =  Correction  factor  for  air  flow  rate 
fpm  =  Air  velocity  to  coil 
C-pC.  j  =  Constants 

Chilled-water.  Water  chilling  evaporators,  as  discussed  in  Chapter  2, 
are  used  to  cool  a  water  supply  which  is  circulated  to  other  coils  to 
absorb  the  cooling/dehumidification  load.  Chilled  water  evaporator 
performance  may  be  described  by  the  equations  [2]. 

Qe  =  4,  +  d2(GPM)  +  d3(GPM)2  +  d4(dT)  +  d5(dT)2  +  d6(GPM)(dT) 

+  dy(GPM)2(dT)  +  d8(GPM)(dT)2  +  dg(GpM)2(dT)2        (3>g) 


23 

dT  =  (inlet  water  temperature )-( saturated  evaporating  temperature)  (3.10) 
where  Q  =  Evaporator  load 

GPM  =  Water  flow  rate 
d.— dq  =  constants 

Expansion  Valve 

The  function  of  the  expansion  valve  is  to  throttle  the  flow  of 
refrigerant  to  the  evaporator  and  decrease  the  pressure  and  temperature  of 
the  liquid  refrigerant  to  a  level  at  which  it  will  have  the  ability  to 
vaporize  easily.  This  is  illustrated  in  Figure  2.5  as  process  7-8.  For 
this  throttling  process,  very  little  heat  transfer  occurs  at  the  expansion 
device  compared  to  the  rest  of  the  system.  Therefore,  an  equation 
describing  the  process  is  not  necessary,  provided  that  the  valve  was 
properly  sized  to  provide  the  evaporator  with  a  good  supply  of  liquid 
refrigerant  for  a  wide  range  of  operating  conditions  [4]. 

Miscellaneous 

Refrigerant  Piping  Pressure  Drops.  With  the  flow  of  refrigerant  in 
pipes  between  the  primary  components  of  a  refrigeration  system,  it  is 
inevitable  that  some  pressure  drop  will  occur  in  the  refrigerant.  This 
process  is  illustrated  in  Figure  2.5  as  process  1-2  and  3^4.  The  effect  of 
these  pressure  drops  is  for  the  compressor  to  operate  with  a  suction 
pressure  lower  than  that  at  the  evaporator  exit  and  a  discharge  pressure 
higher  than  that  at  the  condenser  inlet.  Because  the  compressor  is  modeled 
in  terms  of  saturation  temperatures  (and  therefore  saturation  pressures), 
each  pressure  drop  can  be  expressed  in  terms  of  an  Equivalent  Temperature 
Change  [8].  Thus,  the  compressor  suction  temperature  would  be  the 
evaporator  temperature  minus  the  equivalent  temperature  change  of  the 
suction   line  and  the  discharge  temperature  would  be  the  condenser 
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temperature  plus  the  equivalent  temperature  change  for  the  discharge  line. 

The  actual  magnitude  of  the  pressure  drop  depends  on  the  diameter  and 
length  of  the  piping  and  the  flow  rate  of  the  refrigerant.  Therefore,  the 
effect  of  piping  pressure  drops  on  system  performance  may  vary  widely  from 
one  application  to  the  next  for  a  particular  system  capacity  size. 
Generally,  systems  with  components  spaced  closely  together  would  not  appear 
to  suffer  appreciably  from  piping-  pressure  drops  and  probably  would  not 
require  the  use  of  equivalent  temperatures  in  modeling.  However, 
applications  which  require  long  refrigerant  lines  between  system  components 
must'  account  for  pressure  drops  when  predicting  system  performance  under 
various  operating  conditions. 
3.3.  Subsystem  Modeling  -  Full  Load  Operation 

Many  of  the  refrigeration  components  on  the  market  (e.g.  compressors, 
condensers,  and  evaporators)  are  available  on  an  individual  basis  or  as 
prepackaged  units  which  combine  one  or  more  of  the  units  together.  This 
packaging  of  components  into  a  subsystem  unit  allows  the  manufacturer  to 
optimize  the  various  components  around  each  compressor.  For  this  case,  it 
is  possible  to  generalize  the  overall  performance  of  such  units  by  a  set  of 
algebraic  equations  without  the  need  to  model  each  of  the  individual 
refrigeration  components,  as  performed  by  Leverenz  and  Bergan  [6]  for  a 
water  chiller.  Because  of  the  versatility  of  the  Leverenz-Bergan  Model,  it 
was  used  as  a  basis  for  the  subsystem  modeling  in  this  study. 

A  brief  summary  of  the  Leveren-Bergan  Model  is  presented  next, 
followed  by  the  subsystem  model. 
Leverenz-Bergan  Chiller  Model 

The  performance  of  a  water  chiller  system  depends  strongly  on  the 
pressure-temperature  conditions  of  the  refrigerant  at  the  compressor. 
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These  conditions,  in  turn,  are  closely  tied  to  the  water  temperatures  at 
both  the  evaporator  and  the  condenser.  Because  of  this  dependency,  system 
capacity  and  power  requirements  may  be  described  in  terms  of  leaving 
evaporator  water  temperature  (T  )  and  entering  condensing  water 
temperature  (T  J  instead  of  refrigerant  evaporating  (Te)  and  condensing 
(T  )  temperatures.  Also,  this  dependency  eliminates  the  need  to  directly 
account  for  variations  in  water  flow  rates  for  the  model  (the  secondary 
support  systems  which  determine  T"cw  and  Tcond,  such  as  cooling  coil  and 
cooling  tower  models,  are,  however,  dependent  on  the  water  flow  rates;  as 
will  be  discussed  in  Chapter  5). 

The  model  utilizes  a  set  of  quadratic  equations  to  directly  relate  the 
full  load  capacity  and  power  requirements  of  the  system  to  the  two  water 
temperatures  (T  and  T  .).  These  equations  were  developed  utilizing 
normalized  variables  based  on  a  typical  operational  condition  of  the 
system,  T  =  7.8  °C  (46  °F)  and  TCQnd  =  29.4  °C  (85  °F),  and  the  capacity 
and  power  consumed  at  these  temperatures.  Using  these  conditions,  Leverenz 
and  Bergan  showed  that  the  model  is  capable  of  predicting  the  full  load 
capacity  and  power  requirements  of  a  chiller  system  within  +  3%  of 
manufacturers  catalog  data  for  the  normal  range  of  condenser  and  evaporator 
water  temperatures  typically  encountered  in  actual  operation.  For  the 
model,  only  the  tabulated  capacity  and  power  requirements  of  the  chiller  at 
the  above  mentioned  typical  operational  condition  is  required  as  input  to 
indicate  system  size.  The  rest  of  the  modeling  constants  are  appropriate 
for  all  chillers  regardless  of  size. 

The  chiller  model  was  validated  by  Leverenz  and  Bergan  against 
performance  data  of  a  20- ton   reciprocating  chiller  operating  in  an 


A  ton  of  refrigeration  is  equal  to  12,000  BTU  per  hour  of  operation. 
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experimental  facility  designated  to  simulate  a  conventional  four-zone  HVAC 
system.  The  difference  between  the  model  and  measured  chiller  performance 
was  less  than  5.4%  for  all  operational  points.  However,  it  should  be  noted 
(for  reasons  to  be  discussed  later)  that  chiller  was  operated  only  at 
relatively  high  water  flow  rates. 
Subsystem  Model 

Four  basic  subsystem  component  groupings  were  examined: 

1.  Compressor 

2.  Compressor-water  cooled  condenser  (water  cooled  condensing  unit) 

3.  Compressor-chilled  water  evaporator  (compressor-chiller  unit) 

4.  Compressor-chilled  water  evaporator-water  cooled  condenser 
(chiller  unit). 

Manufacturers  catalog  data  for  a  subsystem  unit  is  generally  presented 
in  terms  of  refrigerant  temperatures  and/or  water  temperatures  (entering  or 
leaving),  with  the  water  flow  rates  varied,  to  achieve  a  10  F  water 
temperature  rise  (or  drop)  for  each  operating  condition.  For  this 
situation,  the  use  of  either  entering  or  leaving  water  temperatures  in  the 
subsystem  model  is  appropriate  in  reproducing  catalog  data,  as  performed  by 
Leverenz  and  Bergan.  However,  as  can  be  seen  in  Figure  3.2,  the  inlet 
water  temperature  varies  significantly  with  a  change  in  water  flow  rate 
while  the  outlet  water  temperature  remains  relatively  stable.  Because  the 
subsystem  model  (and  the  Leverenz-Bergan  model)  do  not  take  into  account 
variations  in  water  flow  rates,  leaving  water  temperatures  must  always  be 
utilized  instead  of  entering  water  temperatures,  else  substantial  errors  in 
performance  predictions  may  result  if  a  water  flow  rate  is  used  which  does 
not  provide  a  10  F  water  temperature  swing  at  full  load  conditions. 
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refrigerant 


1 


-39  gpi 


inlet 


condenser  outlet 


Plot  of  condenser  refrigerant  and  water  temperatures  at  the 

inlet  and  outlet  af  a  15-ton  water  cooled  condensing  unit 

$  evap.teep  =  40  F  ,      cond.teip  =  105  F 


Figure  3.2  :     Variation  of  water  temperature  with  flow 
rates  for  a  water  cooled  condensing  unit. 


For  each  subsystem,  a  typical  operational  condition  was  selected  for 
use  as  the  basis  for  the  normalized  variables  in  the  model.  Constants  were 
then  derived  (using  a  least  squares  fit)  for  each  subsystem  using  the 
catalog  data  for  a  15  or  20  ton  unit  (Table  3. A).  Utilizing  these 
constants,  the  performance  of  each  subsystem  was  checked  against  catalog 
data  for  the  original  unit  size,  a  30-ton  unit,  and  a  80-ton  unit  (Appendix 
B).  This  comparison  showed  all  of  the  performance  predictions  (e.g.  the 
capacity  and  power  requirements)  for  each  subsystem  to  be  within  £  3-1/2% 
of  the  catalog  valves,  regardless  of  system  size  (Appendix  B).  It  should 
be  noted,  however,  that  because  of  the  variations  which  exist  in  component 
design  between  the  different  manufacturers,  some  subsystem  units  may  not 
fall  within  the  above  mentioned  modeling  accuracy  of  3-1/2%.  For  this 
situation,  a  new  set  of  modeling  constants  may  be  required  for  the  unit(s) 
in  question.  For  most  subsystem  units,  however,  the  parameters  presented 
in  Table  3. A  are  appropriate. 
Summary 

The  use  of  a  subsystem  model  for  predicting  overall  refrigeration 
performance  reduces  the  number  of  components  requiring  modeling, 
simplifying  the  iteration  process.  Also,  the  generalized  nature  of  the 
subsystem  model  minimizes  the  amount  of  input  required  from  catalog  data  to 
describe  subsystem  size  and  operational  characteristics,  while  still 
providing  reasonable  full  load  performance  estimates.  The  use  of  the 
subsystem  model  is  the  subject  of  Chapter  5. 


Catalog  data  based  on  16°F  of  subcooling  for  those  subsystems  without 
condensers,  actual  subcooling  for  those  with  condensers.   15  F  of  super- 
heat for  all  subsystems. 
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CHAPTER  IV 
VAPOR  COMPRESSION  SYSTEM  MODELING  ^  PARTIAL  LOAD  OPERATION 

Typically,  refrigeration  systems  are  selected  for  an  installation 
based  on  the  ability  of  the  system  to  handle  a  particular  operational 
(design)  load.  This  load,  however,  generally  occurs,  on  an  annual  basis, 
as  only  a  small  percentage  of  the  system  operational  time.  For  the 
majority  of  the  time,  the  system  encounters  a  load  condition  which  is  far 
less  severe  than  the  design  load.  Thus,  the  manner  in  which  a  system 
handles  these  off-design  conditions  has  a  strong  influence  on  the  system 
annual  energy  consumption  and  must  be  included  in  the  overall  system 
model. 

The  natural  response  for  a  refrigeration  system  to  a  drop  in 
evaporator  load  is  to  decrease  the  pressure  and  temperature  of  the 
refrigerant  in  the  evaporator.  These  lower  conditions,  in  turn,  reduce  the 
compessor  capacity  and  power  requirements  until  the  new  load  is  matched. 
However,  the  lower  temperature  of  the  evaporator  results  in  over-cooling  of 
the  air  (or  water)  and  problems  such  as  coil  frost  formation  (or  chilled 
water  freeze  up)  may  occur.  Therefore,  a  method  to  control  the  capacity  of 
the  system  is  usually  incorporated  into  initial  system  design. 

4.1  Refrigeration  System  Capacity  Control 

Capacity  control  for  a  refrigeration  system  may  be  achieved  by 
modifying  the  performance  characteristics  of  one  or  more  of  the  components 
of  the  system  (e.g.  compressor,  condenser,  and/or  evaporator)  to  match  the 
varying  load.   Because  the  compressor  is  the  major  power  consuming 
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component  of  a  system,  control  methods  which  reduce  its  capacity  and 
therefore  the  power  requirements  are  by  far  the  most  widely  utilized  in 
industry.  Condenser  and  evaporator  control  methods,  on  the  other  hand,  are 
not  as  widely  utilized  due  to  their  limited  abilities  to  reduce  system 
power  requirements  with  a  reduced  load  .  For  this  reason,  only  compressor 
capacity  control  methods  were  considered  in  this  study  as  the  primary  means 
of  system  capacity  control. 

Compressor  Capacity  Control 

Compressor  capacity  control  may  be  obtained  by  using  one  of  the 
following  methods:  cycling,  hot  gas  bypass,  back-pressure  regulator,  and 
cylinder  unloading.  All  of  these  methods,  in  effect,  reduce  the  amount  of 
refrigerant  vapor  supplied  by  the  compressor  to  the  condenser.  This 
reduced  supply,  in  turn,  results  in  less  liquid  refrigerant  at  the 
evaporator  and  a  lower  system  capacity.  A  brief  explanation  of  each  of  the 
methods  follows. 

With  cycling,  the  compressor  is  started  and  stopped  as  needed  to 
provide  the  necessary  refrigerant  flow  to  match  the  evaporator  load. 
Because  the  compressor  is  not  operating  continuously,  substantial  power 
consumption  savings  are  possible.  However,  frequent  cycling  of  a  system 
shortens  the  life  of  the  compressor,  motor,  and  starter  system.  For  this 
reason,  cycling  is  generally  used  only  with  small  durable  compressors  and 
also  with  systems  when  frequent  cycling  of  a  particular  compressor  is  not 
likely  to  occur. 


Proper  expansion  valve  operation  requires  a  minimum  pressure  differential 
across  the  valve.  To  insure  this,  condensing  temperature  is  generally  not 
allowed  to  fall  below  a  preset  level.  This  is  usually  accomplished  by 
condenser  fan  cycling  or  dampers  and  also  by  not  allowing  the  entering 
condenser  water  to  fall  below  a  minimum  temperature. 
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With  hot  gas  bypass,   refrigerant  from  the  compressor  bypasses  the 

condenser  through  a  pressure  reduction  valve  which  feeds  the  evaporator  (or 

the  compressor  suction  line).   This  method  generally  provides  for  precise 

capacity  control.   However,   the  refrigerant  which  the  compressor  must 

process  remains  close  to  that  of  full  loading,   resulting  in  near  full  load 

power  consumption  regardless  of  the  actual  evaporator  load. 

With  a  back-pressure  regulator,  the  flow  of  refrigerant  between  the 
evaporator  and  the  compressor  is  throttled,  resulting  in  the  compressor 
operating  at  a  lower  suction  pressure-temperature  level  than  the 
evaporator.  The  lower  suction  pressure-temperature  level  causes  a 
reduction  in  the  mass  flow  rate  to  the  compressor  and,  therefore,  a  reduced 
power  requirement.  The  lower  suction  level,  however,  does  result  in  a  lower 
coefficient  of  performance  (COP)  for  the  system. 

The  most  widely  utilized  method  of  compressor  capacity  control  is  the 
unloading  of  one  or  more  cylinders  in  a  multi-cylinder  compressor  by 
holding  the  suction  valve  of  the  cylinder(s)  open  throughout  the  entire 
motion  of  the  piston(s).  With  the  unloading  of  a  cylinder,  a  step 
reduction  in  compressor  refrigerant  mass  flow  rate  occurs,  resulting  in 
lower  compressor  capacity  and  power  requirements.  While  some  losses  occur 
because  of  friction  and  the  pumping  of  refrigerant  into  and  out  of  the 
unused  cylinders,  the  resulting  COP  of  a  cylinder  unloading  system  is  still 
greater  than  that  of  a  back-pressure  regulated  system,  as  shown  by  Garland 
[9].  For  the  above  reasons,  cylinder  unloading  was  chosen  as  the  primary 
means  of  system  capacity  control  for  this  study. 
4.2  Cylinder  Unloading  ^  Component  Model 

As  discussed  in  section  4,1,  a  step  reduction  occurs  in  both  capacity 
and  power  requirements  with  the  unloading  of  each  cylinder.   With  system 
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partial  load  operation,  it  is  unlikely  for  the  load  to  perfectly  match  the 
step  reduction  in  capacity  available  with  the  unloading  of  a  cylinder. 
Instead,  the  compressor  would  be  required  to  fluctuate  between  an  unloading 
step  of  under  capacity  and  one  of  over  capacity  to  enable  it  to  match  the 
total  overall  load.  For  this  situation,  Hittle  [5]  proposed  the  use  of  a 
quadratic  equation  fit  to  the  unloading  step  data  to  estimate  compressor 
performance  (equation  4.1).  Also,  a  linear  equation  may  be  used  to 
represent  the  condition  of  compressor  cycling  between  the  last  unloading 
position  and  the  off-position  (equation  4.2).  A  plot  of  these  equations  is 
presented  in  Figure  4.1. 

FFL  =  C1  +  C2(PLR)  +  C3(PLR)2  (4.1) 

FFL  =  C4(PLR)  (4.2) 

where  C.  =  0.1456 
C2  =  .9555 
C3  =  -0.1048 

C,  =  Constant  dependent  on  number  of  cylinders  allowed  to  unload 
PLR  =  Ratio  of  actual  capacity  to  the  full  load  capacity  (also  the 
ratio  of  number  of  cylinders  loaded  to  the  total  number  of 
cylinders 
FFL  =  Ratio  of  actual  power  to  the  full  load  power 
Because  of  the  normalized^ nature  of  the  variables,   the  C.   through  C_ 
constants  should  be  appropriate  for  all  compressors  while  the  C,   constant 
is  dependent  only  on  the  PLR  and  FFL  which  occurs  at  the  last  cylinder 
unloading  step  before  shut  down  (example  below). 
Example 

For  a  4  cylinder  compressor  which  can  unload  down  to  1  cylinder, 

PLR  =  1/4 


1  .  0 
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Figure  4.1       Plot  of  Fraction  of  Full   Load   (FFL)  versus 
Partial   Load  Ratio   (PLR)  for  a  four-cylinder 
compressor. 
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and  from  Equation  4.1, 

FFL  =  0.378 
therefore  C,  =     =  1.512  for  this  compressor. 


Equation  Summary 


PLR  =  (0  )/(Q         ) 

partial  load     full  load 


If  PLR  greater  than  that  of  the  last  unloading  step 
Then     FFL  =  C,  +  C2(PLR)  +  C3(PLR)2 
Else     FFL  =  C4(PLR) 


where  FFL  =  (Ppartial  load)/(Pfull  load> 

4.3  Cylinder  Unloading  -  Subsystem  Model 

The  subsystems  modeled  in  the  previous  chapter  are  all  composed  of 
refrigerant  components  grouped  around  a  compressor.  As  such,  the  capacity 
and  power  requirements  of  each  subsystem  will  decrease  with  the  unloading 
of  each  compressor  cylinder.  For  this  situation,  the  equations  outlined  in 
section  4.2  (Equations  4.1  and  4.2)  are  applicable  for  each  subsystem. 
However,  the  lower  compressor  capacity  resulting  from  the  unloading  of  a 
cylinder  requires  a  reduction  in  the  temperature  difference  between  the 
refrigerant  and  the  water  at  the  evaporator  (or  condenser).  For  a  given 
leaving  water  temperature,  this  requires  a  higher  evaporator  (or  lower 
condenser)  refrigerant  temperature.  This,  in  turn,  reduces  the  losses  in 
compressor  efficiency  for  the  subsystem  which  normally  occurs  with  cylinder 
unloading.  To  account  for  this  in  the  subsystem  model  ,  a  new  set  of 
constants  for  equation  4.1  were  developed  for  each  type  of  subsystem. 


Because  the  compressor  in  the  component  and  the  compressor  subsystem 
models  are  described  only  in  terms  of  refrigerant  temperatures,  the 
constants  developed  by  Hittle  are  adequate  for  these  two  situations. 
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These  constants  (Table  4.D)  and  the  associated  development  procedures  are 
presented  on  the  following  pages. 

4 . 3a  Subsystem  Cylinder  Unloading  Constants 

In  developing  a  set  of  constants  for  each  type  of  subsystem,  a  method 

of  determining  net  subsystem  efficiency  at  partial  loading  was  required. 

To  do  this,  each  subsystem  was  described  using  both  the  subsystem  and  the 
component  models.  Then,  the  full  load  output  of  both  models  and  the  partial 

load  output  of  the  component  model  (using  Hittle's  constants  for  Equation 

4.1)  were  calculated  for  a  variety    of   positions  through   the  normal 

operational  temperature  range  of  the  subsystem.  From  this,  the  PLR's   and 

FFL's  required  by  each  subsystem  model  to   generate  the  partial    load 

output  of  the  component  models  was  found.   These  data  were  then  used  to 

generate  the  constants  required  by  each  subsystem  for  use  in  the  fraction 

of  full  load  equation  (Equation  4.1). 

Component  Model 

The  subsystems  to  be  modeled  are: 

1)  compressor-water  cooled  condenser  (water-cooled  condensing  unit) 

2)  compressor  -  chilled  water  evaporator  (compressor-chiller  unit) 

3)  compressor-chilled  water  evaporator-water  cooled  condenser 
(chiller  unit) 

The  subsystems  modeled  were  Trane  15-ton  units  (RWUB  C15M,  CCAB  C15M,  and 

CGWB  C15M)  which,  like  most  available,  use  the  same  compressor  (and  also 

condenser  and  evaporator  where  applicable)  for  all  three  units.  Because  of 

this,  the  problem  of  breaking  such  subsystems  down  into  individual 

components  was  greatly  simplified.   The  performance  data  of  each  component 


15-ton  subsystem  units  from  Trane  Air  Conditioning  were  utilized  in  this 
process. 
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COMPRESSOR 

8e=sub'-dt(Al*A2tTe+A3UeA2*A4Uc»A5»TcA2tA6neITctA7tTeA2»Tc*A8«TetTcA2*A9tTeA2tTcA2) 
P=Bl*82tTe+B3»TeA2m<T:tB5trcA2+B6JTeirc*B71TeA2*Tc*B8»TeircA2+B9tTeA2»TcA2 


»1« 

17.3676431428 

61« 

-13.61523303 

A7= 

-8.48496697674E-2 

B2= 

1.321157156 

M« 

7.79868275349E-3 

B3= 

-.02127362189 

A4= 

-.157051968372 

B4= 

.455646584 

A5= 

4.15010446512E-4 

B5= 

-.002043495701 

u> 

5.41307441B6E-3 

86= 

-.0249864533 

A7= 

-9.55775851163E-5 

B7= 

.000393825407 

Afl= 

-2.76754412093E-5 

B8= 

1.281939044E-4 

A9= 

3.79140383256E-7 

39= 

-1.834600B45E-6 

WATER-COOLED  CONDENSER 

Qc=9e+P 
Qc=ZI(DT) 
DT=Tc-Tw  in 
Z=Cl+C2t(GP!1]+C3»(6PMA2) 
subcool=Dl+OZt (BT) +D3t (DT"2) 


Cl=. 09348214359  Dl=. 98716936174 

C2=2.25321428653E-2  D2*.  00461114837 

C3=-1.4833333318E-4  D3=-3.2717021235E-5 


CHILLED-WATER  EVAPORATOR 

DT=T»  in-Te 
'=FltF2«8Pn+F3t6PNA2tF4IDTtF5tOTA2tP6«6P«5DTtF7t6P«A2»llTtF8t6PNtDT'2tF9lBP«-'2tOTA 


F1M4. 0665439 
F2=-2.07091539 
F3=. 0241712729 
F4=-4. 22431777 
F5=. 11196747 
F6=. 231325698 
F7=-. 00256334483 
F8=-. 0053363574 
F9=. 00006085379 


TABLE  4. A         15-ton  refrigeration  system 

component  equations  and  constants. 


...out..  ....in re-frig 

Tel      Tel  Tel  Tel       Te  Tc  subcl  Be.cat  Qe.pred  I 

40.0     85.0  50.0  75.0  31.3  98.5  16.5  15.10  15.08  -.12 

4O.0      90,0  50.0  30.0  31.6  103.3  16.4  14.70  14.68  -.15 

40.0      95.0  50.0  85.0  31.8  108.1  16.3  14.20  14.23  .21 

40.0  100.0  50.0  90.0  32.1  113.0  16.2  13.80  13.80  .01 

40.0  105.0  50.0  95.0  32.4  117.8  16.1  13.40  13.36  -.27 

42.0      85.0  52.0  75.0  32.9  98.3  16.8  15.40  15.41  .09 

42.0      90.0  52.0  80.0  33.2  103.7  16.7  15.20  15.19  -.07 

42.0      95.0  52.0  84.9  33.4  108.4  16.5  14.70  14.74  .24 

42.0  100.0  52.0  90.0  33.8  113.3  16.4  14.30  14.30  -.02 

42.0  105.0  52.0  95.0  34.1  118.1  14.3  13.90  13.85  -.38 

44.0      85.0  34.0  75.0  34.5  99.2  17.0  16.10  16.15  .30 

44.0      90.0  54.0  30.0  34.8  104.1  14.9  15.70  15.71  .07 

44.0     95.0  54.0  85.0  35.2  108.9  14.8  15.30  15.27  -.22 

44.0  100.0  54.0  9O.0  35.4  113.7  16.7  14.80  14.79  -.04 

44.0  105.0  54.0  95.0  35.8  118.5  16.5  14.40  14.34  -.29 

44.0      35.0  54.0  75.0  34.2  99.4  17.3  14.40  14.49  .54 

44.0     90.0  56.0  80.0  34.5  104.4  17.1  14.20  14.25  .33 

44.0     95.0  54.0  85.0  34.8  109.2  17.0  15.80  15.80  -.03 

44.0  100.0  54.0  90.0  37.1  114.1  14.9  15.30  15.33  .18 

44.0  105.0  54.0  95.0  37.4  118.8  14.8  14.90  14.34  -.25 

48.0      85.0  58.0  75.0  37.8  100.0  17.5  17.20  17.25  .29 

48.0      90.0  58.0  80.0  38.1  104.8  17.4  14.80  14.79  -.07 

48.0     95.0  58.0  35.0  38.4  109.4  17.3  14.30  14.33  .14 

48.0  100.0  58.0  90.0  38.3  114.5  17.1  15.80  15.35  .29 

48.0  105.0  58.0  95.0  39.1  119.2  17.0  15.40  15.37  -.22. 

30.0      85.0  40.0  75.0  39.4  100.4  17.8  17.70  17.79  .50 

30.0     90.0  60.0  30.0  39.7  105.2  17.4  17.30  17.32  .13 

50.0      95.0  60.0  35.0  40.0  110.0  17.5  16.30  14.34  .23 

50.0  100.0  40.0  90.0  40.4  114.8  17.4  14.40  14.37  -.14 

50.0  105.0  40.0  95.0  40.7  119.4  17.3  15.90  15.88  -.12 


WHERE 
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P.eit 

P.pred 

I 

SPMc 

SFMe 

— — 

--" 

■"■■" 

14.40 

14.41 

.04 

46.03 

36.20 

14.90 

14.89 

-.09 

45.39 

35.23 

15.30 

15.38 

.51 

44.65 

34.15 

13.80 

15.39 

.40 

43.97 

33.12 

14.40 

16.43 

.14 

43.23 

32.07 

14.40 

14.58 

-.14 

47.43 

37.47 

15.10 

15.08 

-.10 

44.75 

36.45 

15.40 

15.59 

-.07 

44.01 

35.37 

14.10 

16.13 

.17 

43.32 

34.31 

14.60 

16.68 

.49 

44.42 

33.23 

14.80 

14.76 

-.27 

48.83 

38.75 

15.30 

15.29 

-.10 

48.14 

37.71 

15.30 

15.33 

.17 

47.44 

34.64 

14.30 

16.37 

.42 

44.68 

35.51 

14.90 

16.95. 

.32 

46.04 

34.44 

15.00 

14.94 

-.39 

50.26 

40.04 

15.50 

15.48 

-.13 

49.57 

39.01 

14.00 

16.04 

.27 

48.36 

37.91 

14.60 

16.62 

.14 

48.13 

36.79 

17.10 

17.21 

.44 

47.41 

35.67 

15.10 

15.13 

.17 

51.72 

41.40 

15.70 

15.49 

-.04 

51.00 

40.29 

16.20 

16.27 

.43 

50.29 

39.13 

16.80 

14.87 

.41 

49.54 

33.03 

17.40 

17.47. 

.43 

48.31 

34.38 

15.30 

15.30 

.03 

53.14 

42.49 

15.90 

15.88 

-.10 

52.41 

41.58 

16.40 

14.49 

.54 

51.63 

40.41 

17.00 

17.12 

.49 

50.98 

39.30 

17.60 

17.74 

.80 

50.23 

33.11 

Tel  =  EVAPORATOR  BATES  TEMPERATURES 

Tel  ■  CONDENSER  WATER  TEMPERATURES 

Te  =  EVAPORATOR  REFRIGERANT  TEflP. 

Tc  =  CONDENSER  REFRIGERANT  TENP. 

SPMe  =  EVAPORATOR  WATER  FLOW  RATE 

SPNc  =  CONOENSBR  WATER  FLOW  RATE 

subel  •   DECREES  OF  SUBCCQUN6 

9e  cat  =  CATALOG  SYSTEM  CAPACITY 

dtjrtd"  PREDICTED  SYSTEM  CAPACITY 

P  cat  *  CATALOG  SYSTEM  POWER 

P'pred  =  PREDICTED  SYSTEM  POWER 


TABLE  4.B       Comparision  of  predicted  performance  versus  catalog  data 
for  the  15-ton  chiller  unit  component  model. 
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was  directly  available  for  the  compressor  (CVWB  unit)  and  the  condenser 

(RWVB  unit),  while  the  data  for  the  evaporator  was  found  by  comparing  the 
compressor  (CUWB  unit)  and  the  compressor-chilled  water  evaporator  (CCAB 
unit)  data  to  find  all  of  the  necessary  variables  (Qe,  GPM,  Tg,  Twater  in). 
Using  these  data  as  a  base,  a  least-squares  fit  was  performed  for  each 
equation.  The  resulting  equations  (Table  4. A),  predict  performance  values 
well  within  1-1/2%  of  catalog  values.  As  a  farther  check,  the  equations 
were  used  to  predict  the  performance  of  the  chiller  unit  through  its  entire 
operational  temperature  range  (Table  4.B).  Examination  of  these  results 
showed  all  of  the  predicted  capacity  values  (column  9)  to  be  within  0.56% 
(column  10)  of  the  catalog  capacity  values  (column  8).  Also,  the  predicted 
power  consumption  (column  12)  was  within  0.80%  (column  13)  of  the  catalog 
power  consumption  (column  11)  for  all  of  the  operational  temperatures. 
Thus,  these  equations  should  provide  reasonable  estimates  of  performance 
for  all  three  subsystems. 
Subsystem  Model 

The  Trane  15-ton  subsystem  units  were  originally  used  to  derive  the 
constants  for  the  subsystem  models.  Thus,  the  constants  represented  in 
Table  3. A  are  appropriate  for  the  subsystems  being  examined  in  this 
chapter.  The  base  capacity  and  power  is  presented  below  for  each 
subsystem. 


T  .  ,  D 

el 

T  .  ,  D 
cl 

CAP,  D 

FLP,  D 

water-cooled 
condensing  unit 

4.44  °C 
(40  °F) 

35  °C 
(95  °F) 

59.1  kW 
(16.8  tons) 

16.5  kW 

compressor- 
chiller  unit 

7.78  °C 
(46  °F) 

40.56  °C 
(105  °F) 

56.6  lcW 
(16.1  tons) 

15.6  kU 

chiller  unit 

7.78  °C 
(46  °F) 

35  °C 
(95  °F) 

55.5  kW 
(15.8  tons) 

16  kW 

Table  4.C 
15-Ton  Subsystem  Base  Capacity  and  Power  Requirements 


40 
Model  Comparison 

As  a  basis  for  comparison,  the  operational  temperature  range  of  each 
subsystem  was  assumed  to  be  the  same  as  those  utilized  in  the  initial 
calculation  of  the  subsystem  modeling  constants.  These  temperature  ranges 
were  broken  down  into  the  following  intervals: 


Water  Cooled 
Condenser  Unit 

Compressor 

el 

-Chiller  Unit 

T 
c 

Chiller 

Unit 

T 
e 

cl 
85  F 

el 

T  i 
cl 

30  F 

40  F 

85 

F 

40  F 

85  F 

40  F 

95  F 

45  F 

95 

F 

45  F 

90  F 

50  F 

105  F 

50  F 

105 

F 

50  F 

95  F 

115  F 

115 

F 

100  F 

125  F 

125 

F 

105  F 

For  each  interval,   partial  load  performance  was  calculated  utilizing 

the  component  model  at  100,  80,  60,  40,  and  20  percent  of  cylinder  loading. 

Also,   the  full  load  capacity  and  power  predictions  of  the  subsystem  models 

were  found.  From  these  data,  the  PLR's  and  FFL's  required  by  the  subsystem 

models  to  provide  the  same  partial  load  performance  values  as  the  component 

models  were  calculated  using  the  following  equations: 

PLR  .       =  (component  model  partial  load  capacity)/ 
subsystem 

(subsystem  model  full  load  capacity) 

FFL  ,       =  (component  model  partial  load  power)/ 
subsystem 

(subsystem  model  full  load  power) 
New  constants  for  equation  4.1  were  derived  from  .these  data  for  each  of  the 
subsystems  and  are  presented  in  Table  4.D.   Also,  plots  of  the  PLR  and  FFL 
data  points  are  illustrated  in  Figures  4.2  through  4.8. 
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cl  = 

Compressor 

Wate 
Cond 

r  Cooled 
ensing  Unit 

Compressor- 
Chiller  Unit 

Chiller 
Unit 

0.1456 

0.129 

0.1545 

0.1334 

c2  = 

0.9555 

0.8849 

0.8323 

0.7439 

c3  = 

-0.1048 

-0.01302 

0.01672 

0.1209 

Table  4.D 
Subsystem  Model  Cylinder  Unloading  Constants  (Equation  4.1) 
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PLR 


1   ,  0 


Figure  4.2      Plot  of  Fraction  of  Full   Load   (FFL)  versus 

Partial   Load  Ratio   (PLR)  for  the  Water  Cooled 
Condensing  Unit  of  Section  4.3. 
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Figure  4.3       Least  Squares  Fit  of  FFL  versus  PLR  data 
for  the  Water  Cooled  Condensing  Unit 
of  Section  4.3. 


43 


1  .  0 


3  - 


FFL 

.  4 


8  .  0  • » 


i  o  m  p  r 
i  h  i  1  1 


e  r   unit 


J . I , I , I , I , I , I ,_!_ 


PLR 


1  ,  0 


Figure  4.4   Plot  of  Fraction  of  Full  Load  (FFL)  versus 
Partial  Load  Ratio  (PLR)  for  the  Compressor- 
Chiller  Unit  of  Section  4.3. 
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Figure  4.5   Least  Squares  Fit  of  FFL  versus  PLR  data 
for  the  Compressor-Chiller  Unit 
of  Section  4.3. 
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Figure  4.6   Plot  of  Fraction  of  Full  Load  (FFL)  versus 

Partial  Load  Ratio  (PLR)  for  the  Chiller  Unit 
of  Section  4.3. 
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Figure  4.7   Least  Squares  Fit  of  FFL  versus  PLR  data 
for  the  Chiller  Unit  of  Section  4.3. 
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Figure  4.8       Plot     of  Partial   Load  Performance   (FFL  vs.  PLR) 
for  a  Compressor, Compressor-Chi Her  Unit, 
Water  Cooled  Condensing  Unit,  and  a  Chiller  Unit. 
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CHAPTER  V 

APPLICATIONS 

5.1  Secondary  Support  Models 

The  performance  of  a  refrigeration  system  is  strongly  influenced  by 
the  operational  characteristics  of  its  support  systems  (cooling  towers, 
cooling  coil  controls,  etc.).  Because  of  this,  models  of  these  support 
systems  must  be  included  in  the  overall  model  used  to  predict  refrigeration 
system  performance. 

5.1A  Cooling  Tower  Model 

When  recirculated  water  is  used  in  liquid-cooled  condensers,  it  is 
generally  cooled  by  being  routed  through  a  cooling  tower.  At  the  cooling 
tower,  some  of  the  water  is  allowed  to  evaporate  into  the  surrounding  air. 
The  energy  required  for  evaporation  is  drawn  from  the  unevaporated  water, 
resulting  in  a  drop  in  the  temperature  of  the  remaining  liquid  water.  The 
cooled  water  is  then  returned  to  the  condenser  to  absorb  more  energy  from 
the  refrigerant  undergoing  condensation. 

The  model  for  counterflow  cooling  tower  performance  developed  by  Webb 
and  Villacres  [10]  provides  excellent  correlation  with  manufacturers  data 
and  was  used  as  the  basis  for  the  cooling  tower  model  in  this  study. 

The  Webb-Villacres  model  is  capable  of  handling  three  operational 
variations: 

1.   Fixed  air  and  water  flow  rates  with  leaving  water  temperature 
allowed  to  vary  with  load. 
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2.  Fixed  water  flow  rate  and  leaving  water  temperature  with  a 

variable  air  flow  rate. 

3.  Fixed  air  flow  rate  and  leaving  water  temperature  with  a 
variable  water  flow  rate. 

Because  most  cooling  towers  are  not  designed  for  a  significant  variation  in 
water  flow  rate  [10],  only  the  first  two  cooling  tower  operational  cases 
were  included  in  the  overall  model.  Also,  some  cooling  towers  provide  a 
fixed  leaving  water  temperature  by  fan  cycling.  To  approximate  the 
operational  time  of  these  fans  in  the  overall  model,  the  ratio  of  the 
required  temperature  drop  to  the  temperature  drop  of  first  operation  case 
was  used  (Equation  5.1). 

%  fan     =  [(T   -  T  .   .  )/(T   -  T  . )]*100       (5.1) 
n        _.   ,  »,.      L  co    ci.min'   co    ci/J 
Operational  Time 

where   T   =  leaving  condenser  water  temp, 
(entering  cooling  tower) 

T  .   .   =  minimum  entering  condenser  water  temp. 


ci,min 


(leaving  cooling  tower) 


T  .  =  leaving  cooling  tower  water  temp, 
of  operational  case  #1 


Pitfall  of  Webb-Viilacres  Model  -  The  iteration  technique  utilized 
in  the  Webb-Villacres  model  was  the  bisection  method  (where  high  and  low 
estimates  of  a  variable  are  used  to  determine  its  actual  value).  However, 
because  of  the  nature  of  the  NTU  Equation  in  the  model,  the  use  of  too  low 
of  an  initial  estimate  of  the  variable  may  result  in  the  cooling  tower 
model  providing  an  incorrect  leaving  water  temperature  or  air  flow  rate 
(Figure  5.1).  To  prevent  this  from  occurring,  the  variable  should  first  be 
evaluated  at  the  high  estimate  (T  .  equal  to  T   or  tower  air  flow  equal  to 

•  CI  CO  M 

design  air  flow)  and  then  have  the  estimate  reduced  in  small  increments 
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correct 
Tci  or  CFM 


Tci 


Tco 


Tci  -  Twb 


Figure  5.1   Typical  plot  of  NTU  versus  cooling  tower  temperature 
range  for  a  low  ambient  wet-bulb,  (normal  operation 
should  ALWAYS  be  to  the  left  of  the  discontinuity)  . 
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until  a  low  estimate  is  available.   At  this  point,  the  bisection  method  may 

be  used  to  find  the  final  value  of  the  variable.   This  approach,  while 

requiring  more  iteration  steps,   was  found  to  always  provide  the  correct 

solution  provided  small  enough  increment  steps  were  used. 

5. IB.  Cooling  Coil  Model 

Dehumidification  by  refrigeration  requires  the  removal  of  energy  from 
the  air-water  vapor  mixture,  resulting  in  a  reduction  in  air  temperature 
and  the  condensation  of  the  undesired  water  vapor  from  the  mixture. 
Cooling  coils  for  this  purpose  tend  to  be  of  multiple-row  construction 
which,  in  dehumidifying,  provide  leaving  air  conditions  at  or  very  near 
that  of  saturation  (Figure  5.2)  [11].  These  coils,  because  of  the 
relatively  large  amount  of  heat  transfer  surface  available,  are  able  to 
operate  at  higher  refrigerant  and/or  chilled  water  temperatures  than  coils 
of  less  depth  for  dehumidifying  air  to  a  given  moisture  content  level  [12]. 
Thus,  coil  frost  formaton  and/or  chilled  water  icing  are  less  likely  to 
occur  with  'deep'  coils.  Also,  the  higher  refrigerant  temperature  improves 
compressor  efficiency  and  reduces  power  consumption.  For  these  reasons, 
the  cooling  coil  used  in  this  model  was  assumed  to  be  of  a  'deep'  multiple- 
row  construction. 

As  was  stated  above,  'deep'  multi-row  cooling  coils  provide 
dehumidified  air  at  or  very  near  that  of  saturaton.  To  simplify  coil 
modeling,  the  dehumidified  air  leaving  the  coil  was  assumed  to  be  at 
saturation.  In  doing  so,  the  need  for  equations  describing  the  variation 
of  the  leaving  air  conditions  with  coil  loading  was  eliminated.  It  should 
be  noted,  however,  that  this  assumption  will  provide  for  a  slight  over 
estimation  of  the  actual  cooling  coil  load  required  to  provide  a  given 
leaving  air-moisture  content  level  (Figure  5.3). 
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AIR  DRY-BULB  TEMPERATURE 

Figure  5.2   Psychrometric  chart  for  dehumidification 
by  a  deep,  multiple-row  cooling  coil. 


AIR  DRY-BULB  TEMPERATURE 


Figure  5.3  Psychrometric  chart  for  dehumidification 
for  the  cooling  coil  model. 
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Cooling  Load 

The  load  to  be  absorbed  at  the  cooling  coil  is  given  by 


%   "  Vr^l  "  V  (5-2) 


where  h.  =  entering  air  enthalpy 

h«  =  leaving  air  enthalpy 

Q  =  coil  load 
e 


u  .   =  mass  flow  rate  of  the  air 
air 


and 


h  =  0.24(Tdb)  +  W(1061  +  0.444(Tdb))  (5.3) 


where   h  =  enthalpy  of  the  air 
T,,  =  dry  bulb  temperature 
W  =  moisture  content  of  air 
Air  enthalpy,  as  can  be  seen  in  Equation  5.3,  is  a  function  of  both  dry- 
bulb  temperature  and  the  moisture  content  of  the  air.   However,   coil 
conditions   are  frequently  known  in  terms  of  dry-bulb  and  wet-bulb 
temperatures.   For  this  situation,   air  enthalpy  may  be  estimated  from  the 
wet-bulb  temperature  utilizing  the  subroutines  developed  by  Webb  and 
Villacres  [9].   Thus,  air  enthalpy  may  be  estimated  by  either  of  the  above 
methods  for  most  of  the  circumstances  which  a  cooling  coil  would  normally 
encounter.   The  calculation  procedure  used  to  determine  coil  load  (Q  )  is 
outlined  below. 

Calculation  Procedure 


1)  Calculated  entering  air  enthalpy  by  either 

a)  Eqn.  5.3  (Tjb,  W) 

b)  Wet-bulb  (Tdb,  Twfa) 

2)  If  only  entering  dry-bulb  and  wet-bulb  specified,  calculate 
entering  moisture  content  from  Equation  5.3 
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3)  For  dehuraidification,  leaving  wet-bulb  is  approximately  equal  to 

leaving  dry-bulb,  calculate  enthalpy  based  on  this  wet-bulb 

4)  From  Eqn.  5.3,  calculate  saturated  moisture  content  of  Step  #3 

5)  If  saturated  moisture  content  is  less  than  inlet  moisture 
content,  then  leaving  air  enthalpy  is  enthalpy  of  Step  #3 

If  saturated  moisture  content  is  greater  than  inlet  moisture 
content,  then  use  inlet  moisture  content,  outlet  dry  bulb,  and 
Equation  5.3  to  calculate  leaving  air  enthalpy 

6)  Calculate  coil  load  (Q  )  using  Equation  5.2 

5.2  Generalized  Model-Refrigeration/Dehumidif ication  System 

For  the  purposes  of  this  study,  four  basic  component  variations  of  a 
refrigeration/dehumidif ication  system  were  examined  (Figures  5.4  through 
5.7): 

1.  Compressor  unit  -  DX-coil  -  air  cooled  condenser  system 

2.  Water  cooled  condensing  unit  -  DX-coil  -  cooling  tower  system 

3.  Compressor  chiller  unit  -  chilled  water  coil  -  air  cooled 
condenser  system 

4.  Chiller  unit  -  chilled  water  coil  -  cooling  tower  system 

These  systems  were  all  incorporated  into  a  generalized  computer  program. 
This  program,  which  is  listed  in  Appendix  C,  is  described  below. 

5 . 2A  System  Controls 

DX-coil  systems.  As  outlined  in  Chapter  3,  DX-coil  evaporator  load  is 
a  function  of  the  refrigerant  temperature  and  of  the  inlet  air  conditions 
(and  therefore  outlet  air  conditions).  For  a  DX-coil  to  provide  a  given 
cooling/dehumidif ication  load  (cooling  coil  model),  the  coil  must  operate 
with  a  specific  refrigerant  temperature  (as  defined  by  Eqns.   3.7,   3.8). 
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Figure  5.4       Compressor  Unit~DX  Coil—Air  Cooled  Condenser  system. 


C WATER 


WATER  COOLED 

CONDENSING 

UNIT 


?3 


.COOLING 
TOWER 


■f.h 


(REFRIGERANT 


■f  f- 


.EXPANSION 
VALVE 

DX  COIL 
EVAPORATOR 


AIR   FLOW 


Figure  5.5  Water  Cooled  Condensing  Unit— DX  coil— Cool  ing  Tower  system. 
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Figure  5.6     Compressor  Chiller  Unit—Chilled  Water  Coil- 
Air  Cooled  Condenser  system. 
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Figure  5.7       Chiller  Unit—Chilled  Water  Coil—Cooling  Tower  system. 
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55 
This  temperature,   in  turn,  varies  with  coil  loading  conditions   (e.g. 

maintaining  of  fixed  evaporator  temperature  would  result  in  the  over 
cooling/dehumidifying  of  the  air  stream  for  a  reduction  in  the  inlet  air 
wet-bulb  temperature  or  air  flow  rate).  For  a  DX-coil  system,  the  variable 
evaporator  temperature  and  system  capacity  required  for  the  desired  air 
outlet  conditions  may  be  provided  by  regulating  the  cylinder  unloading 
process  (such  as  by  a  thermostat  in  the  outlet  air  duct).  This  was  the 
procedure  used  in  this  model. 

Chilled-water  systems.  The  schematic  used  for  modeling  chilled-water 
systems  is  illustrated  in  Figure  5.8.  As  with  the  DX-coil  systems,  coil 
temperature  of  a  chilled-water  system  must  be  allowed  to  vary  with  inlet 
air  conditions  to  provide  the  desired  levels  of  cooling/dehumidifying.  To 
accomplish  this,  a  thermostat  in  the  air  stream  is  generally  used  to 
control  a  valve  at  the  coil  to  regulate  the  chilled-water  flow  rate  through 
the  coil.  Thus,  cooling  coil  chilled-water  flow  rate  will  vary  directly 
with  the  coil  load.  Most  systems,  however,  maintain  a  constant  chilled- 
water  flow  rate  at  the  evaporator.  To  allow  for  this,  a  chilled  water- 
cooling  coil  bypass  circuit  is  generally  included  in  these  systems;  with 
either  a  bypass  circuit  at  each  coil  or  a  single  bypass  circuit  located 
near  the  evaporator.  The  use  of  these  circuits  eliminates  the  need  for 
modeling  individual  chilled-water  coil  performance  (e.g.  coil  water  flow 
rates  for  each  loading).  Instead,  the  load  absorbed  at  the  coil  (from  the 
cooling  coil  model)  and  an  energy  balance  at  the  evaporator  (Equation  5.4) 
may  be  used  to  predict  evaporator  performance. 

Q  =  m  C  (T     -  T    )  (5.4) 

e    m  p  cw,r    cw,s 

where  Q  =  evaporator  load  (cooling  coil  model) 
m  =  water  flow  rate  through  evaporator 
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Figure  5.8   Chilled  water  flow  route  schematic. 
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T     =  return  mixed  water  temperature 
cw,r 

T     =  chilled  water  supply  temperature 
cw,s 

C  =  specified  heat  of  water 
P 

For  these  systems,  cylinder  unloading  is  typically  used  to  match 
system  capacity  (and  therefore  evaporator  capacity)  to  the  desired  load  by 
providing  either  a  preset  supply  or  return  mixed  chilled  water  temperature. 
Thus,  as  coil  load  increases  (or  decreases),  the  number  of  compressor 
cylinders  loaded  must  increase  (or  decrease)  and  the  water  temperature  drop 
across  the  evaporator  must  change  (as  given  by  Equation  5.4)  for  the  system 
to  maintain  the  preset  water  temperature.  For  a  system  with  return  water 
temperature  control,  the  response  to  an  increase  coil  load  is  to  provide 
colder  chilled  water.  For  a  system  which  controls  supply  water 
temperature,  the  return  water  temperature  will  increase  with  coil  loading. 

Air-cooled  condenser  systems.  As  briefly  mentioned  in  Chapter  4,  to 
maintain  proper  expansion  valve  operation,  the  condensing  temperature  is 
generally  not  allowed  to  fall  below  a  preset  level.  To  accomplish  this, 
the  air  flow  rate  across  the  condenser  is  generally  reduced  at  low  ambient 
temperatures  by  either  the  use  of  dampers  or  by  cycling  the  fans.  For  this 
study,  a  minimum  condensing  temperature  was  included  in  the  model.  A 
method  of  estimating  fan  cycling  time,  however,  is  not  currently  part  of 
the  model. 

Water-cooled  condenser  system.  As  mentioned  in  the  previous  section, 
cooling  towers  do  not  operate  efficiently  at  water  flow  rates  other  than 
those  of  design.  Because  of  this,  the  water  flow  rate  of  the  condenser  and 
the  cooling  tower  were  assumed  to  be  the  same  at  all  operating  conditions 
(e.g.  no  bypass  circuits).  However,  proper  expansion  valve  operation 
requires  that  entering  condensing  water  temperature  does  not  fall  below  a 
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preset  level  (as  would  be  the  case  for  low  ambient  wet-bulb  temperatures). 

For  this  situation,   the  leaving  cooling  tower  water  temperature  is 

maintained  at  the  minimum  entering  condensing  water  temperature  by  either 

cooling  tower  dampers  or  fan  cycling  —  as  outlined  in  section  5.1. 

5.2B.  Use  of  the  Generalized  Model 

The  flow  chart  for  the  computer  program  listed  in  Appendix  C  is 

illustrated  in  Figure  5.9.   The  use  of  this  model  may  be  divided  into  four 

steps: 

1.  Determining  system  equipment  constants. 

2.  Calculating  T  .  of  subsystem  model 

el 

3.  Iterating  for  T  ,  of  subsystem  model 

4.  Performance  output 

System  equipment  constants.  Constants  for  15-ton  systems  of  all  four 
of '  the  component  variations  handled  by  this  model  are  included  in  the 
constants  subroutine  of  the  program.  The  procedures  for  finding  the 
constants  for  the  appropriate  system  is  as  follows: 

a)  subsystem  unit  —  input  catalog  data  for  capacity  and  power  at  the 
base  temperatures  as  outlined  in  Chapter  3. 

—  input  appropriate  subsystem  constants  and 
cycling  PLE  as  outlined  in  Chapters  3  and  4. 

b)  DX-coil  —  input  constants  for  equation  fitting  of  Equations  3.7, 

3.8  (constants  given  for  6-row  Aerofin  coils) 

—  input  face  area  of  coil. 

c)  Chilled-water  coil  —  select  either  return  or  supply  water 

temperature  control 

—  input  water  temperature  chosen  for  the  above 

—  input  chilled-water  flow  rate 
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d)  Air-cooled  condensers  —  input  constants  for  equation  fit  of 

Equations  3.4,  3.5  (constants  given  for  Equation  3.5  for  all 
Trane  air-cooled  condensers). 

—  input  minimum  condensing  temperature 

e)  Cooling  towers  —  input  Webb-Villacres  cooling  tower  constants 
(program  of  Appendix  G). 

—  input  cooling  tower  design  air  and  water  flow  rates 

—  input  minimum  condensing  water  temperature 

rer 


the  desired  outlet  conditions,   the  cooling  coil  load  may  be  calculated  (Q 

—  cooling  coil  model).   For  a  DX-coil,  equation  3.7  and  3.8  may  be  solved 

directly  utilizing  the  load  and  inlet  conditions  to  find  T  .   For  a 

chilled-water  system,   equation  5.4  may  be  solved  directly  for  the  leaving 

chilled  water  temperature  (T  ).   The  appropriate  temperature  is  then  used 

as  the  T  ,  variable  in  the  subsystem  model. 

Iterating    T  ,.   Unlike  the  T  ,  variable,   the  T  ,  must  be  found 
cl  el  cl 

through  an  iteration  process.   First,   an  estimate  of  the  T  ,  variable  is 

o  i-  cl 

made.   Utilizing  this  estimate  of  T  .  and  the  known  T  , ,   the  full  load 
s  cl  el' 

capacity  and  power  of  the  system  is  calculated  (subsystem  model).   Using 

these  values  and  the  cooling  coil  load,   the  partial  load  power  is 

calculated  (from  Equations  4.1,   4.2).   The  partial  load  power  and  cooling 

coil  load  is  summed  to  find  the  condenser  load.  The  condenser  load  is  then 

used  to  find  a  new  condensing  temperature  (T  )  for  an  air-cooled  condenser 

(Eqns.   3.4,   3.5)  or  a  new  leaving  condensing  water  temperature  (T  )  for 

water-cooled  condenser  (cooling  tower  model  and  Eqn.   3.6).   These  values 

are  then  used  as  the  new  estimates  of  T  , ;  with  the  orocess  beins  repeated 

cl  '  s   * 

until  the  convergence  criteria  is  met. 
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figure  5.9   Generalized  model  flowchart 
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Performance  Output .   At  this  stage,   the  partial  load  performance  for 

the  cooling  tower  -  if  desired  -  may  be  calculated  (reduced  cfm,   fan 

cycling).   The  performance  data  for  the  system  at  this  steady-state 

condition  may  then  be  multiplied  by  the  number  of  hours  of  occurrence  of 

this  condition  for  inclusion  in  the  total  annual  performance  estimate  . 

The  process  is  then  repeated  starting  with  Step  #2  (calculating  T  .)  until 

all  of  the  likely  operational  coil  and  ambient  conditions  have  been 

accounted  for. 

5.3  Typical  Annual  Performance 

The  annual  energy  consumption  of  a  refrigeration  system  is  dependent 
on  equipment  size,  location  (for  ambient  conditions),  and  the  variable 
annual  cooling/dehumidification  load.  As  such,  energy  requirements  of  a 
system  may  vary  widely  from  one  application  to  another  and  should  therefore 
be  calculated  on  an  individual  basis. 

5.3A  Single  Coil  Systems 

Load.  To  illustrate  the  use  of  the  generalized  model  for  annual 
energy  consumption  estimates,  the  performance  of  each  of  the  systems 
included  in  the  model  was  calculated  for  the  following  two  load  conditions: 

1)  Use  of  the  system  for  pre-cooling  ambient  air  to  a  present  level. 

2)  Use  of  the  system  for  cooling/dehumidifying  recirculated  air  to  a 
preset  level. 


It  should  be  noted,  however,  that  this  includes  only  compressor  energy 
requirements.  The  annual  energy  requirements  of  the  fans  and  pumps, 
because  they  operate  generally  at  constant  flow  rate,  may  be  approximated 
by  multiplying  operational  time  by  the  design  power  consumption  of  each  of 
these  items. 
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For  both  loads  single  bin  ambient  data  for  New  Orleans,  LA  was  utilized. 

Also,   both  loads  set  55  F  as  the  temperature  to  which  the  air  was  to  be 

processed. 

System.  The  15-ton  refrigeration  systems  utilized  in  the  development 
of  the  modeling  constants  for  this  study  were  also  used  for  these 
performance  predictions.  The  equipment  modeled  is  as  follows: 

1.  Trane  CUAB  C15M  compressor  -  6  row  Aerofin  DX  coil  - 
Trane  CAUA  200  air  cooled  condenser. 

2.  Trane  CUAW  C15M  water  cooled  condensing  unit  -  6  row  Aerofin 
DX  coil  -  Baltimore  Aircoil  VXT-15  cooling  tower. 

3.  Trane  CCWB  C15M  compressor  chiller  unit  -  Trane  CAVA  200  air 
-cooled  condenser. 

4.  Trane  CGWB  C15M  chiller  unit  -  Baltimore  Aircoil  VXT-15  cooling 
tower . 

The  minimum  condensing  temperature  and  entering  condensing  water 
temperature  were  set  at  85  F  and  75  F  (which,  from  manufacturers  catalogs, 
appear  to  be  typical  values),  respectively.  Also,  the  chilled  water  was 
maintained  at  either  a  supply  temperature  of  45  F  or  a  return  temperature 
of  55  F,  (depending  on  the  system).  The  water  flow  rates  for  evaporator 
(34.5  gpm)  and  the  condenser  (44.8  gpm)  were  selected  to  provide  a  10  F 
water  temperature  swing  at  the  largest  loading  condition.  Also,  the  face 
area  for  the  DX  coil  (5  sq.  ft.)  was  chosen  to  provide  a  coil  face  velocity 
of  400  fpm  at  this  condition.   All  of  these  conditions  and  associated 


Single  bin  data  is  typically  presented  in  terms  of  dry-bulb  temperature 
(usually  grouped  into  5  degree  intervals),   the  number  of  hours  of 
occurrence  of  each  temperature  interval,  and  the  mean  coincident  wet-bulb 
temperature  of  each  interval. 
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equipment  modeling  constants  are  available  in  the  constants  subroutine  of 
the  computer  program  listed  in  Appendix  C. 

Performance.  The  performance  calculations  for  the  systems  were 
presented  in  Appendix  D.  A  summary  of  the  annual  compressor  power 
consumption  estimates  for  the  pre-cooling  load  is  presented  below. 

System 

Compressor-DX  coil-Cooling  Tower 

Compressor-DX  Coil-Air  Cooled  Condenser 

Chiller  Unit  -  Cooling  Tower 
(55  F  return) 

Compressor  Chiller  Unit-Air  Cooled 
Condenser  (55  F  return) 

Chiller  Unit  -  Cooling  Tower 
(45  F  supply) 

Compressor  Chiller  Unit-Air  Cooled 
Condenser  (45  F  supply) 

As  expected,   the  annual  compressor  power  requirements  for  DX  coil  systems 

are  lower  than  for  chilled  water  systems  (due  to  higher  evaporating 

temperatures  for  each  loading  condition).   Also,   systems  which  control 

chilled  water  return  temperature  are  more  energy  efficient  than  systems 

which  maintain  a  fixed  supply  temperature   (again  because  of  higher 

evaporator  temperatures).   Also,   systems  utilizing  cooling  towers  operate 

at  lower  annual  power  consumption  levels  than  systems  with  air  cooled 

condensers  due  to  lower  condensing  temperatures  (this,   however,   may  vary 

from  location  to  location  depending  on  ambient  dry-bulb  and  wet-bulb 

conditions).   The  second  set  of  performance  calculations  followed  similar 

trends  as  the  first  set.   It  should  be  noted,  however,  that  the  coil  load 


Annual  Power 

48300 

(kW-hrs 

0 

48690 

49480 

50090 

50390 

51790 
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conditions  for  this  calculation  run  were  arbitrarily  chosen  to  illustrate 
use  of  the  model.  As  such,  the  annual  power  consumption  values  calculated 
may  or  may  not  be  typical  of  refrigeration  systems  used  for 
cooling/dehumidifying  recirculating  air  at  this  location.  These 
performance  predictions  including  a  listing  of  temperatures,  air  and  water 
flow  rates,  and  load  and  power  values  for  each  temperature  bin  may  be  found 
in  Appendix  D. 

5.3B  Multiple  Coil  (Chilled  Water)  Systems 

As  was  discussed  in  Chapter  2,  water  chiller  systems  are  used 
primarily  for  multiple  coil  applications.  Because  the  chilled  water  system 
modeled  earlier  utilized  an  energy  balance  at  the  evaporator  to  predict 
performance,  the  generalized  model  can  easily  be  adapted  to  handle  as  many 
chilled  water  cooling  coils  as  are  desired.  For  this  situation,  the 
evaporator  load  to  be  used  in  equation  5.4  is  simply  the  sum  of  all  of  the 
individual  cooling  coil  loads  (from  the  cooling  coil  model).  From;  this, 
system  performance  can  be  calculated  in  the  same  manner  as  before. 

To  illustrate  the  use  of  the  model  for  multiple  coil  applications,  the 
generalized  computer  program  of  Appendix  C  was  modified  to  handle 
specifically  a  two  coil-chiller  unit-cooling  tower  system  (Appendix  E). 
This  program  was  then  utilized  to  make  performance  predictions  for  the 
refrigeration  system  while  providing  pre-cooling  of  ambient  make  up  air  and 
after-cooling  of  air  from. a  HONEYCOMBE  dehumidif ier;  with  the  conditions  to 
and  from  each  coil  being  provided  to  this  program  from  a  model  prediction 
for  that  specific  brand  of  dehumidif ier  [1]  (eliminating  the  need  to 
include  the  dehumidif ier  in  the  program).  The  schematic  for  the  system  is 
presented  in  Figure  5.10. 
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Load.   Performance  calculations  were  generated  for  two  different  data 

sets;  with  the  first  data  set  based  on  single  bin  ambient  data  and  the 
second  set  based  on  double  bin  ambient  data  (both  for  New  Orleans,  LA). 

System.  Based  on  the  largest  total  cooling  load,  a  20-ton  Trane 
chiller  unit  and  a  20-ton  Baltimore  AIRCOIL  cooling  tower  were  used  in  the 
computer  program.  An  evaporator  water  flow  rate  of  42.9  gpm  and  condenser 
water  flow  rate  of  55.4  gpm  were  selected.  Also,  chilled  water  supply 
temperature  was  held  at  45  F  and  the  minimum  entering  condensing  water 
temperature  at  75  F. 

Performance.  The  performance  calculations  (along  with  coil-air 
conditions)  are  presented  in  Appendix  F.  A  summary  of  system  annual 
performance  values  for  each  data  sets  is  presented  below. 

Annual 


Single  Bin 

Double  Bin 

123.63 

133.66 

838.88 

858.72 

962.51 

992.33 

76801 

79689 

8417 

8537 

Pre-cool  coil  load  (MMBTU) 
After-cool  coil  load  (MMBTU) 
Total  coil  load  (MMBTU) 
Compressor  power  (KW-HR) 
Hours  of  operation 

Examination  of  this  table  shows  that  all  of  the  annual  performance  values 

for  the  single  bin  are  consistently  lower  than  the  double  bin  performance 

values.  The  single  bin  annual  power  consumption,  however,  was  within  4%  of 

the  double  bin  value.   Thus,   the  use  of  single  bin  ambient  data  (which 

requires  less  processing  time  due  to  fewer  data  points  than  with  double  bin 

data)   may   provide  adequate  power  consumption  estimates   for   most 

refrigeration  system  applications. 


Double  bin  data  is  typically  presented  in  terms  of  dry-bulb  temperature 
(usually  grouped  into  5  degree  intervals)  and  the  hours  of  occurrence  of 
each  wet-bulb  temperature  (usually  grouped  into  2  degree  intervals)  for 
each  dry-bulb  temperature  interval. 
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CHAPTER  VI 
CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  Summary  and  Conclusions 

The  main  objective  of  this  study  was  to  develop  a  generalized  model 
for  predicting  the  annual  energy  requirements  of  reciprocating  vapor- 
compression  refrigeration/dehumidif ication  systems.  The  annual  performance 
of  a  system  is  directly  dependent  on  location  and  the  variable  load  to 
which  the  system  is  to  be  subjected.  The  model  developed  is  capable  of 
handling  the  following  systems: 

1.  Compressor  -  DX-coil  -  air  cooled  condenser  systems. 

2.  Water  cooled  condensing  unit  -  DX-coil  -  cooling  tower  systems. 

3.  Compressor  chilled  unit  -  chilled  water  coil(s)  --air  cooled 
condenser  systems. 

4.  Chiller  unit  -  chilled  water  coils(s)  -  cooling  tower  systems. 

For  each  of  these  systems,   cylinder  unloading  was  used  for  partial  load 
system  capacity  control.   The  model  agrees  with  catalog  data  (+3-1/2%)  for 
a  wide  range  of  equipment  variations  and  sizes.   Also,   the  use  of 
subsystem  units  in  the  model  greatly  reduces  the  amount  of  data  required 
for  individual  system  performance  modeling. 

6.2  Recommendations 

It  is  recommended  that: 

1.  A  procedure  be  developed  for  including  system  annual  energy 
requirements  and  costs  into  the  equipment  selection  procedure. 
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2.  The  generalized  model  be  extended  to  other  refrigeration  systems 
utilizing  screw  and  centrifugal  compressors. 

3.  The  model  be  verified  experimentally  to  determine,  if  any,  the 
shortcomings  of  the  model  in  making  annual  energy  consumption 
estimates  for  actual  operation  systems. 
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Equation  Summary  -  Component  Model  -  Full  Load  Operation 

Compressor 

%   "  al  +  a2Te  +  a3Te2  +  a4Tc  +  a5TC2  +  a6TeTc  +  a7Te\ 

7  2   2 

+  aaT  T  *   +  aQT  ZT 
sec     9  e  c 


P  "  bl  +  b2Te  +  b3Te2  +  b4Tc  +  b5Tc2  +  Wc 
+  b7Te2Tc  +  b8TeTc2  +  b9Te2Tc2 


Condenser 

Qc  =  Qe  ♦  P 

%  •   Fl  +  ¥TC  "  V 
subcooling  =  F3  +  F4(Tc  -  7 J   +  F5(T.  -  Tm)2 

air-cooled  T  =  ambient  dry  bulb  temperature 

m 

water-cooled  T  =  inlet  water  temperature 

n 

Evaporator 

DX-Coil 

Qg  =  (CF)(AREA)[C1  +  CjT,  ♦  C3Te2  +  C4Twb  +  C5Twb2  ♦  C^T^  +  C^h^ 

CDT  T  .  2  +  CQT  2T  2] 
8  e  wb     9  e  wb 

CF  =  C10  +  C.j  fpm  +  C12  fpm2 
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Chilled-water 


Q  =  d1  +  d2(GPM)  +  d3(GPM)2  +  d4(dT)  +  d5(dT)2  +  d&(GPM)(dT) 

+  d?(GPM)2(dT)  +  dg(GPM)(dT)2  +  dg(GPM)2(dT)2 
dT  =  (inlet  water  temperature)-(saturated  evaporating  temperature.) 


where  a,.~ag  • 

bl-b9 

FrF5 
crci2 
dr  d9 


=  constants  to  be  found  for  each  system 


(generally  by  a  least-squares  fit  of  each 

equation  to  manufacturers  catalog  data) 


• 
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Equation  Summary  -  Leverenz-Bergan  Chiller  Model  -  Full  Load  Operation 


T   =  (T   ,  -  T   .  -)/T  _.  -  (T   -  T   „) 
eq     cond    cond.D'   ratio     cw    cw,D' 

ADCR  =  a,  +  a,(T  )  +  a,(T  )2 
1    2V  eq'    3V  eq' 

CAP  =  (ADCR)(CAP,D) 

COPC  =  (T   ,  +  273.2)/(T   .  -T  ) 
cond  cond   cw' 

COPCD  =  (T   ,  „)  +  273.2)/T   .  -  -  T   n) 
cond.D'  cond.D    cw,D 

COPCN  =  (COPC) /(COPCD) 

n 
COPAN  -  b,  +  b2(C0PCN)  +  b3(C0PCN) 

COPAD  =  (CAP,D)/(FLP,D) 

COPA  =  (C0PAN)(C0PAD) 

FLP  =  (CAP)/(C0PA) 

where 

all  input  and  output  in  Kilowatts  and  degrees  Celsius 

input  - 

T   =  leaving  chilled  water  temperature 

T   ,  =  entering  condensing  water  temperature 

T  ,  .  ,  T   _,  T   j  n>  a^aii  b-.'^b^   =  constants  for  all  chiller 
ratio   cw,D'   cond,D   1   3'   1   3 

systems 

CAP,D  =  System  capacity  at  T   „  and  T    ,  _  temperatures 

FLP,n  =  System  power  at  T   _  and  T    ,  „  temperatures 
D    '     r        cw,D      cond,D    r 

output  - 

CAP  =  predicted  system  capacity 

FLP  =  predicted  system  power 

(for  definitions  of  T  ,  ADCR,  See  subsystem  model  -  equation  summarv) 

eq  '  n  '' 
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Equation  Summary  -  subsystem  model  -  Full  Load 


eq 


ADCR 


(T 


el 


Tcl,D^/Tratio  "  'Tei 


el  ,D; 


^  *  A2(Teq)  *  A3(T 


eq' 


CAP  =  (ADCR)(CAPD) 
COPC  -  (T ,  +  273.2)/(Tcl  -  Tel) 
COPCD=(TcljD+273.2)/(Tcl)D-Tel>D) 
COPCN  =  COPC/COPCD 
COPAN  =  B1  +  B2(C0PCN)  +  B3(C0PCN)2 
COPAD  =  (CAP,D)/(FLP,D) 


COPA  =  (COPAN) (COPAD) 
FLP  =  CAP/COPA 
where 

compressor 


el 


el,D 


cl 


cl,D 


where 


evap.  temp 

base 
evap.  temp 

cond.  temp 


base 
cond.  temp 


water-cooled 
condensing  unit 


evap .  temp 

base 
evap.  temp 


leaving  cond. 
water  temp. 

base 


leaving  cond.   cond.  temp 
water  temp. 


compressor- 
chiller  unit 

leaving  evap . 
water  temp. 

base 
leaving  evap. 
water  temp 

cond.  temp. 


base 


chiller  unit 


leaving  evap. 
water  temp . 


base 
leaving  evap. 
water  temp 

leaving  cond. 
water  temp. 

base 
leaving  cond. 
water  temp 


all  input  and  output  in  Kilowatts  and  degrees  Celsius 
input  -  T  ,  and  T  ,  temperatures 

Tratio'  Tel  D'  Tcl  D*  ^1^3'  B1~B3  =  constants  for  appropriate 

system  (Table  3. A) 
CAP,_  =  System  capacity  at  T  .  .  and  T  .  _  temperatures 
FLP,D  =  System  power  at  T^p  and  T^  D  temperatures 
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output  -  CAP  =  predicted  system  capacity 

FLP  =  predicted  system  power 

where 

T   =  Equivalent  operational  temperature 
eq 

T      =  ratio  of  change  in  condenser  temperature  (refrig.  or  water) 
ratio 

required  for  a  given  change  in  evaporator  temperature  (refrig. 
or  water)  which  maintains  a  given  capacity. 
ADCR  =  Actual  versus  Design  (base)  Capacity  Ratio 
COPC  =  Coefficient  of  Performance,  Carnot 
COPCD  =  Coefficient  of  Performance,  Carnot  at  Design  (base)  temperatures 
COPCN  =  Normalized  Carnot  Coefficient  of  Performance 

COPA  =  Coefficient  of  Performance,  Actual 
COPAD  =  Coefficient  of  Performance,  Actual  at  Design  (base)  temperatures 
COPAN  =  Normalized  Actual  Coefficient  of  Performance 
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APPENDIX  B 


COMPARISON  OF  SUSBSYSTEM 
MODELING  ACCURACY 


CONPRESSOR  CHARACTERISTICS  77 


BASE 

EVAPORATOR  TEMP. 

(deg.F) 

40.0 

BASE 

CONDENSER 

TEHP. 

(deg.F) 

105.0 

BASE 

CONPRESSOR 

LOAO 

(tons) 

17.3 

BASE 

CONPRESSOR  POWER 

■■•»s! 

15.9 

Tel 

Id 

Qejred 

Qe  cat 

dif. 

Pored 

P  cat 

dif. 

. .  I'deg.F) .. 

....(tans).... 

1 

..(k«s>... 

I 

30.0 

85.0 

16.2 

16.0 

1.1 

13.2 

13.1 

.8 

35.0 

85.0 

17.9 

17.8 

.8 

13.4 

13.4 

.3 

40. 0 

85.0 

19.8 

19.7 

.4 

13.6 

13.6 

.1 

45.0 

85.0 

21. a 

21.7 

.5 

13.8 

13.8 

-.2 

50.0 

85.0 

23.9 

23.9 

.1 

14.0 

13.9 

.4 

30.0 

105.0 

13.9 

14.0 

-.5 

15.2 

14.9 

1.3 

35.0 

105.0 

15.5 

15.4 

-.4 

15.6 

15.4 

1.6 

40.0 

105.0 

17.3 

17.3 

-.3 

16.0' 

15.9 

.9 

45.0 

105.0 

19.1 

19.2 

-.6 

16.4 

16.4 

-,2 

50.0 

105.0 

21.0 

21.1 

"■%) 

16.6 

16.9 

-1.7 

30.0 

125.0 

11.9 

11.9 

.1 

16.7 

16. a 

-.6 

35.0 

125.0 

13.4 

13.3 

.4 

17.4 

17.6 

-1.0 

40.0 

125.0 

14.9 

14.9 

.1 

18.1 

18.3 

-1.1 

45.0 

125.0 

16.4 

14.6 

-.1 

18.7 

19.0 

-1.5 

50.0 

125.0 

18.4 

18.3 

.4 

19.2 

19.7 

-2.4 

nhere 


Tel  -  saturated  eyanerating  teioerature 

Tel  ■  saturated  condensing    teeoerature 

Qejred  =  capacity  predicted  by  lodel 

3e_:at  "  lanutacturers  catalog  data 

PJired  ■  poner  requiretents  predicted  by  lodel 

P.cat  •  poiier  requireients  froi  catalog 


constants 


Tratio  ■  2.9 
Al  =  .9973 
A2  =-.03694 
A3  =  .0004335 
Bl  =  -.4141 
B2  =  1.7717 
B3  ■  -.3694 


COMPRESSOR  CHARACTERISTICS 


BASE 

EVAPORATOR  TEMP. 

(deq.F) 

40.0 

BASE 

CONDENSER 

TEHP. 

(deg.F) 

105.0 

BASE 

COMPRESSOR  LOAD 

(tons) 

34.1 

BASE 

COMPRESSOR  PONER 

!k«s) 

29.3 

Tel 

Tel 

Qeored 

Be  cat 

dif. 

Pjrsd 

P  cat 

dif. 

..(de 

g.F).. 
35.0 

....(ton 
31.? 

It.... 

I 
-1.0 

,.(tM 

>)7.. 

24.0 

I 

30.0 

32.2 

24.3 

3.1 

35.0 

85.0 

35.4 

35.6 

-.7 

25.2 

24.6 

2.4 

40.0 

85.0 

39.1 

39.2 

-.4 

25.5 

25.2 

1.3 

45.0 

85.0 

43.0 

43.0 

-.0 

25.8 

25.6 

.9 

50. 0 

85.0 

47.1 

47.0 

.3 

26.2 

26.1 

.3 

30.0 

105.0 

27.5 

27.9 

-1.6 

28.4 

28.1 

1.2 

35.0 

105.0 

30.6 

31.0 

-1.2 

29.3 

29.0 

1.1 

40.0 

105.0 

34.0 

34.1 

-.3 

30.1 

29.3 

.9 

45.0 

105.0 

37.4 

37.5 

.3 

30.7 

30.5 

.6  - 

50.0 

105.0 

41.5 

41.0 

1.1 

31.1 

31.1 

.1 

30.0 

125.0 

23.5 

23.7 

-.9 

31.3 

31.6 

-.9 

35.0 

125.0 

24.3 

26.3 

.1 

32.7 

32.9 

-.7 

40.0 

125.0 

29.4 

29.2 

.7 

33.9 

34.1 

-.5 

45.0 

125.0 

32.7 

32.1 

1.9 

35.1 

35.2 

-.4 

50.0 

125.0 

36.2 

35.3 

2.6 

36.0 

36.3 

-.7 
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BASE 

COMPRESSOR    CHARACTERISTICS 

40.0 

EVAPORATOR  TEMP. 

(deg.F) 

BASE 

CONDENSER 

TEMP. 

(deg.F) 

105.0 

BASE 

COMPRESSOR  LOAD 

(tons) 

91.1 

BASE 

COMPRESSOR  POKER 

(k«s) 

a4.7 

Tel 

Tel 

Sepred 

Qe  cat 

dif. 

Pared 

P  cat 

dif. 

..  (deq.F).. 

....(tons)7... 

I 

..(k«s>:.. 

Z 

30.0 

85.0 

75.3 

74.9 

;.: 

53.7 

52.5 

2.4 

35.0 

85.0 

84.1 

83.3 

.9 

54.7 

53.5 

2.2 

40.0 

85.0 

92.9 

92.1 

.3 

55.4 

54.4 

i.a 

45. 0 

85.0 

102.2 

101.5 

.7 

54.0 

55.0 

!.? 

50.0 

85.0 

112.1 

111.3 

.7 

54.9 

55.4 

2.5 

30.0 

105.0 

65.3 

65.3 

.0 

41.9 

41.2 

.9 

35.0 

105.0 

72.3 

73.0 

-.2 

43.7 

63.1 

.7 

40.0 

105.0 

90.9 

81.1 

-.3 

45.3 

64.7 

.? 

45.0 

105.0 

89.5 

89.4 

-.1 

44.6 

66.0 

.9 

50.0 

105.0 

98.4 

93.7 

-.1 

67.6 

67.1 

.7 

30.0 

125.0 

55.9 

55.0 

1.6 

68.0 

68.5 

-.8 

35.0 

125.0 

62.6 

42.3 

.5 

70.9 

71.6 

-.9 

40.0 

125.0 

69.9 

49.4 

.5 

73.7 

74.5 

-1.1 

45.0 

125.0 

77.8 

77.4 

.5 

76.1 

77.3 

-1.5 

50.0 

125.0 

B6.2 

95.4 

.0 

78.3 

79.9 

-2.1 

COMPRESSOR-CONDENSER  CHARACTERISTICS 
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BASE  EVAPORATOR  TEHP.  (deg.F)  =  40.0 

BASE  CONDENSER. HATER  TENP.  fdeg.F)  =  95.0 

BASE  CONPRESSOR  LOAD  (tons)  =  23.7 

BASE  COMPRESSOR  POKER  (k»s)  =  17.3 


Tel 

Tel 

Qe  ored 

Qe  cat 

Hi. 

Pored 

P  cat 

dif . 

..(deg 

.F).. 

....(tors).... 

1 

..(knsl... 

1 

30.0 

85.0 

20.4 

20.5 

-.5 

15.2 

14.9 

1.9 

35.0 

85.0 

22.9 

23.0 

-.6 

15.6 

15.4 

1.5 

40.0 

85.0 

25.5 

25.6 

-.4 

16.0 

15.8 

1.1 

45.0 

85.0 

28.3 

23.4 

-.3 

16.2 

16.2 

.2 

50. 0 

85.0 

31.3 

31.4 

-.3 

16.4 

16.4 

.2 

30.0 

95.0 

18.8 

18.9 

-.5 

16.1 

16.0 

.3 

35.0 

95.0 

21.1 

21.2 

-.3 

16.8 

16.7 

.4 

40.0 

95.0 

23.7 

23.7 

-.2 

17.3 

17.3 

.1 

45.0 

95.0 

26.3 

26.3 

1 

17.7 

17.3 

-.3 

50.0 

95.0 

29.2 

29.2 

A 

18.1 

18.2 

-.8 

30.0 

105.0 

17.3 

17.2 

.7 

16.9 

16.9 

-.1 

35.0 

105.0 

19.5 

19.4 

.6 

17.7 

17.3 

-.5 

40.0 

105.0 

21.9 

21.3 

.5 

18.5 

18.6 

-.3 

45.0 

105.0 

24.5 

24.3 

.7 

19.1 

19.3 

-1.1 

50.0 

105.0 

27.2 

27.0 

.8 

19.6 

19.9 

-1.5 

•here 


Tel  ■  saturated  evaperating  teiperature 

Tel  ■  leaving  condensing  »ater  teiperature 

Qejred  =' capacity  predicted  by  lodel 

Qe  cat  ■  lanufacturers  catalog  data 

Pjred  =  poner  requiretents  predicted  by  lodel 

P_cat  ■  poner  requirements  fro«  catalog 


constants 


Tratio  =  2.9 
Al  =  .9932 
A2  =-.03952 
A3  =  .004931 
Bl  =  -.2244 
B2  =  1.5438 
B3  •    -.3219 


COMPRESSOR-CONDENSER  CHARACTERISTICS 


BASE 

EVAPORATOR 

TEMP. 

(deg.F) 

=  40.0 

BASE 

CONDENSER  WATER  TEMP 

.   (deg.F) 

=  95.0 

BASE 

COMPRESSOR  LOAD 

(tons) 

=  38.1 

BASE 

COMPRESSOR 

POSER 

<kns) 

=  28.4 

Tel 

Tci 

Oejred 

Qe  cat 

dif. 

Pjred 

P  cat 

at. 

..  Ideg.F).. 

....(tons).... 

I 

..(km 

i)7.. 

i 

30.0 

85.0 

32.8 

33.2 

-1.2 

24.9 

24.5 

1.7 

35.0 

S5.0 

34.7 

36.9 

-.4 

25.7 

25.2 

1.8 

40.0 

85.0 

41.0 

40.9 

.2 

26.2 

25.3 

1.7 

45.0 

85.0 

45.5 

45.2 

.7 

26.6 

26.3 

1.3 

50.0 

85.0 

50.3 

49.8 

1.1 

27.0 

26.6 

1.4 

30.0 

95.0 

30.2 

30.8 

-1.8 

26.5 

26.6 

-.5 

35.0 

95.0 

34.0 

34.3 

-.9 

27.5 

27.5 

.1 

40.0 

95.0 

38.0 

38.1 

-.2 

2B.4 

28.4 

.1 

45.0 

95.0 

42.4 

42.2 

.4 

29.1 

29.1 

.1 

50.0 

95.0 

47.0 

46.5 

1.0 

29.6 

29.8 

-.4 

30.0 

105.0 

27.  a 

28.5 

-2.3 

27.7 

28.5 

-2.8 

35.0 

105.0 

31.4 

31.9 

-1.6 

29.1 

29.7 

-2.1 

40.0 

105.0 

35.2 

35.4 

-.5 

30.3 

30.3 

-1.7 

45.0 

105.0 

39.3 

39.3 

.1 

31.3 

31.9 

-1.8 

50.0 

105.0 

43.8 

43.3 

1.1 

32.2 

32.8 

-1.9 

COMPRESSOR-CONDENSER  CHARACTERISTICS 


BASE  EVAPORATOR  TEMP.  (deg.F)  =  40.0 

BASE  CONDENSER  WTER  TEMP.  (deg.F)  =  95.0 

BASE  COMPRESSOR  LOAD  (tons)  *  88.4 

BASE  COMPRESSOR  POKER  (kits)  *  66.6 


Tel 

Tci 

Qe  jred 

3e  cat 

dif. 

Pored 

P  cat 

dif. 

. . :Jeg 

.Fl.. 

....(tons).... 

I 

..(k»s)r.. 

I 

30.0 

35.0 

76.1 

77.7 

-2.1 

58.4 

57.9 

.9 

35.0 

35.0 

85.3 

36.3 

-1.2 

60.2 

59.4 

1.3 

40.0 

B5.0 

95.1 

95.4 

-.3 

61.5 

60.3 

1.2 

45.0 

35.0 

105.6 

105.1 

.5 

42.5 

41.3 

1.1 

50.0 

85.0 

116.8 

115.4 

1.2 

43.2 

62.4 

1.0 

30.0 

95.0 

70.2 

71.9 

-2.4 

42.1 

42.5 

-.4 

35.0 

95.0 

78.9 

79.9 

-1.3 

44.6 

44.7 

-.2 

40.0 

95.0 

38.2 

83.4 

-  j 

44.6 

44.6 

.1 

45.0 

95.0 

98.3 

97.4 

.7 

68.3 

68.2 

.1 

50.0 

95.0 

109.0 

107.3 

1.6 

69.5 

69.5 

-.0 

30.0 

105.0 

64.6 

66.2 

-2.5 

45.0 

44.7 

-2.4 

35.0 

105.0 

72.8 

73.7 

•1.2 

48.2 

49.5 

-1.9 

40.0 

105.0 

81.7 

81.3 

-.1 

71.0 

72.0 

-1.3 

45.0 

105.0 

91.3 

90.3 

1.1 

73.5 

74.2 

-1.0 

50.0 

105.0 

101.5 

99.5 

2.0 

75.5 

74.1 

-.8 

COMPRESSOR  CHILLER  UNIT  CHARACTERISTICS 
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BASE 

CHILLED  MATER  TOP. 

I4m.F1 

=  46.0 

BASE 

CONDENSER 

TEMP. 

(deg.F) 

=105.0 

BASE 

COMPRESSOR  LOAD 

(tons) 

«  16.1 

BASE 

COMPRESSOR  POKER 

IkHS) 

=  15.6 

Tel 

Td 

Sepred 

Be  cat 

dif. 

Pjred 

P  cat 

dif. 

..(deg.F).. 

....(tons).... 

I 

..(k»s)... 

I 

40.0 

35.0 

16.3 

16.2 

.7 

13.2 

13.2 

.2 

44.0 

85.0 

17.4 

17.3 

.6 

13.3 

13.3 

.2 

46.0 

S5.0 

18.0 

17.9 

.3 

13.4 

13.4 

-.1 

50.0 

85.0 

19.1 

19.1 

-.1 

13.5 

13.5 

.3 

40.0 

95.0 

15.4 

15.3 

.6 

14.3 

14.1 

1.2 

44.0 

95.0 

16.5 

16.5 

-.2 

14.4 

14.3 

.7 

46.0 

95.0 

17.0 

17.0 

.0 

14.4 

14.4 

.3 

50.0 

95.0 

18.1 

18.2 

-.5 

14.5 

14.7 

-1.2 

40.0 

105.0 

14.5 

14.5 

-.0 

15.3 

15.0 

l.S 

44.0 

105.0 

15.5 

15.6 

-.4 

15.5 

15.4 

.5 

46.0 

105.0 

16.1 

16.1 

-.2 

15.6 

15.6 

-.3 

50.0 

105.0 

17.2 

17.2 

"■%) 

15.7 

15.9 

-1.4 

40.0 

115.0 

13.6 

13.6 

.1 

16.1 

16.1 

.3 

44.0 

115.0 

14.6 

14.7 

-.4 

16.5 

16.5 

-.3 

46.0 

115.0 

15.2 

15.2 

-.3 

16.6 

16.7 

-.7 

50.0 

115.0 

16.2 

16.3 

-.5 

16.8 

17.1 

-1.7 

Tel  =  leaving  chilled  »ater  tesperature 

Td  ■  saturated  condensing    teiperature 

Qejred  =  capacity  predicted  oy  todel 

3e_cat  ■  •anufacturers  catalog  data 

Pjred  •  po»er  reguire«ents  predicted  by  »odel 

P  cat  =  po«er  reguireients  froi  catalog 


constants 


Tratio  =  2.9 
Al  =  .9982 
A2  =-.02990 
A3  =  .0001711 
Bl  =  -.2484 
B2  =  1.6600 
B3  =  -.4105 


COMPRESSOR  CHILLER  UNIT  CHARACTERISTICS 
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BASE 

CHILLED  KATER  TEMP. 

(deg.F) 

=  44.0 

BASE 

CONDENSER 

TEMP. 

(deg.F) 

=105.0 

BASE 

COMPRESSOR  LOAD 

(tons! 

=  32.0 

BASE 

COMPRESSOR  POKER 

(kws) 

=  29.3 

Tel 

Tel 

Qepred 

Qe  cat 

dif. 

Pjred 

P  cat 

dif. 

.. (dsg.Fl.. 

....(tans):... 

I 

..(kNS)... 

I 

40.0 

85.0 

32.4 

32.5 

-.2 

24.9 

24.1 

3.1 

44.0 

85. 0 

34.4 

34.7 

-.3 

25.0 

24.5 

2.2 

46.0 

85.0 

35.7 

35.8 

-.3 

25.1 

24.4 

2.2 

50.0 

85.0 

37.9 

33.1 

-.4 

25.4 

25.0 

1.8 

40.0 

95.0 

30.4 

30.7 

-.3 

24.8 

24.3 

1.9 

44.0 

95.0 

32.7 

32.8 

-.3 

27.0 

24.8 

.9 

14.0 

95.0 

33.3 

'33.9 

-.3 

27.1 

27.0 

.5 

50.0 

95.0 

34.0 

36.1 

-.3 

27.3 

27.4 

-.5 

40.0 

105.0 

28.8 

28.9 

-.3 

28.7 

28.4 

1.0 

44.0 

105.0 

30.9 

30.9 

-.0 

29.1 

29.0 

.2 

44.0 

105.0 

31.9 

32.0 

-  2 

29.2 

29.3 

-.3 

50.0 

105.0 

34.1 

34.1 

-,0 

29.5 

29.8 

-1.2 

40.0 

115.0 

27.1 

27.1 

-.1 

30.3 

30.5 

-.4 

44.0 

115.0 

29.1 

29.0 

.3 

30.9 

31.2 

-.9 

44.0 

115.0 

30.1 

30.0 

.4 

31.2 

31.5 

-1.1 

50.0 

115.0 

32.2 

32.0 

.7 

31.4 

32.1 

-1.4 

COMPRESSOR  CHILLER  UNIT  CHARACTERISTICS 


BASE 

CHILLED  HATER  TEMP. 

(deg.F) 

=  44.0 

BASE 

CONDENSER  TEMP. 

(deg.F) 

=105.0 

BASE 

COMPRESSOR  LOAD 

(tons) 

=  78.0 

BASE 

COMPRESSOR  POKER 

(kNS) 

=  44.1 

Tel 

Tc: 

Qeored 

fle  cat 

dif. 

Pored 

P  cat 

dif. 

..(deg.F).. 

....(tons).... 

I 

..(kxsi;.. 

Z 

40.0 

35.0 

79.0 

79.4 

.a 

24. 4 

53.0 

2.6 

44.0 

85.0 

34.3 

84.0 

.4 

54.8 

53.4 

2.2 

44.0 

35.0 

87.0 

84.9 

.1 

55.0 

53.9 

2.0 

50.O 

35.0 

92.5 

92.7 

-.3 

55.7 

54.4 

2,3 

40.0 

95.0 

74.4 

74.3 

.4 

58.7 

57.8 

1.3 

44.0 

95.0 

79.7 

79.7 

.1 

59.2 

58.7 

.3 

44.0 

95.0 

32.4 

82.5 

-.1 

59.3 

59.0 

.4 

50.0 

95.0 

37.7 

88.2 

-.5 

59.7 

59.7 

-.1 

40.0 

105.0 

70.2 

70.1 

.2 

42.7 

42.4 

.5 

44.0 

105.0 

75.3 

75.3 

-.0 

63.4 

63.6 

-.0 

44.0 

105.0 

77.9 

78.0 

-.2 

43.9 

64.1 

-.3 

50.0 

105.0 

83.1 

83.5 

-.5 

64.4 

65.1 

-1.0 

40.0 

115.0 

66.0 

45.7 

.4 

44.3 

66.3 

-.7 

44.0 

115.0 

70.9 

70.3 

5 

47.4 

68.4 

-1.1 

44.0 

115.0 

73.4 

73.4 

!i 

48. 2 

69.1 

-1.3 

50.0 

115.0 

78.6 

73.7 

-.2 

W.l 

70.5 

-2.0 
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BASE 

CHILLER 

CHARACTERISTICS 

Fl 

=  46.0 

CHILLED  HATER  TEPtP. 

<deg 

BASE 

CONDENSER  HATER  TEHP 

idea 

F) 

■  95.0 

BASE 

CONPRE5S0R  LDAD 

(tans) 

=  15.8 

BASE 

COMPRESSOR  POWER 

(Ins) 

=  16.0 

Tel 

Tel 

Dejred 

Be  cat 

dif. 

Pjred 

P  cat 

tit. 

..MH.FI.. 

....(tons) 7... 

I 

..(km 

)7.. 

1 

40.0 

85.0 

15.1 

15.1 

14.7 

14.4 

2.0 

44.0 

85.0 

16.1 

16.1 

14.9 

14.3 

.7 

16.0 

85.0 

16.7 

16.6 

15.0 

15.0 

.1 

50.0 

85.0 

17.7 

17.7 

15.3 

15.3 

-.0 

40.0 

95.0 

14.2 

14.2 

15.6 

15.3 

2.0 

44.0 

95.0 

15.3 

15.3 

- 

15.9 

15.8 

.4 

46.0 

95.0 

15.8 

15.8 

- 

16.0 

16.0 

-  <J 

50.0 

95.0 

16.8 

16.8 

16.2 

16.4 

-i!3 

40.0 

105.0 

13.4 

13.4 

- 

16.4 

16.4 

T 

44.0 

105.0 

14.4 

14.4 

- 

16.8 

16.9 

-!6 

46.0 

105.0 

14.9 

14.9 

- 

17.0 

17.1 

-.9 

50.0 

105.0 

15.9 

15.9 

17.2 

17.6 

-2.2 

Tel  =  leaving  evaporator  «ater  teiperature 

Tel  =  leaving  condenser  water  temperature 

Qejred  ■  capacity  predicted  by  iodel 

8e_cat  ■  lanufaclurers  catalog  data 

P_pred  =  power  requireients  predicted  by  iodel 

P_cat  =  power  requireients  froi  catalog 


constants 


Tratio  =  2.9 
Al  =  .9978 
A2  =-.02943 
A3  =  .0001472 
Bl  =  -.1119 
B2  =  1.5815 
B3  ■  -.4702 


CHILLER 

CHARACTERISTICS 

BASE 

CHILLED  HATER  TEMP. 

(d«9, 

,F)  =  44, 

,0 

BASE 

CONDENSER 

HATER  TEMP 

Idea, 

,F)  =  95.0 

BASE 

COMPRESSOR  LOAD 

(tone)  =  31, 

,2 

BASE 

COMPRESSOR  PGMER 

(k»s)  =  30 

,1 

Tel 

Tel 

Qe  pred 

Qe  cat 

di-f. 

Pored 

P  cat 

dif. 

..(deg.F).. 

....(tons):... 

I 

.. (k*s)7.. 

I 

40.0 

95.0 

29.8 

30.1 

-.9 

27.4 

27.1 

2.0 

44.0 

85.0 

31.9 

32.1 

-.7 

28.0 

27.9 

.9 

46.0 

85.0 

32.9 

33.1 

-.4 

28.3 

29.1 

.4 

50.0 

85.0 

35.0 

35.1 

~.3 

28.9 

28.7 

.3 

40.0 

95.0 

29. 1 

2B.4 

-1.0 

29.4 

29.0 

1.3 

44.0 

95.0 

30.1 

30.3 

-.4 

29.9 

29,8 

.2 

46.0 

95.0 

31.1 

31.2 

-.2 

30.1 

30.1 

_  T 

50.0 

95.0 

33.2 

33.2 

-.0 

30.4 

30.9 

-U2 

40.0 

105.0 

24.4 

24.6 

-.6 

30.9 

30.9 

.1 

44.0 

105.0 

28.4 

28.4 

-.0 

31.6 

31.7 

-.3 

44.0 

105.0 

29.4 

29.4 

-.0 

31.9 

32.1 

-.7 

50.0 

105.0 

31.4 

31.2 

.7 

32.4 

32.9 

-1.6 

BASE 

CHILLER 

CHARACTERISTICS 

=  46.0 

CHILLED  HATER  TEMP. 

(deq.F! 

BASE 

CONDENSER  HATER  TEMP 

,   (deq.F) 

=  95.0 

BASE 

COMPRESSOR  LOAD 

(tons) 

=  75.9 

BASE 

COMPRESSOR  POHER 

(tM) 

=  65.0 

Tel 

Tel 

fleored 

Se  cat 

dif.    P  ored 

P  cat 

itf, 

..(deg.F).. 

....(tons)7... 

I 

..(k»s)... 

z 

40.0 

95.0 

72.4 

72.4 

.3 

59.7 

58.6 

1.1 

44.0 

95.0 

77.5 

77.5 

.0 

60.6 

59.9 

1.3 

44.0 

95.0 

80,1 

80.2 

-.2 

61.0 

60.4 

1.0 

50. 0 

95.0 

85.2 

85.6 

-.5 

62.1 

61.4 

1.2 

40.0 

95.0 

68.4 

68.4 

.0 

63.4 

62.8 

1.0 

44.0 

95.0 

73.3 

73.4 

-.2 

64.4 

64.3 

,2 

44.0 

95  J) 

75.7 

75.9 

-.2 

64.9 

65.0 

_  2 

50.0 

95.0 

90.8 

81.1 

-.4 

65.7 

66.4 

-ill 

40.0 

105.0 

64.3 

64.3 

.0 

66.8 

66.9 

-.2 

44.0 

105.0 

69.1 

69.0 

.1 

68.2 

69.7 

-.7 

46.0 

105.0 

71.5 

71.5 

-.0 

68.9 

69.6 

-1.1 

50.0 

105.0 

76.4 

76.5 

-.1 

69.9 

71.4 

-2.0 
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APPENDIX  C 


GENERALIZED 
RECIPROCATING  REFRIGERATION/DEHUMIDIFICATION 
SYSTEM  COMPUTER  PROGRAM 


10 
20 
SO 

JO 
:0 
60 
70 

ao 
M 
100 
110 

ne 

120 

no 

140 

ISO 

160 
!70 
1B0 
190 
200 
210 
220 

no 

240 
250 
260 
270 
2S0 
290 
300 
310 
320 
330 
340 
350 
340 
370 
380 
390 
400 
410 
420 
430 
440 
450 
440 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
590 
590 
600 
610 
620 
630 
640 
650 
660 
470 
490 
700 
710 
720 

no 

740 


rem 

REM 

REN 

REN 

REN 

REN 

REN 

REN 

REN 
REN 
REM 

REN 


GENERALIZED  PROGRAM  FOR  REFRIGERATION  GYSTEN  PERFORMANCE  PREDICTIONS  86 

FOR  THE  FOLLOHINS  SYSTENS: 

C0MPRE5SOR--0I  COIL-RENOTE  AIR  COOLED  CONDENSER 
HATES  COOLED  CONDENSING  UNIT—OX  COIL-COOLING  TOMER 
CONPRESSOR  CHILLER  UNIT-CHILLED  BATES  COIL-RENOTE  SIR  COOLED  CONDENSES 
CHILLER  UNIT-CHILLED  WATER  COIL-COOLING  TOMES 
KITH  CYLINDER  UNLOADING  MATCHED  TO  HOLD  DESIRED  DX-COIL  LEAVING  AIR  TEMP         -OR- 
HITH  CYLINDER  UNL0ADIN6  MATCHED  TO  HOLD  DESISED  EVAPORATOR  LEAVING  OR  ENTERIN6  SATES  TEMP 
CONSTANT  KATES  FLOW  RATE  AT  EVAPORATORS, CONDENSERS  (HATER  BYPASS — IF  NECE3SARY--AT  COIL! 

THE  CONSTANTS  DESCRIBING  EQUIPMENT  SIZE  ARE  ALREADY  INCLUDED  IN  THE  CONSTANT  SUBROUTI 

THE  COIL  CONDITIONS  ARE  PROVIDED  TO  THE  PROGRAM  IN  THE  COIL  CONDITIONS  SUBROUTINE 


DIM  Osens(20O), Slat (200), Bb(ZOO) ,P(200) ,brs(200) ,TOB(2001 .THBC200) 

DIM  Tdb(200) ,«b.»i (200) ,CFNe(200l ,Ts(200l ,UNITS*(200) 


INDICATE  THE  TYPE  OF  SYSTEM  TO  BE  UTILIZED  BY  SETTING  SYSTEM  E3UAL  TO  THE  APPROPRIATE  VALUE 

1  =       COMPRESSOS  UNIT DX  COIL BEMOTE  AIR-COOLED  CONDENSES 

2  =  NATES  COOLED  CONDENSING  UNIT— OX  COIL COOLING  TOMER 

3  =  COMPRESSOS  CHILLES  UNIT WATER  COIL REMOTE  AIR-COOLED  CONDENSER 

4  =         CHILLES  UNIT HATER  COIL COOLING  TOMER 


SYSTEMS 


SELECTED  FROM  THE  ABOVE  VALUES 


SOSUB  CONSTANTS  ! 
SOSUB  COIL  CONDITIONS  ! 
SOSUB  HEADING  ! 
GOSUB  heading:  ! 

FOR  J«l  TO  It  ! 

GOSUB  COIL  DATA  ! 
GOSUB  COIL"  ! 

IF  evaporatorl=-air"  THEN  GOSUB  DX  EVAP    ! 
IF  evaporator*='«ater'  THEN  SOSUB  HATER  EVAP 
Td=95  ! 

I 

INTESATION: 

GOSUB  UNIT  MODEL  ! 

GOSUB  PLR  "! 

Qc=Qe+P/3.514  ! 

IF  condenser^'air"  THEN  GOSUB  AIR  COND    ! 

IF  condensers-Hater"  THEN  GOSUB  KATES  COND     ! 

IF  ABS  (Tc2-Tc!K.l  THEN  GOTO  CHECK     !" 

Tc!=Tc2 

SOTO  INTERATIQN  ! 

CHECK* 

IF  condenserJ="air'  THEN  GOTO  570 

IF  Tc  in=Tc  tin  THEN  SOSUB  FORCED  Tci  ! 

GOSUB'QUTPUT  ! 

NEXT  J  ! 

GOSUB  output:  ! 


END 


input  eguip«ent  constants 
input  all  coil  loading  conditions 
svstei  heading 
data  heading 

loop  for  eacn  systei  loading 

individual  coil  data 
find  cooling  load 
finds  Tel=Te 
finds  Tel=Tc» 
initial  estiiate  of  Tc! 


finds  full  load  capacity, poller 
finds  partial  load  poiier 
calculation  of  condenser  load 
find  lien  Tel  (air-cooled  condenser) 
find  nei  Tc!  (cooling  to»er! 
check  for  convergence 

continue  interation  of  Tel 


find  reduced  cooling  toner  cfa— if  necessary 
systei  steady  state  perfor»ance  output 

next  load  condition  for  processing 
annual  load  out out 


THEN  SOTO  net  bulb 
THEN  GOTO  grains 


COIL: 

IF  UNITGJ="KB 

IF  UNITSS'S- 

■et  bulb:  ! 

Tnb=nb  »i 

Tesp=Tnb 

GOSUB  HVTA 

Hin=H 

Hi=IHin-.24»Tdb)/(1061+.444fTdb) 


inlet  enthalpy  froi  net-bulb 


inlet  air  coisture  content 


inlet  enthalpy  froa  aoisture  content 


checking  tor  dehuaidification 


finding  enthalpy  for  no  dehuaidification 


enthalpy  for  dehuaidification 


total  load  in  Btuh 
sensible  load  in  Btuh 
latent   load  in  Btuh 
load  in  Ions 
load  in  Tons 
load  in  Tons 


convert  deg.F  to  deg.C 
convert  deg.F  to  deg.C 


capacity  in  kits 
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750  GOTO  dehuaid.ck 

770  grains:  ! 

730  8i=»b  Ki/7000 

790  BOSUB'KETbulb 

800  Hin=.  24*T(fb+Ki  *  ( 1061+.  444*Tdh) 

810  BQTQ  dehuaid  ck 

930  dehuaid  ck:  T 

840  Teag=Ts~ 

850  SOSUB  HVTfl 

860  IF  H5>Ui  THEN  SOTO  no  dehuaid  ELSE  SOTD  dehuaid 

380  no  dehuaid: 

890  KS=»i 

TOO  Hout=.24»TstKitll04IM44*Ts) 

710  SOTO  loading 

930  dehuaid: 

940  Hout=H  ! 

950  SOTO  loading 

970  loading: 

980  IF  Tdb<Ts  THEN  Hout-Hin 

990  3e=CFI(e*(Hin-Hout)t4.5  ! 

1000  Qlat=(lli-«S)KI041t.444«Tdb!»CFNe*4.5  ! 

1010  Qsens=Se-fllat  ! 

1020  "sens=asens/!2000  ! 

1030  "lat=01at/12000  ! 

1040  8e=Qe/12000  ! 

1050  RETURN 

1060  ! 

1070  ! 

1080  UNIT  MODEL: 

1090  Tcl=TTcl-32)*5/9  ! 

1100  Tel=(Tel-32)»5/9  ! 

1110  Teg=(Tcl-Td  DES)/Tratio-(Tel-Tel  SES1 

1120  «DCR=AUA2»Tiq+A3*Teg'2 

1130  CAP=AOCR*CAP  SE5  ! 

1140  C0PC=(Tcl+273.2) / (Tel-Tell 

1150  C0PCD=(Tcl  0ES+273.2)/(Tcl  DES-Tel  OES) 

1160  C0PCN=C0PC7CaPC3 

1170  C0PflN=Sl+B2»C0PCN+B3»C0PCN-2 

1180  COPflD=CflP  OES/FIP  OES 

1190  C0Pfl=C0Pflff*C0PflD  " 

1200  FlP=CAP/C0PA  ! 

1210  Tcl=Tcl*9/5+32  ! 

1220  Tel=Tel*9/5+32  ! 

1230  CAP=CAP/3.5I6  ! 

1240  RETURN 

1250  ! 

1260  ! 

1270  PLR: 

1280  IF  condenser^'nater"  THEN  SOTD  1310 

1290  CAP=CAP/subd*subcooI  ! 

1300  ! 

1310  PLR=Qe/CAP 

1320  IF  PLR>CYCLE  THEN  SOTO  1350 

1330  FFL=PLR»C4  ! 

1340  SOTO  1360 

1350  FFL=C1+C2«PLRK3»PLRA2  ! 

1360  P=FFUFLP 

1370  RETURN  , 

1380  ! 

1390  ! 

1400  DI  EVAP:  ! 

1410  Te?40 

1420  fpa=CFNe/AREA 

1430  CF=el+e2*fpa+e3«fpaA2 

1440  E';=fle-AREA»CF«(El+E2»T«btE3»T«hA2tE4»Te+E5«Te«Tiib+E6tTe*T»b"2tE7tTeA2tE8»T»b«Te"2tE9»T»b"2»Te',2) 

1450   0ERIV=-(AREAtCF*(E4+E5»T»b+E6tTiib*2t2tE7*Te+2tEatT»b'Te+2»E9»Ti«b-2*Te)) 

1460   Te=Te-EV/DERIV 

1470  IF  ABS  !EV/DERIV1>.001  THEN  SOTO  1440 

1480  Tel=Te 

1490  RETURN 

1500  ! 

1510  ! 

1520  WATER  EVAP:  !  finds  chilled  nater  teag. 

1530  IF  »tf,cntrU="return"  THEN  SOTO  1560 

1540  Tc*=Tcntrl_teap  !  supply  »ater  teap  contra! 


poller  in  kns 
convert  deg.C  to  deg 
convert  deg.C  to  deg 


convert  kits  to  tons  (CAP) 


correcting  for  subcooling  of  air  cooled  condenser 
(already  included  in  Hater  cooled  condenser  systeasl 


cycling  operation 
continuous  operation 


find  evap.  teap.  (Nenton-Ralphsonl 


1550  SOTO  1570 

1540  Tca=Tcntr!  teip-Qe*24/6Plle  !  return  alter  te»o  control 

1570  Tel=Tca  !  "  chilled  alter  teno 

1580  RETURN 

1590  ! 

1600  ! 

1610  AIR  COND:  !  finds  condensing  teap.  and  subcaoling 

1620  DTc=(Qc-Fl]/F2 

1630  subcooI=Dl+D2»DTctD3tDTcA2 

1640  Tc=TDB+DTc 

1650  IF  Tc<Tc  «in  THEN  Tc=Tc  «in  !  holding  tin  cond.  teip 

1660  Tc2=Tc  " 

1670  RETURN 

1660  ! 

16W  ! 

1700  MATER  COND:  !  finds  condensing  alter  tenp. 

1710  Tco=Tcl 

1720  60SUB  FREE  Tci 

1730  IF  TciUc  iin  THEN  Tc  in=Tc  iin  aSE  Tc  in=Tci  !  holding  tin.  alter  teip 

1740  Tc  out=Tc~in+9,c»24/6Pf!c  !    "  nea  leiving  alter  te«p 

1750  Tc7=Tc  ouF. 

1760  RETURN" 

1770  ! 

1780  ! 

1790  CD0LIN6  TWl  •   ' 

1300  ! 

1310  FREE  Tci:  !  floating  leaving  aatar  tesn 

1320  no  ansaer=0 

1830  Cfflt=CFH  tar 

1340  R=0 

1850  INTERJ="initiar 

1860  Rlg=8Prtc*8.33/(CFH  tar*. 07121 

1870  NTUA=C»(RlgA(-U+!<8))) 

1880  initial: 

1890  MM 

1900  SOTO  1930 

1910  final: 

1920  R=<RI+Rs)/2 

1930  A=Tco-8-TaB 

1940  SOSUB  COUNT 

1950  HTU=CT/Rlg 

1960  IF  ABS  (NTUA-NTU)/NTUA<.001  THEN  SOTO  2020 

1970    IF  NTUA>NTU  THEN  Rs=R 

1980  IF  NTUA<NTU  THEN  R1=R 

1990  IF  NTUAWU  THEN  INTER»=,fiml" 

2000  IF  INTER$='finir  THEN  SOTO  fiml 

2010  SOTO  initiil 

2020  Tci=Tco-R 

2030  RETURN 

2040  ! 

2050  FORCED  Tci:  !  fixed  leaving  alter  teap 

2060  R=Tco-Tc  in 

2070  ft=Tc  in-TKB 

2080  CFNt^CFH  tar 

2090  CFHlFCFh"  tar 

2100  INTER*=*fnitialIy" 

2110  SOTO  2190 

2120  initially: 

2130  IF  CFHt/CFN  tar>.3  THEN  CFNt=CFNt-.UCF«  tar  ELSE  CF«t=CFNt-.01*CFH  tar 

2140  IF  CFMt=0  THEN  no  ansaer=l 

2150  IF  CFNt=0  THEN  SOTO  2290 

2160  SOTO  2190 

2170  finally: 

2130  CFnXCFtlh+CFHll/2 

2190  Rlg=6PNc*8.33/(.0712»CF11t) 

2200  MTOA=C*(RlQ'(-(l+XN))) 

2210  60SUB  COUNT 

2220  NTU=CT/Rlg 

2230  IF  8BS  (INTUA-NTU1/NTUAK.001  THEN  SOTO  2290 

2240    IF  NTUA>NTU  THEN  CF«h=CFf1t 

2250    IF  NTUA<NTU  THEN  CFHl=CFf!t 

2260    IF  NTUA<NTU  THEN  INTER*="f inally- 

2270      IF  INTER*="finally"  THEN  SOTO  finally 

2280  SOTO  initially 

2290  RETURN 


2300  ! 
2310  ! 

2320  COUNT:  !  cooling  toner  subroutine 

2330  CT=0 

2340  Y(l)».l 

2350  Y(2)=.3 

2340  Y(3)=.2 

2370  Y(41=.3 

2330  Te«p=TMB 

2390  BQ5UB  HVTfl 

2400  Hft=H 

2410  TIKHB+A 

2420  FOR  1=1  TO  4 

2430   T«=TK+Y(I)*R 

2440   Hfl=Hfl+Rlq*Y!l)*fi 

2450   Te«p=T« 

2460   BQSUB  HVTA 

2470   HT=H 

2480   DH=HT-HA 

2490   CT=CT*1/0H 

2500  NEXT  I 

2510  CT=CT»R/4 

2520  RETURN 

2530  ! 

2540  ! 

2550  HVTfl:  !  enthalpy  subroutine 

2540  PB=101325 

2570  SOSUB  PSAT 

2580  IS=.622MPS/(PB-PS1) 

2590  H=.24»Te«ptHS*(1061t.444«Teip) 

2600  RETURN 

2410  ! 

2420  ! 

2430  HETbulb:  !  net-bulb  subroutine 

2440  leap  H=Tdb 

2450  Teip"L=0 

2460  h=.2?tTdb+Hi.(1061+.444Hdb) 

2670  PB=101325 

2680  Teio=(Teto  H+Tetp  Ll/2 

2690  SOSUB  PSflT" 

2700  T«b=Teip 

2710  «S=.622*(PS/(PB-PS1) 

2720  H".24»T«WB»<MM*.4+W»W 

2730  IF  (IBS  (H-HK.005  THEN  SOTO  2770 

2740  IF  H>h  THEN  Teip  H=Te«p 

2750  IF  H<h  THEN  Teip~l=TeiP 

2760  SOTO  2680 

2770  RETURN 

27S0  ! 

2790  ! 

2900  PSflT:  !  saturated  pressure  subroutine 

2S10  Ml. 11=6. 41533947 

2820  A(l,2)=6.41542599 

2830  All, 31=6.44302266 

2840  All, 41=6.53947833 

2850  A(2,l)=. 082103347S 

2860  AI2, 21=. 0724023398 

2870  AI2, 31=. 0698803796 

2880  fl!2,41=.06567255B4 

2890  AI3, 11=-. 000340068991 

2900  AI3, 2)=-. 000264361973 

2910  A(3,31=-2.0606744496E-4 

2920  AI3, 4)=-. 000160854949 

2930  TS=(Teip-32H5/9 

2940  IF  TS>22.0322  THEN  SOTO  2950  ELSE  SOTO  2990 

2950   IF  TS<46.0304  THEN  J=3  ELSE  SOTO  2970 

2960   SOTO  3050 

2970   IF  TS<72.0476  THEN  J=4  ELSE  SOTO  3030 

2980   SOTO  3050 

2990  IF  TS>. 00895  THEN  J=2  ELSE  SOTO  3010 

3000   60T0  3050 

3010   IF  TS>-19.94  THEN  J=l  ELSE  6OT0  3030 

3020   SOTO  3050 

3030  PRINT  ■         TEMPERATURE  IS  OUTSIDE  OF  OUR  RAN6E- 

3040  PAUSE 


3050  PS=€XP  (fl(l,J)+TSt(fl(2,J!+TS*fl(3,J))) 

3040  RETURN 

3070  ! 

3080  ! 

3090  HEADIN6: 

3100  PRINT  CHRJ  (12) 

3110  PRINT  '  SYSTEfl" 

3120  PRINT  ■  


90 


3130  IF  SYSTEMS  THEN  PRINT  USINB  "49X,75fl"  ;  "COMPRESSOR 01  COIL REMOTE  AIR-COOLED  CONDENSER' 

JIM  IF  SYSTEMS  THEN  PRINT  USIN6  "341,75ft-  ;  ■  WATER  COOLED  CONDENSING  UNIT--OX  COIL COOLING  TOWER 

3150  IF  SYSTElt=3  THEN  PRINT  USINB  -34X,75fl-  ;  ■  COMPRESSOR  CHILLER  UNIT WATER  COIL REMOTE  AIR-COOL 

ED  CONDENSER- 

3160  IF  SYSTEHM  THEN  PRINT  USING  ■47X.75A"  ;  'CHILLER  UNIT WATER  COIL COOLING  TONER' 

3170  IF  evaporatar*='\iater'  THEN  PRINT  USIN6  "44X,14A,6A,15fl,DD.D,6A'  ;   'chilled  Hater  ,,rtr.atrtt,'  teip 

.  held  at  ".Tcntrl  teip,'  deq.F" 

3180  PRINT 

3190  PRINT 

3200  RETURN 

3210  ! 

3220  ! 

3230  OUTPUT: 

3240  ' 

3250  (*•"  "*  ' 

3260  F*="  »»*«  ' 

3270  S*='   ♦***♦■ 

32S0  IF  condenser*='iiater'  THEN  GOTO  3310 

3290  PRINT  USING  '30.1,5(40,0) .6D.3I5D.D) ,3(4D.D) ,21,3(6(1) ,3A,2(aA> ,X,3A"  ;  x,TDB,TWB,Tdb,Ti<b,Ts,CFMe,Te!, 

Tcl,PLRHOO,Be,P,ac,D$,D*,D$,Ft,Fi,F*,Dl 

3306  SOTO  33a0 

3310  P  Cl=(Tco-Tc  in)/(Tco-Tci)*100 

3320  P"C2=CFMt/CFH  twHOO 

3330  IF  no  ans«er=0*  THEN  SOTO  3360 

3340  PRINT'USING  '3D,X,5(4D.D) ,6D,3(5D.B),3(4D.D! ,3(4D.D) ,5D.D,X,3(9A)'  !  x,TDB,TWB,Tdb,T»b,Ts,CFMe,Tel,Tc 

l,PLH*100,Qe,P,Qc,Tco,Tci,Tc  in.P  Cl,S$,Bt,F$ 

3350  GOTO  3380  "   " 

3360  PRINT  USING  "3D, I, 5I4D. D) ,6D,3(5D.D) ,6(4D.D), 50.0,1,2(801 ,6D.D"  I  x ,TDB,THB,Tdb .Tab ,Ts,CF!fe,Tal , Tel ,P 

LRM0O,ae,P,ac,Tco.Tci,Tc  in,P  C1,CFM  tw,CFNt,P  C2 

3370  !  -       -  - 

3380  IF  x=l  THEN  LOAD  annual=0  !  annual  variables 

3390  IF  1=1  THEN  Qseni  annual=0 

3400  IF  x=l  THEN  Olat  annual=0 

34t0  IF  x=l  THEN  hrs  annual^ 

3420  IF  x=l  THEN  P  annual=0 

3430  Qsens(x)=8sens*l2 

3440  Qlat(x)=Qlat»12 

3450  Be(x)=8e»12 

3460  P(x)=P 

3470  Osens  annual=8sens  annual+Ssens(x)*hrs(x) 

3480  aiat .annual=81at .annual +31at(x)*hrs(x) 

3490  hrs  annual=brs  annual+hrslx) 

3500  LOAD"  annual  =LOBD  annual*0e!x)»hrs(x. 

3510  P  annual^  annuaT>P»hrs(x) 

3520  RETURN 

3530  ! 

3540  0UTPUT2: 

3550  PRINT 

3560  PRINT 

ANNU 


3570  PRINT  ' 

AL' 

3530  PRINT  ' 

CASE 
1 

STEADY  STATE 

Osensible  Blatant  Qtotal 
(Hbtuh)  

Ol/Qt 

I 

POWER 

llutsl 

3590  PRINT  • 

poiier/yr  " 
3600  PRINT  ! 

hours 
(hrs) 

I    (kx-hrs)    I" 
3610  PRINT  " 

3620  FOR  >1  TO  X 

3630  P  C4=81at(I)/8e<I)HO0 

3440  a5=ae!I)*tirs(I) 

3650  P  C5=85/L0AD  annualHOO 

3660  P5=P(I!thr5(D 

3670  ?  C6=P6/P  annual* 100 

laad/yr 
iMBTU) 


3680  PRINT  USIN6  " 19X , 4D,?D. 0,7D. D, 7D.D,fiD.D,6D.D, 10D, 1 1D,5D. D, I 1D,7TJ.D" 
tI),hrs(I),a5,P  C5,Pi,P  C4 
3690  NEXT  {   " 


I,8sens(I),«at(I),ae(I),P.C4,P 


3700 

3710  PRINT  USING  '3/' 

3720  PRINT  USIN6  "30X,34fl,5D. 

3730  PRINT  USING  "30X,34A,5D. 

3740  PRINT  USING  "30X,34A,5D. 

3750  PRINT 

3740  PRINT  USING 

3770  PRINT  USING 

3780  RETURN 

3790  ! 

3800  ! 

3810  HEA0IN62: 

3320  PRINT  ' 

3830  PRINT  ■ 


1,18ft" 

],13A" 


30X,34fl,aD,iafl* 

30X,34A,8D,19A" 


'Annual  latent  cooling  load 
'Annual  sensible  cooling  load 
"Annual  total   cooling  load 

'Annual  hours  of  operation 
"Annual  Power  Consuaption 


=  ";Qlat_annual/!000, 
■  "jQsens  annual/ 1000 
=  ■ j  LQflD.Innual / 1000, 


";hrs_annual, 
■;P_annuaI,' 


SKBtu" 
RltBtu' 
URStu' 


1      Hours" 
Kw-hrs" 


SYSTEM  STEADY  STATE  OPERATION 


3840  PRINT  '    .aabient 

...  cool  ing  toner " 

3850  PRINT  'CASE      TUB      TKB      Tdb 
c  in      I      CF»  tar     CFtlt      X" 

3340  PRINT  '  t  " deg.F... 

(flturl   ...cfi...   (CFHtwr) 
3870  PRINT  '— 


...coil.. 
Tub   Ts 


CFH 

eta 


Tel   Tel 
..deg.F.. 


.systea 

PLR   Se 
I   tons 


P 

Sews 


9c 
tons 


o   Tci 
.deg.F. 


3880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3940 
3970 
3980 
3990 
4000 
4010 
4020 
4030 
4040 
4050 
4040 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4140 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4240 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4340 


RETURN 


CONSTANTS: 

ON  SYSTEM  GOTO  UNIT1  ,UNIT2  ,UNIT3  ,UNIT4 


UNIT1;  ! 

condenser*='air"  ! 

evaporatorl="air" 

Tratio=2.9 

Al=.9973 

IB"-.  03494 

O*.  0004335 

Bl=-.4141 

67*1.7717 

B3=-.3494 

Cl=.1454 

C2=.9555 

C3=-.1048 

Tel  DES=(105-32)»5/9 

Tel"DES=(40-32H5/9  ! 


(coapressor  subsytea)  systea 
constants  good  for  all  unit!  systeis- 


teap  at  which  CAP  0ES,RP  DES  selected 
teap  at  which  CAPl)ES,FLP"J}ES  selected 


lines  3950  to  4080 


UNIT  CONSTANTS 

CAP_0ES=17.3»3.514  T     coapressor  capacity  in  kws  (  froa  lach,  data  1 

FLP  DE5*15.9  !         coapressor  power   in  kws  (  froa  aach.  data  ) 

CYCL"E=.25  !  PLR  at  which  systea  cycles 

C4=(C1+C2»CYCLE+C3»CYCLE*2) /CYCLE  !        systea  cycling  constant 

!  DX-COIL  CONSTANTS 

GOSUB  DX  CONSTANTS  !  "  - 
AREA=5  !" 


subcl=1.075 
subcool=subd 
m*.9774  ! 
D2=. 0044485  ! 
03=-. 00002978 
F2=. 94203  ! 
FU-3.4047  ! 
Tc  ain=85  ! 
RETURN 


UNIT2:   ! 

condenser*="water" 

evaporatori="air" 

Tratio=2.9 

Al=.9982 


CONDENSES  CONSTANTS 


loading  OXJVAP  constants 

DX-EVAP  face  area  selected  by  designer (sg. ft.) 

jnit  base  on  14  deg.F  subcooling 
used  to  find  actual  subcooling 
constants  for  subcoolino 


load/dt  of  air-cooled  condenser  (tons/deg.F) 
offset  froa  zero  of  above  curve 
■iniaua  condensing  teaperature  (deg.F) 


(water  cooled  condensing  subsvstea) systea 
constants  good  for  ill  unit!  systeas— lines  4330  to  4440 


91 


4370 
4380 
4370 
4400 
4410 
4420 
4430 
4440 
4450 
4440 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4400 
4410 
4420 
4430 
4440 
4450 
4440 
4470 
4480 
4490 
4700 
4710 
4720 
4730 
4740 
4750 
4740 
4770 
4780 
4790 
4800 
4810 
4820 
4S30 
4840 
4850 
4840 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
3050 
5040 
5070 
3080 
5090 
5100 
5110 


A2=-.  03952 

O*.  004931 

81=-.2244 

82=1.5438 

B3=-.3219 

Cl=.129 

C2=.8849 

C3=-  0130? 

Tel  DES=(95-32)*5/9 

Tel"DES=(40-32)«5/9 
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teap  at  which  CAP  DES.FLP  DE5  selected 
teip  at  which  CAP~DES,FLPT)ES  selected 


UNIT  CONSTANTS 

CAPJE5=16. 8*3.516  T      coipressor  capacity  in  kws  (  froi  >ach.  data) 

FLP  DES=Si.i  !         compressor  power   in  kws  (  froi  »ach.  data) 

CYCCE=.25  !  PLR  at  which  systei  cycles 

C4= (CI +C2*CYCLE+C3tCYCLEA2) /CYCLE  !  systei  cycling  constant 

!  DX-EVAP  CONSTANTS 

SOSUB  OX  CONSTANTS  !     loading  OJ-EVAP  constants 

AREA=5  !  BX-EVflP  face  area  selected  by  designer (sq.ft. ) 

!  C00L_TKR_C0NSTANTS 

C=l. 23144  !  Hebb-Villacres  cooling  tower  constant 

XN=.9193  !  Hebb-Vi 11  acres  cooling  tower  constant 

CF"._twr=4800  !  air  flow  rate  of  cooling  tower  (  cfi  ) 

BP".c=44.9  !  water  flow  rate  of  tower, condenser  I  gps  ) 

Tc  »in=75  !  liniiui  entering  condensing  water  tesp(deg.F) 

RETURN 


UNIT3:  ! 

condenser*="air'  ! 

evaporatorl='water" 

Tratio=2.9 

Al=.9982 

A2=-.0299 

A3=. 0001711 

Bl=-.2484 

B2=1.66 

83=-. 4105 

Cl=.1545 

C2=.8323 

C3=.  01672 

Tel  DES=(105-32)*5/9  ! 

Tel"DES=<46-32)»5/9  ! 


(coipressor  chiller  subsystei) systei 

constants  good  for  all  unit3  systeis— lines  4630  to  4810 


teip  at  which  CAP  OES.RP  DES  selected 
teip  at  which  CAPT)E5,FLPT)E5  selected 


UNIT  CONSTANTS 

CAP  J)ES=16. 1*3.516  T    coip-chiller  capacity  in  kws  (  fro«  lach.  data  I 

FLP  DES=15.6  !        coip-chiller  power  ■  '  in  kws  (  froi  tach.  data  ) 

CYCL"E=.25  !  PLR  at  which  systei  cycles 

C4=<C1+C2*CYCLE+C3»CYCLEA2) /CYCLE  !  systei  cycle  constant 


subcl=1.075 
subcool=subcl 
01=. 9774  ! 
02=.00446S5  ! 
03=-. 00002978 
F2=. 94203  ! 
Fl=-3.4047  ! 
Tc  «in=85  ! 


CONDENSER  CONSTANTS 


unit  based  of  16  deg.F  subcooling 

used  to  find  actual  subcoolina 

constants  for  subcoolinolal!  TRANE  condensers) 

load/dt  of  air-cooled  condenser  (tons/deg.F) 
offset  froi  zero  of  above  curve 
•iniiui  condensing  tetperature  (deg.F) 


wtr  cntrl*='supply" 
Tcnfrl  teip=45 
GPNe=3?.5  ! 
RETURN 


UNIT4:  ! 

condenserl="water' 

evaporator<="water" 

Tratio=2.9 

Al=.9978 

A2=-. 02943 


COIL  CONSTANTS 


chilled  water  teip  control— return  or  supply 
teip  at  which  above  water  is  held  at  (deg.F) 
chilled  water  flow  rate   (gpi) 


(chiller  substei)  systei 
constants  good  for  all  unit4  systeis— lines  5070  to  5200 


5120 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 
5260 
5270 
5230 
5270 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570 
5580 
5590 
5600 
5610 
5620 
5630 
5640 
5650 
5660 
5670 
5680 
5690 
5700 
5710 
5720 
5730 
5740 
5750 
5760 
5770 
5780 
5790 
5800 

saio 

5820 
5830 
5840 
5850 
5860 


A3*. 0001472 

B1=-.1U9 

B2=1.5815 

B3=-.4702 

Cl=. 1334 

C2=.743? 

C3=.1209 

Tc1'dES=(95-32)*5/7 

Tel"DES= (44-32) «S/7 
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teap  at  which  CAP  DES.FLP  DEB  selected 
teap  at  nhich  CAPlES, FIPTIS  selected 


!  UNIT  CONSTANTS 

CAP  DES=15. 8*3.516  T    chiller  capacity 

FLP"DES=14  !         chiller  Doner 

CYC[E=.25  ! 

C4=(Cl*C2fCYCLE+C3»CYCLEA2) /CYCLE  ! 


in  kns  (  froa  tach.  data  ) 
in  kns  (  fro*  aach.  data  ) 
PLR  at  which  systea  cycles 
systea  cycling  constant 


!          COIL  CONSTANTS 

ntr  cntrlt=*supply"  !    chilled  nater  teip  control— return  or  supply 

Tcnfr!  teip=45  !  teap  at  nhich  above  nater  is  held  at  Ideg.F) 

SPtle=34.5  !  chilled  nater  flon  rate  (gp«) 


C=t. 23146 
XN=.?173  ! 
CFK  tnr=4800  ! 
BPBc=44.8  ! 
Tc  iin=75  ! 
RETURN 

I 

DX  CONSTANTS!  ! 
El**. 08552 
E2=-. 0844875 
E3=. 000813752 
E4=-. 242335 
E5=.0O21154 
E6=.  0000 1460892 
E7«. 0018496 
EB=-. 00001131545 
E9=-. 000000323206 
el*. 0911429 
e2=. 002245425 
•3=-. 000000855422 
RETURN 


COOL  TKR  CONSTANTS 


Webb-Vi II acres  cooling  toner  constant 
Kebb-Villacres  cooling  toner  constant 
air  flon  rate  of  cooling  toner  (  cfi  ) 
nater  flon  rate  of  toner, condenser  I  gpa  1 
iiniaua  entering  condensing  nater  teip  (deg.F) 


constants  for  all  6-ron  AEROFIN  'net'  DX-coil  evaporators 


COIL  CONDITIONS: 


EI 


aibient 


coil  conditions 

inlet  outlet 


dry-bulb  , net-bulb,   dry-bulb, net-bulb  OR  grains, cfi  ,  drv-bulb 
DATA    95   ,   81   ,    95   ,   75, KB,       2500,   55 
DATA    95   ,   81   ,    95   ,         120,8,  2500,   55 


w 


READ  hrs(I) , TUB IX) .TKBIX) .Tdb(X) ,nb  ni (X) ,UNIT5J(X) ,CF«e(X) ,Ts(X) 
IF  TDB(X)=0  THEN  SOTO  5720 


X=X+! 
SOTO  5630 
W-l 

RETURN 
i 

DATA  7, 97, 79. 97, 77, SB. 2000, 55 
DATA  H2,92,78,92,78,NB,2000,55 
DATA  574,87.77.87.77, KB, 2000,55 
DATA  1008, 82, 75, 85,74, KB, 2000, 55 
DATA  1526, 77, 73, SO, 74, KB, 2000, 55 
DATA  1325, 72, 49,76, 72, KB, 1850, 55 
DATA  925, 67, 63, 74, 70, KB, 1611. 55 
DATA  77? , 42 , 5B , 70 , 68 , KB , 1 500  i  55 
DATA  701, 57, 53, 48, 64, KB, 1500, 55 
DATA  0,0, 0,0,0, KB, 0,0' 

COIL  DATA: 


5870  T0B=TBB(x)  n. 

5880  THB-THB(x)  94 

5890  TdtJ=Tdbtx> 

5900  nil  «i=»t]  Kill) 
5910  UNrTS»=U(JlTSt(x) 
5920  CFHe=CFHe<x) 
5930  Ts=Ts(x] 
5940  RETURN 


APPENDIX   D. 

TYPICAL  ANNUAL  PERFORMANCE  OUTPUT 
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APPENDIX   E 


RECIPROCATING  CHILLER— TWO  COIL- 
COOLING  TOWER  SYSTEM  COMPUTER  PROGRAM 


THIS  PROGRAM  IS  FOR  A  chiller-cooling  toner-tnin  chilled  nater  coil-svstei 

THE  CONSTANTS  DESCRIBING  EOUIPRENT  SIZE  ARE  ALREADY  INCLUDED  IN  THE  CONSTANT  SUBROUTINE 

CYLINDER  UNL0ADIN6  UTILIZED  TD  CONTROL  CHILLED  WATER  TENPERATURE 

CONSTANT  WATER  FLOW  RATES  ASSUMED  AT  THE  EVAPORATOR  4  CONDENSER— CHILLED  HATER  BYPASS 

THE  COIL  CONDITIONS  ARE  AVAILABLE  TO  THE  PROSRAN  THROUGH  THE  COIL  CONDITION 


THEN  GOTO  CHECK 


10  REH 
20  REH 

22  REN 

23  REH 
30  REH 
40  ! 
50  I 
40  ! 

80  DIN  Qsens<200) ,Qlat(200) ,0e(200) ,P(200) ,hrs(200) ,TDB(200) ,TWB(200) 

90  DIN  Tdb  1(200), nil  ni  I(200),CFNe  l(200l,Ts  1(200), UNITS  1*<200) 

100  DIN  Tdf/  2(200), nS  nf  2(200), CFHe  2(2001, Ts  21200)  .UNITS'  21(200) 

110  DIN  COIL"  1(2001, CdIL'21200)         " 

120  ! 

150  ! 

160  SOSUS  CONSTANTS  ! 

170  S0SU8  COIL  CONDITIONS 

ISO  SOSLIB  HEADING  ! 

190  SOSUB  HEADING2  ! 

200  ! 

210  FOR  x=l  TO  J  ! 

220  ! 

230  60SUB  COIL  ! 

240  SOSUB  WATER  EVAP  ! 

250  Td=95  ! 

260  ! 

270  ITERATION: 

230  GOSUB  UNIT  HOEEL  ! 

290  GOSUB  PLR  "! 

300  Qc=0e+P/3.516  ! 

310  GOSUB  WATER  COND     ! 

320  IF  ABS  (Tc2rTcl)< 

330  Td«Ic2 

340  SOTO  INTERATION  ! 

350  ! 

360  CHEEK* 

370  IF  Tc'in=Tc  «in  THEN  SOSUB  FORCED  Tci 

380  60SUBHUTPUT  ! 

390  ! 

400  NEXT  x  ! 

410  SOSUB  0UTPUT3  ! 

420  SOSUB  0UTPUT2  ! 

430  END 

440  ! 

450  ! 

440  COIL: 

470  COIL  1=0  ! 

480  COIL  2=0 

490  Se=0~ 

500  Qlat=0 

510  Qsens=0 

520  ! 

530  FOR  DATA  coil-1  TO  2  ! 

540  SOSUB  COIL  DATA  ! 

550  IF  CFHe=0  THEN  60TO  950  ! 

540  ! 

570  IF  UNITSI="»B-  THEN  GOTO  net  bulb 

580  IF  UNITSJ="6-  THEN  SOTO  grains 

400  Mt  bulb:  ! 

410  Tnb=»b  Hi 

620  Teno=T»b 

430  SOSUB  HVTA 

440  Hin=H 

450  Ni=(Hin-.24*Tdb)/(1041+.444fTdb)  ! 

440  SOTO  dehuiid  ck 

400  grains: 

690  8i=nb  wi/7000 

700  Tnb=Teip 

710  Hin=.24*Tdb+Wi»(1061t.444*Tdb)  ! 

720  GOTO  dehuiid  ck 

740  dehuiid  ck:  T 

750  Teip=Ts~ 

760  SOSUB  HVTA 

770  IF  NS>«i  THEN  SOTO  no  dehuiid  ELSE  SOTO  dehuiid 

790  no  dehuiid: 


input  equipient  constants 
input  all  coil  loading  conditions 
systei  heading 
data  heading 

loop  for  each  systei  loading 

finds  cooling  load 

finds  Tel=Tcn 

initial  estisate  of  Tel 


finds  full  load  capacity, poner 
finds  partial  load  power 
calculation  of  condenser  load 
find  new  Tel  (cooling  toner) 
check  for  convergence 

continue  interation  of  Tci 


find  reduced  cooling  toner  cf» — if  necessary 
systei  steady  state  perforaance  output 

next  load  condition  for  processing 
coil  condition  output 
annual  load  output 


initializing  variables 


tno  coil  systei 
data  for  each  coil 
for  coil  bypass 


inlet  enthalpy  froi  net-bulb 


inlet  air  loisture  content 


inlet  air  enthalpy  froi  loisture  content 
checking  for  dehuiidification 
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finding  enthalpy  for  no  dehuaidif ication 


enthalpy  for  dehu«idification 


total   load  in  Btuh 
sensible  load  in  Btuh 
latent   load  in  Btuh 
II  coil  load 
<2  coil  load 


load 
load 
load 
load 


n  Tons 
n  Tons 
n  Tons 
n  Tons 


300  US-Hi 

810  Hout=.24*Ts*HiM1061+.444»Ts] 

820  6QTD  loading 

340  dehuiid: 

850  Hout=H  ! 

860  SOTO  loading 

3S0  loading: 

890  IF  Tdb(Ts  THEN  Hout=Hin 

900  Qe=fle+CFBe«IHin-Hout)M,5  ! 

910  Qlat=giat*(Ni-HS)*(1061*.444Hdb)*CFHe*4.5 

920  Qsens=5e-01at  ! 

930  IF  DATA  coil=l  THEN  COIL  l=Qe  ! 

940  IF  DATA~coil=2  THEN  C0IL~2=0e-C0IL  1  ! 

950  NEXT  DATA  coil 

940  ! 

970  Qsens=dsen5/ 12000  ! 

980  Qlat=Qlat/I2000  ! 

990  Qe=fle/12000  ! 

iooo  coil  wan.  1/12000  i 

1010  C0IL~2=C0IL~2/t2000  ! 

1020  RETURN 

1030  ! 

1040  ! 

1050  UNIT  H0DEL: 

1060  Td=TTcl-32)«5/9  ! 

1070  Tel=(Tei-32!»5/9  ! 

1080  Teg=(Td-Td  DE5)/Tratio-(Tel-Tel  DES) 

1090  A0CR=Al+A2*TeVA3»TegA2 

1100  CAP=ADCR*CAP  DES  ! 

1110  C0PC=<Tc!*27T.2)/(Tcl-Te!) 

1120  C0PCD=(Tcl  DES+273.2)/ITcl  DES-Tel  OES) 

1130  C0PCN=C0PC7C0PCD 

1140  COPAN=81+B2*C0PCN+B3»COPCNA2 

1150  C0PAD=CAP  DES/FLP  OES 

1160  C0PA=C0PAR*C0PAD  " 

1170  RP=CAP/C0PA  ! 

1130  Td=TcU9/5+32  ! 

1190  Tel=TeU9/5*32  ! 

1200  CAP=CAP/3.516  ! 

1210  RETURN 

1220  ! 

1230  ! 

1240  PLR: 

1250  PLR=8e/CAP 

1260  IF  PLR>CYCLE  THEN  SOTO  1290 

1270  FFL=PLR*C4  ! 

1280  SOTO  1300 

1290  FFL=C1*C2»PLR+C3«PLR"2  ! 

1300  P=FFL»FLP 

1310  RETURN 

1320  ! 

1330  ! 

1340  HATER  EVAP: 

1350  IF  rtr  cntrlJ="return'  THEN  SOTO  1330 

1360  Tc»=Tcntr!  tetp  ! 

1370  SOTO  1390  " 

1380  Tc»=Tcntrl  teip-Qe*24/6Plte  ! 

1390  Tel=Tc»  !  ' 

1400  RETURN 

1410  ! 

1420  .' 

1430  HATER  COND: 

1440  Tco=Tcl 

1450  SOSUB  FREE  Tci  ! 

1460  IF  Tci<Tc  iin  THEN  Tc  in=Tc  iin  ELSE  Tc  in=Tci 

1470  Tc  out=Tc"in+Qc»24/6PHc  !  " 

1480  Tc2=Tc  ouf. 

1490  RETURN" 

1500  ! 

1510  ! 

1520  C00LIN6  THR:  !  SLIBHTLY  MODIFIED  HEBB-VILLACRES  COCilNS  TOHER 

1530  ! 

1540  FREE_Tci:   !  floating  leaving  nater  te«p. 

1550  no  ans*er=0 

1560  CFn*t=CFH  tw 
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load  in  Tons 


convert  deg.F  to  deg.C 
convert  deg.F  to  deg.C 


capacity  in  k«s 


poiier  in  k»s 
convert  deg.C  to  deg.F 
convert  deg.C  to  deg.F 
convert  k»s  to  tons  (CAP) 


cycling  operation 
continuous  operation 


supply  »ater  tent  control 

return  Hater  teip  control 
chilled  xater  teip 


cooling  toner  for  Tci 
holding  iin.  Mater  teip 
nen  leaving  natsr  teip 


1570  8=0  , ,, 

1580  !NTERl="initiar  '  '  ' 

1590  Rlg=SPflc*a.33/(CFH  Urf.0712) 

1600  ™=C«(Rlg'(-!l+X(J))) 

1610  initial: 

1620  R=R+2 

1630  SOTO  1660 

1640  final: 

1650  R=(Rl+Rs)/2 

1660  A=Tco-R-T«B 

1670  60SUB  COUNT 

1680  NTU=CT/Rlg 

1690  IF  ABS  (NTUA-NTU)/NTUA<.001  THEN  60TD  1750 

1700    IF  NTUA>NTU  THEN  Rs=R 

1710     IF  NTUA<NTU  THEN  R1=R 

1720     IF  NTUA<NTU  THEN  IMTER*="f inal " 

1730      IF  INTERS="final'  THEN  SOTO  final 

1740  SOTS  initial 

1750  Tci=Tco-R 

1760  RETURN 

1770  ! 

1780  FORCED  Tci:  !  fixed  leaving  nater  te«p 

1790  R=Tco-Tc  in 

1800  A=Tc  in-TNB 

1810  CFHt=CFH  t»r 

1810  CFHh=CFh"tw 

1830  INTER^Tnitially" 

1840  SOTO  1920 

1850  initially: 

1860  IF  CFNt.'CFN  t»r>.3  THEN  CFHt=CF!lt-.UCF!(  t»r  ELSE  CFNt=CF«t-.01»CFN  twr 

1870  IF  CFNt=0  THEN  no  ans«er=l 

1880  IF  CFHt=0  THEN  SOTO  2020 

1890  SOTO  1920 

1900  finally: 

1910  CFHt=fCFhVCFM)/2 

1920  Rlg=6P«c*8.33/1.0712*CFMt) 

1930  NTOA=CKRlgA(-(l+XN>)) 

1940  SOSUB  COUNT 

1950  NTU=CT/Rlg 

1960    IF  ABS  ((NTUA-NTUl/NTUAK.OOl  THEN  SOTO  2020 

1970         IF  NTUA>NTU  THEN  CFHh=CFHt 

1980         IF  NTUA<NTU  THEN  CF«l=CFHt 

1990         IF  NTUA<NTU  THEN  INTERt="f inally" 

2000  IF  INTEM="finally'  THEN  SOTO  finally 

2010  SOTO  initially 

2020  RETURN 

2030  ! 

2040  ! 

2050  COUNT:  !  win  cooling  toner  subroutine 

2060  CT=0 

2070  Y(l)=.l 

2080  Y(2)=.3 

2090  Y(3)=.2 

2100  Y(4)=.3 

2110  Te«p=TNB 

2120  SOSUB  HVTA 

2130  HA=H 

2140  T1KHB+A 

2150  FOR  1=1  TO  4 

2160   TK=TN+Y(I)*R 

2170   HA=HAtRlg*Y(I)*R 

21S0   Te«p=T« 

2190   SOSUB  HVTA 

2200   HT=H 

2210   DH=HT-HA 

2220   CT=CT+1/BH 

2230  NEXT  I 

2240  CT=CT*R/4 

2250  RETURN 

2260  ! 

2270  ! 

2280  HVTA:  !  enthalpy  subroutine 

2290  PB=101325 

2300  SOSUB  PSAT 

2310  HS=.622KPS/!P8-PS1) 


2320  H=.24*Teip**S»!1061+.444fTe«p)  119 

2330  RETURN  '  ' c 

2340  ! 

n350  ' 

2360  HETbulb:  !  net-bulb  subroutine 

2370  Te«p  H=Tdb 

2380  Te«p~L=0 

2390  h=.2?*Tdb+«i*{10iI+.444*Tdb) 

2400  PB=101325 

2410  Te»p=(Teip  H+Teip  D/2 

2420  GOSUB  PSAT" 

2430  T*b=Teip 

2440  HS=,622t(PS/tPB-PS)) 

2450  H=.24*Tnb+HS*(1061+.444*Titb) 

2440  IF  ABS  (H-DK.005  THEN  SOTO  2500 

2470  IF  H>h  THEN  Tetp  H=Teip 

2480  IF  H<h  THEN  Teip~L=Teip 

2490  SOTO  2410 

2500  RETURN 

2510  ! 

2520  ! 

2530  PSflT:  !  saturated  pressure  subroutine 

2540  Ml.  11=6.41533947 

2550  Ml, 21=6. 41542599 

2560  All, 31=4. 44302266 

2570  0(1,41=4.53947833 

2580  A(2,ll=.0821033478 

2590  A(2, 2)=. 0724023398 

2600  A(2,31=.0698B03796 

2410  A (2, 4) =.0454725534 

2420  MJ.l)"-.  00034006899! 

2630  MJ, 21=-. 000264361973 

2640  A(3.31=-2.0606744496E-4 

2650  A(3,41=-.000160854949 

2660  TS=(Te»n-32)«5/9 

2670  IF  TS>22.0322  THEN  SOTO  2680  ELSE  SOTO  2720 

2680   IF  TS<46.0304  THEN  J=3  ELSE  SOTO  2700 

2690   SOTO  2780 

2700   IF  TS<72.0476  THEN  J=4  ELSE  SOTO  2760 

2710   SOTO  2780 

2720  IF  TS>. 00895  THEN  J=2  ELSE  BOTO  2740 

2730   SOTO  2780 

2740   IF  TS>-19.94  THEN  J=l  ELSE  SOTO  2760 

2750   SOTO  2780 

2760  PRINT  "         TEMPERATURE  IS  OUTSIDE  OF  OUR  RANGE" 

2770  PAUSE 

2730  PS=EIP  IA(l,J)tTS*(A(2,J)+TS*A(3,Jl)) 

2790  RETURN 

2800  ! 

2310  ! 

2820  HEADIN6: 

2330  PRINT  CHR1  (12! 

2840  PRINT  ■  SYSTEM'  ' 

2350  PRINT  '  : 


2860  PRINT  USING  "451, 70A"  ;  "  CHILLER CHILLED  HATER  COILS COOLING  TONES" 

2870_PRINT  USING  '44X, I4A,6A, 15A,DD. D ,6A"  ;  -chilled  Hater  ■.utr.oitrlV  tetp.  held  at  \Tcntrl_te5p,"  d 

2880  PRINT 

2890  PRINT 

2900  RETURN 

2910  ! 

2920  ! 

2930  OUTPUT: 

2940  P  Cl=(Tco-Tc  inl/(Tco-Tci)»100 

2950  P"C2=CFHt/CFK  txrtlOO 

2960  IF  no  ans*er=0"  THEN  SOTO  2990 

2970  PRINT  USING  "5X ,3D,X,2 (4D.D) ,X,2!6D.D) ,I,2(4D.DD1 ,5D.D,4D.DD,4(5D.0r  ;  s,TDB,TSB,Tel,Tc[,Se,?,?LR»!0 

0,Sc,Tco,Tci,Tc  in,P  CI 

2980  SOTO  30 10  "   " 

2990  PRINT  USING  "5I,3D,I,2!4D.D1 ,I,2(6D.D) ,X,2(4D.DD),5D.0,4D.0D,4(5D.D) ,2!3D) ,6D.D"  1  s ,TBB,T»B,Tel,Tcl, 

Be,P,PLR*100,0c,Tco,Tci,Tc_in,P_cl,CF(l_tMr,CFr1t,P_C2 

3010  iF  x=l  THEN  LOAD_annual=0  !  initializing  annual  variables 

3020  IF  <=1  THEN  P  annual=0 


3030  IF  ««l  THEN  COIL  1  annual=0 

3040  IF  «■!  THEN  C0"IL"2"annual=0 

3050  IF  ««l  THEN  sens~annual=0 

3060  IF  ««!  THEN  lit  annual=0 

3070  IF  ««1  THEN  hrs~annual=0 

3080  COIL  l(j)«COU  TM2 

3090  COIL  2(x)=C0IL~2»12 

3100  flseni(x)=Qsen5?12 

3110  Qlat(x)=Qlat*12 

3120  Qe(x)=9e»12 

3130  P(x)=P 

3140  hrs  annual=hrs  annual+hrslx) 

3150  LOAl  annual =UDSD  annual+Qe(x)*hrs(x) 

3160  P  annual=P  annuaT+PHirsIx) 

3170  CDIL  1  annual=CQIL  1  annual  +COIL  l(x)»hrs(x> 

3180  C0IL"2"annual=C0IL"2"annuaI tC0IL"2(x ) »hrs(x ) 

3190  sens  annual=sens  annual+Qsens(x)Jhrs(x) 

3200  lat  annual=lat  annual+61at(x)*hr5(x] 

3210  RETURN 

3220  ! 

3230  0UTPUT2: 

3240  PRINT  USIN6  "2/" 

3250  PRINT  ■ 

ANNUAL' 
3260  PRINT  '  
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STEADY  STATE 


3270  PRINT 
!oad/yr 
3280  PRINT 
It)   I 
3290  PRINT 


'CASE  TDB  TK8 
po»er/yr 
1  4   .deg.F, 
(few-hrs)    V 


BUM  CQILI2  Qsensible  Qlatsnt  Qtotal   Bl/Qt   POWER 


.  (Kbtutt). 


!kss) 


hours 
!hrs) 


(BBT 


3300  FDR  1=1  TO  I 

3310  lpc*«lat(I)/Qe(I)«100 

3320  aa=ae(I)thrs(I) 

3330  Qoc=3a/L0A8  annualMOO 

3340  Pa=P(I)thrsTll 

3350  Ppc=Pa/P  annualHOO 

3360  B.  1=CDIC  1(11 

3370  CL"2=CH1L"2(I1 

3380  PRTNT  USING  "4DJ,2(4D1 ,5(7D.D),2(6D.D),10D,11D,4D.DD,11D,6D.DD 

.aiat(I).Qe(I),Lpc,Pill,hrs(I),aa,flpc,Pa,Ppc 

3390  NEIT  I 

3400  PRINT  USIN6  "3/" 

3410  PRINT  USING  "30X,35A,5D.DD,18A" 

3420  PRINT  USING  '30)! , 35A.5D.DD, lafl- 

3430  PRINT 

3440  PRINT  USIN6  "30X,35A,5D.DD,18A" 

3450  PRINT  USIN6  "JOI,3M,».80,1M" 

3460  PRINT 

3470  PRINT  USING  "301 .35A.5D.0D, 18ft" 

3480  PRINT  USING  '30X.33A, SUM'  i 

3490  PRINT  USING  "J0I,3H,HUM' 

3500  RETURN 

3510  ! 

3520  ! 

3530  HEADIN62: 

3540  PRINT  " 

3550  PRINT  "  


I,TuB(I),T«BU]Ia._:,a._2,SsstisU) 


'Annual 
'Annual 


cooling  coil  I   1 
cooling  coil  t  2 


load 
load 


Annual  sensible  cooling  load 
Annual  latent   cooling  load 

Annual  cooling  load 
Annual  Hours  of  operation 


";COIl  1  annual/1000,' 
■;C0Ii;2>nual/!000,' 

>';uns  annual/1000." 
=";lat.annual/1000,' 

=';LDAD_anr.ual/1000,' 
;hrs_annual,"  Hours' 


fflBtu" 
IMBtu" 


mto 

IMBtu" 


NltBtii" 


"Annual  Chiller  Poner  consuaption  =';P  annual,"  Ka-hrs" 


SYSTEM  STEADY  STATE  OPERATION 


3560  PRINT  " 

oxer ' 

3570  PRINT  "    CASE 

CFN  t«r  CFflt  V 
3580  PffINT  "  1  .. 
r)  ...cfi...  (CrUtw) 
3590  PRINT  ' 


.aibient. 
TDB  TUB 


Tc» 
deg.F... 


Tc! 


..systei cooling  t 

Be    P    PLR   Qc    Tco  Tci   Tcjn   I 
tons   <»5    I   tons   ....deg.F (QtK 


3600  RETURN 

3610  ! 

3620  '. 

3630  0UTPUT3: 

3640  PRINT  USING  'If 

3650  PRINT  " 

3660  PRINT  ' 


COIL  CONDITIONS" 


equipaent  constants 


3670  PRINT 

3480  PRINT 

CHI" 

3490  PRINT 

eft)" 

3700  PRINT 


3710  ! 

3720  FOR  1=1  TO  I 

3730  PRINT  USING  -20X,3D,X,5(5D.D) ,7D,4X,3(5D.D) ,7D" 

D.Ttt  2(11, nil  »i  2(11, Ts  2(I),CFNe  2(11 

3740  NfjT  I 

3750  ! 

3740  RETURN 

3770  ! 

3780  ! 

3790  CONSTANTS:  ! 

3800  ! 

3810  Tratio=2.9 

3820  A1=.997B 

3830  fa--. 02943 

3S40  O*. 0001472 

3850  81=-. 1119 

3840  82=1.5315 

3870  B3=-.4702 

3880  Cl=.1334 

3890  C2=.7439 

3900  C3=.I209 

3910  CYCLE=.25  ! 

3920  C4=(Cl+C2tCYCLE+C3fCYCLEA2) /CYCLE 

3930  Tel  0ES=(95-32)»5/9  ! 

3940  Tel"DES= (46-32) t5/9  ! 

3950  !  " 

3960  !  UNIT  CONSTANTS 

3970  CAP  DES=21. 2*3.514  T 

3980  FLP~DES=20.3  ! 

3990  !  ' 

4000  !  COIL  CONSTANTS 

4010  «tr  cntrlJ="supply'~! 

4020  Tcnfrl  tetp=45  ! 

4030  GP".e=42\9  ! 

4040  ! 

4050  !  COOL  TNR  CONSTANTS 

4060  C=2. 15171  !      "   " 

4070  XN=.  99777  ! 

4080  CFH  t»r=4630  ! 

4090  6Pnc=55.4  ! 

4100  Tc  iin=75  ! 

4110  RETURN 

4120  ! 

4130  ! 

4140  ! 

4150  COIL  DATA: 

4140  TD8=TDB(xl 

4170  T«8=T»B(x) 

4180  IF  DATA  coil=2  THEN  SOTO  4250 

4190  Tdb=Tdb"l(x)  ! 

4200  ab  ai=no"  ai  l!x) 

4210  UNrTS»=UN~IT5  lt(x! 

4220  Ts=Ts  l(x)     " 

4230  CF«e=CF«e  Kxl 

4240  RETURN      " 

4250  Tdb=Tdb  2(x)  ! 

4260  «b  ai=nff  Hi  2(x) 

4270  UNrTS*=USlTS  2l(x) 

42S0  Ts=Ts  2(x)  " 

4290  CFNe=CF«e  2(x) 

4300  RETURN   " 

4310  ! 

4320  ! 

4330  COIL  CONDITIONS: 

4340  !   " 

4350  !  EI    hours    a«bient 


.aibient. 
TDB   TUB 

coil  t  1. 

Tdb   Tan   Ts 

CFM 

...coil  t  2 

CASE 

Tdb 

Grains  Ts 

> 

(cfa) 

deg.F 

fill    deg.F  1 











I,TOB(!) ,T»B(I) ,Tdb_l (I) ,»b_ai_l (I) ,Ts_l!I) ,CF*e_l ( 


PLR  at  ahich  systet  cycles 

syste«  cycling  constant 

tesp  at  »hicb  CAP  DES.RP  DES  selected 

te«p  at  which  CAP"DES,FLP"DES  selected 


chiller  capacity  in  kas  (  troi  iach.  data  ) 
chiller  potier    in  kas  (  froi  iach.  data  ) 


chilled  aater  teip  control— return  or  suaoly 
teip  at  which  above  Hater  is  held  at  (deg.F) 
chilled  aater  flow  rate     (  gpa  ) 


Hebb-Villacres  cooling  toaer  constant 
Bebb-Villacres  cooling  to»er  constant 
air  Ho*  rate  of  cooling  toner  (  cfi  ) 
aater  Hon  rate  of  toner, condenser  1  gpa  ) 
■iniiui  entering  condensing  Hater  teip  (deg.F) 


first  coil 


second  coil 


inlet 


.coil. 


outlet 
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dry-bulb- wet-bulb  dry-bulb, net-bulb  OR  grains, cf«  dry-bulb 
DATA   504    95   ,81     95   ,  75, KB,        ,2500   55 
DATA   504    95   ,81     95   ,         120,6  ,2500   55 


hrs  .33b.  ...tl  coil »2  coil. 


10o4,si,75.ai,75,llB.55,4oi, 82. 04,^55, 4497 
1524, 77, 73, 77, 73, KB, 55, 400, 84. 1,18. 6, 8,55,4026 


4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430  DATA  7, 97, 79.97, 79, KB, 55, 400, 78.4, 23, 6, 55, 6710 

4440  DATA  142, 92, 7a, 92, 78, IB, 55, 400, 79. 7, 22. 1,6.55, 6039 

4450  DATA  574, 87, 77, 87, 77, KB, 55, 400, 80. 9, 21, S, 55. 5343 

4460  DATA 

4470  DATA 

4480  DATA  1325,72,49;72;69,KB,55,400.35.7;i7.5;5,55,3355 

4490  DATA  923, 67, 43, 67, 43, KB, 55, 400, 87.7, 16.7,6, 55, 2684 

4500  DATA  779, 62, 58, 62, 58, KB, 55, 400, 92.2, 14. 2, 6, 55, 2013 

4510  DATA  701, 57, 53, 57, 53,NB,55, 400, 98.9, 8. 9, 3,55,1342 

4520  DATA  624, 52, 48,0, J, KB, 0,0, 99. 5, 6.1,6,41,1125 

4530  DATA  466, 47, 43, 0,0, KB, 0,0, 95. 6, 5.8, 8,70, 1125 

4540  DATA  342,42, 39, 0,0,KB,0, 0,93.4, 5.4, 6,79, 1125 

4550  ! 

4560  ! 

4570  1=12  !  nuiber  of  perfonance  calculations  above 

4580  ! 

4590  FDR  1=1  TO  J 

4400  READ  hrs(I) ,TDB(I) ,TKB(I) 

4410  READ  Tdb  1(1), nb  Hi  1(1), UNITS  U(I),Ts  lUKCFKe  1(1) 

4420  READ  Tdb"2(!) ,«b"»i  2(1) ,UNITS~2*(!) ,Ts "2(11 ,CFHe'2U) 

4430  NEXT  I   "     "  " 

4440  ! 

4450  ! 

4440  RETURN 
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APPENDIX   F 

ANNUAL  PERFORMANCE  DATA  FOR  THE  CHILLER  UNIT- 
TWO  COIL-COOLING  TOWER  SYSTEM 
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CHILLER CHILLED  HATER  COILS CO0LIN6  'ONER 

chilled  nater  supply  -sap.  held  at  45.0  deq.F 


SYSTEil    STEAOY    STATE    OPERATION 


:sse 

TIB      TUB           Ten 

P 

PLS 

Ten      Tci      T 

:pa: 

_in 

1       CF«  S»r 

CFKt 

Tel 

B 

9c 

1 

i 

97.0 

79.0         43.0 

94.5 

tons 
17.24 

kits 
17.13 

I 
82.4 

tons 
22.11 

94.5      84.9 

(9tw> 

.r.eri 

i 

.FHtar 

4.= 

100. 0 

4430 

4630 

194,0 

2 

92.0 

78.0         45.0 

92.9 

16.25 

15.92 

77.2 

20. 7B 

92.9      83.3 

3.8 

100.0 

4430 

4430 

100.0 

3 

87.0 

77.0         45.0 

91.2 

15.18 

14.67 

71.4 

19.35 

91.2      82.8      ' 

2.5 

100.0 

4430 

4430 

100.0 

4 

82.0 

75.0         45.0 

88.5 

13.90 

13.17 

44.4 

17.65 

83.5      80.3      80.3 

100.0 

4430 

4430 

100.0 

5 

77.0 

73.0         45.0 

85.7 

12.62 

11.75 

57.4 

15.96 

85.7      78.8-    78.3 

100.0 

4430 

4430 

100.0 

6 

72.0 

69.0         45.0 

80.9 

10.81 

9.SB 

43.0 

13.62 

30.9      74.9 

3.3 

98.2 

4430 

4321 

97.7 

7 

67.0 

63.0         45.0 

79.8 

8.74 

3.31 

38.4 

11.10 

79.3      71.9      75.0 

40.9 

4430 

2595 

34.1 

8 

62.0 

53.0         45.0 

78.9 

7.05 

7.10 

31.0 

9.07 

73.9      69.4 

3.) 

42.2 

4430 

1752 

37.3 

9 

57.0 

53.0         45.0 

73.1 

5.39 

5.31 

23.4 

7.04 

7S.1      67.5      75.0 

23.9 

4430 

1181 

35.3 

10 

52.0 

48.0         45.0 

77.2 

3.90 

4.18 

17.0 

5.09 

77.2      65.5 

3.3 

18.9 

4630 

767 

16.  o 

11 

47.0 

43.0         45.0 

76.5 

2.60 

2.77 

11.3 

3.33 

74.5      63.3 

3.3 

12.1 

1430 

483 

10.3 

12 

42.0 

39.0         45.0 

73.9 

1.46 

1.55 

4.3 

1.90 

75.9      62.4      73.0 

6.4 

4630 

259 

5.3 

.ill 
CASE         TBB 

COIL  CONDITIONS 

lent. 
TUB 

1  t  1.. 

"H 

Tdb      Srains 

Ts 

Tdb 

Tup 

Ts 

a 

CFB 

4 

1        97.0 

79.0 

.Idea. 
97.0 

Fl 

left! 

400 

deq.F    gr/lb 

:ea 

F    (cfil 

79.0 

55.0 

78.6      23.0 

35.0      4710 

2       92.0 

78.0 

92.0 

73.0 

55.0 

400 

79.7      22.1 

55 

0      4039 

3       87.0 

77.0 

87.0 

77.0 

55.0 

400 

80.9      21.0 

55.0      5348 

4       82.0 

75.0 

82.0 

75.0 

55.0 

400 

82.3      20.0 

55.0      4497 

5       77.0 

73.0 

77.0 

73.0 

55,0 

400 

34.1      13.6 

55.0      4024 

6       72.0 

69.0 

72.0 

69.0 

55.0 

400 

35.7      17.5 

55 

0      3355 

7       67.0 

65.0 

67.0 

63.0 

55.0 

400 

87.7      14.7 

53.0      2434 

8       62.0 

58.0 

62.0 

58.0 

55.0 

400 

92.2      14.2 

55.0      2013 

9       57.0 

53.0 

57.0 

53.0 

55.0 

400 

9B.9     '3.9 

55.0      1342 

10       52.0 

48.0 

0.0 

0.0 

0.0 

0 

99.5       6.1 

61 

0      1125 

11        47.0 

43.0 

0.0 

0.0 

0.0 

0 

95.6       5.3 

70 

0      1123 

12       42.0 

39.0 

0.0 

0.0 

1.0 

0 

93.4       5.4 

79.0      1123 

TUB  TWS 

STEADY  STATE 

ANNUAL 

CASE 

COILil    C0ILI2 

sensibl 

Olatent    Ototal 

ll/9t 

PONS          hours 

load 

vr 

p 

Mr/yr 

t 

.deq.F. 

▼ 

Ikusl           tin 

(MBTU) 

I 

ikn-nrs 



— 

— 

i        -.__ 



___ 

. 



„ 

1 

37 

79 

34.3        172.1 

190.5 

16 

.4 

204.9 

7.9 

17.1 

1413 

.13 

120 

16 

2 

n 

73 

32.9        162.0 

173.3 

16 

.7 

95.0 

3 

15.9             142 

27688 

.38 

2241 

H 

3 

37 

77 

31.1        151.0 

165.0 

17 

.1 

82.1 

9 

14.7            574 

104532    1 

.36 

3413 

it. 

:4 

4 

32 

75 

27.6       139.2 

151.1 

13 

.7 

66.8 

9 

13.2           1008 

68161     1 

.47 

13274 

17. 

3: 

S 

77 

73 

24.3       127.2 

136.3 

14 

.6 

51.4 

9 

11.8           1326 

231043    2 

.00 

17934 

4 

n 

69 

13.0        111.9 

119.2 

10 

•  5 

:°.i 

3 

9.9           1325 

"1919    1 

.34 

13097 

1;# 

04 

7 

a 

43 

9.6         95.2 

100.5 

4 

04.8 

4 

3.3           923 

94730    10.05 

7471 

9, 

99 

9 

a 

53 

3.4         81.2 

84.3 

,3 

34.6 

7.1             779 

43900 

.33 

20 

9 

u 

53 

.9         63.S 

64.7 

0.0 

64.7 

0 

0 

5.3             701 

15322 

.71 

1073 

j. 

30 

10 

52 

43 

0.0         46.9 

44.9 

0.0 

46.9 

0 

0 

4.2            624 

29234 

.04 

2611 

3. 

40 

11 

17 

43 

0.0         31.2 

31.2 

0.0 

31.2 

0 

0 

2.3            466 

14317 

.31 

1291 

\t 

o3 

12 

12 

a 

0.0         17.3 

17.3 

0.0 

I7.J 

0 

0 

1.6            342 

3992 

.62 

331 

o9 

Annual    coolinq 

coil  t 

1    load 

,     1 

23.63  imtll 

Annual    coolinq 

coil  t 

2    load 

=      333.88  MBtu  . 

Annual  sensible  c 

polinq 

load 

"      393.37  mBtu 

Annual  latent      coolinq 

load 

= 

63.64  KltBtu 

Annual  coo 

inp  load 

=      962.31  BNBtu 

Annual  Hours  at  coerati 

n 

: 

3417  Hours 

Annual  Chiller  Ppner  consuiption  = 

76801  Kn-hrs 

(single  bin  data  run) 
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.an! 
TDB 

COIL  CONDITIONS 

nent. 
TUB 

.1  *  1.. 

Ts 

1  2 

* 

CSSE 

Tdh 

Tub 

CFH 

Tdb 

Srains 

Ts 

CFH 

1 

97.0 

32.0 

97.0 

,F) 

82.0 

55.0 

;t'i) 

deg.F 

79.5 

gr/lb 

deg.F 

(cfi) 

1 

400 

22.1 

55.0 

6710 

2 

97.0 

78.0 

97.0 

78.0 

55.0 

400 

78.4 

23.2 

55.0 

6710 

3 

92.0 

86.0 

92.0 

84.0 

55.0 

400 

80.6 

21.1 

55.0 

6039 

4 

92.0 

84.0 

92.0 

84.0 

55.0 

400 

80.6 

21.1 

55.0 

4039 

5 

92.0 

32.0 

92.0 

82.0 

55.0 

400 

30.6 

21.1 

55.0 

6039 

6 

92.0 

80.0 

92.0 

80.0 

55.0 

400 

80.3 

21.5 

55.0 

6039 

7 

92.0 

78.0 

92.0 

78.0 

55.0 

400 

79.7 

22.1 

55.0 

6039 

B 

92.0 

76.0 

92.0 

76.0 

55.0 

400 

79.1 

22.7 

55.0 

6039 

9 

92.0 

74.0 

92.0 

74.0 

55.0 

400 

78.6 

23.2 

55.0 

6039 

10 

92.0 

72.0 

92.0 

72.0 

55.0 

400 

78.1 

23.8 

55.0 

6039 

11 

B7.0 

34.0 

87.0 

84.0 

55.0 

400 

81.9 

19.9 

55.0 

5368 

12 

87.0 

82.0 

87.0 

82.0 

55.0 

400 

81.9 

19.9 

55.0 

5368 

13 

37.0 

30.0 

87.0 

80.0 

55.0 

400 

81.9 

20.0 

55.0 

5368 

14 

87.0 

78.0 

87.0 

78.0 

55.0 

400 

81.3 

20.7 

55.0. 

5368 

15 

87.0 

76.0 

87.0 

76.0 

55.0 

400 

80.6 

21.4 

55.0 

5368 

16 

87.0 

74.0 

87.0 

74.0 

55.0 

400 

80.0 

22.0 

55.0 

5368 

17 

87.0 

72.0 

37.0 

72.0 

55.0 

400 

79.5 

22.6 

55.0 

5368 

18 

87.0 

70.0 

87.0 

70.0 

55.0 

400 

79.0 

23.2 

55.0 

5368 

19 

87.0 

68.0 

87.0 

68.0 

55.0 

400 

78.6 

23.7 

55.0 

5368 

20 

87.0 

66.0 

87.0 

66.0 

55.0 

400 

78.0 

24.3 

55.0 

5368 

21 

82.0 

82.0 

82.0 

32.0 

55.0 

400 

33.7 

18.5 

55.0 

4697 

22 

82.0 

80.0 

82.0 

80.0 

55.0 

400 

83.7 

18.5 

55.0 

4697 

23 

32.0 

78.0 

32.0 

78.0 

55.0 

400 

83.4 

18.9 

55.0 

4697 

24 

32.0 

76.0 

82.0 

76.0 

55.0 

400 

82.7 

19.6 

55.0 

4697 

25 

82.0 

74.0 

32.0 

74.0 

55.0 

400 

31.9 

20.4 

55.0 

4697 

26 

82.0 

72.0 

82.0 

72.0 

55.0 

400 

81.3 

21.1 

55.0 

4697 

27 

82.0 

70.0 

32.0 

70.0 

55.0 

400 

80.7 

51.7 

55.0 

4697 

26 

82.0 

68.0 

82.0 

63.0 

55.0 

400 

30.1 

22.3 

55.0 

4697 

29 

82.0 

66.0 

32.0 

66.0 

55.0 

400 

79.7 

23.0 

55.0 

4697 

30 

82.0 

64.0 

82.0 

64.0 

55.0 

400 

79.0 

23.5 

55.0 

4697 

31 

82.0 

62.0 

82.0 

62.0 

55.0 

400 

77.7 

24.4 

55.0 

4697 

32 

78.0 

78.0 

78.0 

78.0 

55.0 

400 

85.4 

17.0 

55.0 

4160 

33 

77.0 

76.0 

77.0 

76.0 

55.0 

400 

35.4 

17.3 

55.0 

4026 

34 

77.0 

74.0 

77.0 

74.0 

55.0 

400 

84.5 

18.2 

55.0 

4026 

35 

77.0 

72.0 

77.0 

72.0 

55.0 

400 

33.7 

19.0 

55.0 

4026 

36 

77.0 

70.0 

77.0 

70.0 

55.0 

400 

83.0 

19.7 

55.0 

4024 

37 

77.0 

68.0 

77.0 

66.0 

55.0 

400 

62.3 

20.5 

55.0 

4026 

38 

77.0 

66.0 

77.0 

66.0 

55.0 

400 

81.7 

21.2 

55.0 

4026 

39 

77.0 

64.0 

77.0 

64.0 

55.0 

400 

31.3 

21.9 

55.0 

4026 

40 

77.0 

62.0 

77.0 

62.0 

55.0 

400 

30.3 

22.6 

55.0 

4026 

41 

77.0 

60.0 

77.0 

60.0 

55.0 

400 

78.9 

23.4 

55.0 

4026 

42 

77.0 

58.0 

77.0 

58.0 

55.0 

400 

7B.2 

23.4 

55.0 

4026 

43 

74.0 

74.0 

74.0 

74.0 

55.0 

400 

87.7 

15.1 

55.0 

3623 

44 

72.0 

72.0 

72.0 

72.0 

55.0 

400 

87.1 

16.1 

55.0 

3355 

45 

72.0 

70.0 

72.0 

70.0 

55.0 

400 

36.2 

17.1 

55.0 

3355 

46 

72.0 

68.0 

72.0 

68.0 

55.0 

400 

35.3 

17.9 

55.0 

3355 

47 

72.0 

66.0 

72.0 

66.0 

55.0 

400 

84.6 

18.8 

55.0 

3355 

48 

72.0 

64.0 

72.0 

64.0 

55.0 

400 

34.0 

19.6 

55.0 

3355 

49 

72.0- 

62.0 

72.0 

62.0 

55.0 

400 

33.5 

20.4 

55.0 

3355 

50 

72.0 

60.0 

72.0 

60.0 

55.0 

400 

32.3 

21.1 

55.0 

3355 

51 

72.0 

58.0 

72.0 

58.0 

55.0 

400 

80.7 

22.1 

55.0 

3355 

52 

72.0 

56.0 

72.0 

56.0 

55.0 

400 

79.9 

22.1 

55.0 

3355 

53 

72.0 

54.0 

72.0 

54.0 

55.0 

400 

79.2 

22.0 

55.0 

3355 

54 

72.0 

52.0 

72.0 

52.0 

55.0 

400 

76.7 

21.9 

55.0 

3355 

55 

70.0 

70.0 

70.0 

70.0 

55.0 

400 

90.8 

12.5 

55.0 

3086 

56 

68.0 

68.0 

68.0 

63. 0 

55.0 

400 

88.9 

15.0 

55.0 

2S!3 

(double  bin  data  run) 
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.aal 
TDB 

COIL  CONDITIONS 

lient. 
TUB 

1  1 1.. 

Ts 

cfh 

1  2 

:ase 

Tdfa 

Tub 

Tdb 

Brains 

Ts 

CF« 

i 

Fl 

icfi) 

deg.F 

gr/lb 

deg.F 

Icfi) 

57 

67.0 

66.0 

67.0 

66.0 

55.0 

400 

89.0 

15.1 

55.0 

2684 

58 

67.0 

64.0 

67.0 

64.0 

55.0 

400 

38.1 

16.2 

55.0 

2684 

59 

67.0 

62.0 

67.0 

62.0 

55.0 

400 

87.4 

17.1 

55.0 

2634 

60 

67.0 

60.0 

67.0 

60.0 

55.0 

400 

86.7 

18.1 

55.0 

2684 

61 

67.0 

58.0 

67.0 

58.0 

55.0 

400 

85.2 

19.0 

55.0 

2684 

62 

67.0 

56.0 

67.0 

56.0 

55.0 

400 

83.4 

20.1 

55.0 

2684 

63 

67.0 

54.0 

67.0 

54.0 

55.0 

400 

82.3 

20.1 

55.0 

2684 

64 

67.0 

52.0 

67.0 

52.0 

55.0 

400 

81.5 

20.0 

55.0 

2684 

65 

67.0 

50.0 

67.0 

50.0 

55.0 

400 

80.9 

19.9 

55.0 

2684 

66 

67.0 

48.0 

67.0 

48.0 

55.0 

400 

80.4 

19.3 

55.0 

2684 

67 

64.0 

64.0 

64.0 

64.0 

55.0 

400 

93.0 

11.7 

55.0 

2281 

a 

62.0 

62.0 

62.0 

62.0 

55.0 

400 

94.1 

11.7 

55.0 

2013 

69 

62.0 

60.0 

62.0 

60.0 

55.0 

400 

93.2 

13.0 

55.0 

2013 

70 

62.0 

58.0 

62.0 

58.0 

55.0 

400 

92.3 

14.2 

55.0 

2013 

71 

62.0 

56.0 

62.0 

56.0 

55.0 

400 

90.3 

15.4 

55.0 

2013 

?2 

62.0 

54.0 

62.0 

54.0 

55.0 

400 

87.8 

16.8 

55.0 

2013 

73 

62.0 

52.0 

62.0 

52.0 

55.0 

400 

86.5 

16.8 

55.0 

2013 

74 

62.0 

50.0 

62.0 

50.0 

55.0 

400 

85.4 

16.7 

55.0 

2013 

73 

62.0 

48.0 

62.0 

48.0 

55.0 

400 

84.7 

16.5 

55.0 

2013 

74 

62.0 

46.0 

62.0 

46.0 

55.0 

400 

84.0 

16.4 

55.0 

2013 

77 

60.0 

60.0 

60.0 

60.0 

55.0 

400 

97.6 

9.1 

55.0 

1744 

73 

58.0 

58.0 

58,0 

58.0 

55.0 

400 

101.3 

7.2 

55.0 

1476 

79 

57.0 

56.0 

57.0 

56.0 

55.0 

400 

102.9 

6.3 

55.0 

1365 

30 

57.0 

54.0 

57.0 

54.0 

55.0 

400 

100.5 

3.0 

55.0 

1342 

31 

57.0 

52.0 

57.0 

52.0 

55.0 

400 

96.9 

10.3 

55.0 

1342 

82 

57.0 

50.0 

57.0 

50.0 

55.0 

400 

95.0 

10.2 

55.0 

1342 

83 

57.0 

48.0 

57.0 

48.0 

55.0 

400 

93.3 

10.1 

55.0 

1342 

34 

57.0 

46.0 

57.0 

46.0 

55.0 

400 

92.3 

9.8 

55.0 

1342 

85 

57.0 

44.0 

57.0 

44.0 

55.0 

400 

91.2 

9.6 

55.0 

1342 

86 

57.0 

42.0 

57.0 

42.0 

55.0 

400 

90.2 

9.4 

55.0 

1342 

37 

54.0 

54.0 

0.0 

0.0 

0.0 

0 

102.6 

6.6 

59.0 

1320 

33 

52.0 

52.0 

0.0 

0.0 

0.0 

0 

102.3 

6.3 

62.0 

1245 

39 

52.0 

50.0 

0.0 

0.0 

0.0 

.  0 

101.6 

6.2 

61.0 

1125 

90 

52.0 

48.0 

0.0 

0.0 

0.0 

0 

99.4 

6.2 

61.0 

1125 

91 

52.0 

46.0 

0.0 

0.0 

0.0 

0 

97.3 

6.2 

61.0 

1125 

92 

52.0 

44.0 

0.0 

0.0 

0.0 

0 

96.1 

5.8 

61.0 

1125 

93 

52.0 

42.0 

0.0 

0.0 

0.0 

0 

94.3 

5.6 

61.0 

1125 

94 

52.0 

40.0 

0.0 

0.0 

0.0 

0 

93.6 

5.3 

61.0 

1125 

95 

50.0 

50. 0 

0.0 

0.0 

0.0 

0 

100.8 

3.3 

65.0 

1200 

96 

48.0 

48.0 

0.0 

0.0 

0.0 

0 

100.5 

5.9 

68.0 

1125 

97 

47.0 

46.0 

0.0 

0.0 

0.0 

0 

98.8 

5.7 

70.0 

1125 

98 

47.0 

44.0 

0.0 

0.0 

0.0 

0 

96.6 

5.3 

70.0 

1125 

99 

47.0 

42.0 

0.0 

0.0 

0.0 

0 

94.9 

5.7 

70.0 

1125 

too 

47.0 

40.0 

0.0 

0.0 

0.0 

0 

JJ.7 

5.5 

70.0 

1125 

101 

47.0 

38.0 

0.0 

0.0 

0.0 

0 

=2.5 

5.2 

70.0 

1125 

102 

47.0 

36.0 

0.0 

0.0 

0.0 

0 

91.4 

5.1 

70.0 

1125 

103 

44.0 

44.0 

0.0 

0.0 

0.0 

0 

97.6 

5.7 

75.0 

1125 

104 

42.0 

42.0 

0.0 

0.0 

0.0 

0 

96.0 

5.4 

79.0 

1125 

105 

42.0 

40.0 

0.0 

0.0 

0.0 

0 

94.1 

5.5 

79.0 

1125 

106 

42.0 

38.0 

0.0 

0.0 

0.0 

0 

92.9 

5.2 

79.0 

1125 

107 

42.0 

36.0 

0.0 

0.0 

0.0 

0 

91.7 

5.1 

79.0 

1125 

108 

42.0 

34.0 

0.0 

0.0 

0.0 

0 

36.2 

4.9 

79.0 

1125 

109 

40.0 

40.0 

0.0 

0.0 

0.0 

0 

94.3 

5.4 

83.0 

1125 

110 

38.0 

38.0 

0.0 

0.0 

0.0 

0 

93.3 

5.2 

86.0 

1125 

111 

37.0 

36.0 

0.0 

0.0 

0.0 

0 

91.8 

5.1 

88.0 

1125 

(double  bin  data  run) 
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5YSTEB 

CHIUEH HATES-COIL — CG0LIN6  TMEB 

chilled  niter  supply  teep.  held  at  '5.0  deg.F 


SYSTEJt    STEflDT    STATE    BPEKflTION 


.aiftieflt.  systet coelinq  toner 

CASE  TUB  T»8  Ten  Tel           0e         P  PUS  Oc  Tco  Tci  Tc  in  I  Cffl  twr  Cfltt  I 

t  dtd.f tons  kin  I  tons  ....deq.F .".  (Star)  ..".c(«...       Iffltirl 

1  97.0  82.0  *S.O  97.4  18.27  18.77  85.1  23.41  97.4  37.4  87.4  100.0  (430  (430  100.0 

2  97.0  78.0  43.0  93.3  14.94  14.48  80.8  21.71  93.3  84.0  84.0  100.0  4430'  4430  100.0 

3  92.0  84.0  45.0  100.4  18.10  19.23  89.8  23.57  100.4  90.5  90.5  100.0  4430  4430  100.0 

4  92.0  84.0  43.0  98.7  17.73  18.48  87.1  22.99  98.7  88.3  8B.8  100.0  4430  4430  100.0 

5  92.0  82.0  45.0  94.9  17.38  17.73  84.4  22.43  94.9  87.2  87.2  100.0  4430  4430  100.0 
4  92.0  80.0  45.0  95.0  14.89  14.9!  81.2  21.70  95.0  S5.4  85.4  100.0  4430  4430  100. 0 

7  92.0  78.0  43.0  92.9  14.23  15.92  77.2  20.78  92.9  33.S  83.8  1OO.0  4430  4430  100.0 

8  92.0  74.0  43.0  90.3  13.42  14.99  73.3  19.89  90.3  82.1  82.1  100.0  4430  4430  100.0 

9  92.0  74.0  43.0  88.3  15.04  14.18  49.9  19.10  88.3  30.4  30.4  100.0  4430  4430  100.0 

10  92.0  72.0  45.0  84.4  14.52  13.40  44.4  13.33  34.4  73.7  78.7  100.0  4430  4430  100.0 

11  87.0  84.0  45.0  98.0  14.30  17.42  82.2  21.74  -  98.0  88.7  83.7  100.0  4430  4430  100.0 

12  87.0  82.0  43.0  94.1  14.44  14.71  79.4  21.21  94.1  37.0  87.0  100.0  4430  4430  100.0 

13  87.0  30.0  45.0  94.4  14.13  14.08  77.2  20.70'  94.4  85.4  35.4  100.0  4430  4430  100.0 

14  87.0  78.0  45.0  92.3  15.52  15.17  73.5  19.83  92.3  33.4  33.4  100.0  4430  4430  100.0 

15  87.0  74.0  45.0  90.1  14.38  14.24  49.4  18.93  90.1  81.9  81.9  100.0  4430  4430  100.0 
14  37.0  74.0  43.0  88.1  14.31  13.44  44.1  18.13  38.1  80.2  30.2  100.0  4430  4430  100.0 
17  87.0  72.0  45.0  84.0  13.79  12.72  43.0  17.41  34.0  78.4  78.4  100.0  4430  4430  100.0 
ia  37.0  70.0  43.0  84.0  13.29  12.04  40.1  14.72  84.0  74.7  74.7  100. 0  4430  4430  10O.0 

19  87.0  48.0  45.0  82.1  12.85  11.50  37.5  14.12  32.1  73.0  75.0  99.3  4430  4594  99.2 

20  87.0  44.0  45.0  81.7  12.35  11.09  33.1  15.51  81.7  74.0  73.0  34.7  4430  3878  33.3 

21  82.0  82.0  43.0  93.4  13.58  15.77  73.1  20.07  95.4  84.3  34.8  100.0  4430  4430  100.0 

22  32.0  80.0  43.0  93.4  15.23  13.15  72.7  19.54  93.4  83.1  35.1  100.0  4430  4430  100.0 

23  82.0  73.0  43.0  91.8  14.31  14.44  49.9  18.92  91. a  33.5  33.5  100.0  4430  4430  100.0 

24  82.0  74.0  43.0  39.4  14.22  13.39  44.3  18.08  89.4  31.7  31.7  100.0  4430  4430  100.0 

23  32.0  74.0  43.0  B7.4  13.59  12.73  42.4  17.22  87.4  79.9  79.9  1OO.0  4430  4430  100.0 

24  82.0  72.0  43.0  85.3  13.07  12.04  59.5  14.50  83.3  78.1  78.1  100.0  4430  4430  1OO.0 

27  82.0  70.0  45.0  83.3  12.53  11.41  54.5  15.80  83.3  74.4  74.4  100.0  4430  4430  100.0 

28  32.0  48.0  45.0  81.4  12.03  10.85  53.7  13.14  81.4  74.7  75.0  94.3  4430  4413  95.3 

29  82.0  44.0  45.0  31.3  11.44  10.53  31.9  14.44  81.3  73.8  75.0  34.1  4430  3740  30.3 

30  32.0  44.0  43.0  31.1  11.20  10.18  49.3  14.09  81.1  72.9  73.0  74.5  4430  3244  70.1 

31  82.0  42.0  45.0  80.3  10.44  9.73  47.3  13.42  80.8  72.0  73.0  43.9  4430  2827  41.1 

32  78.0  78.0  43.0  91.1  14.18  13.80  44.7  18.10  91.1  83.2  83.2  100.0  4430  4430  100.0 

33  77.0  74.0  43.0  89.0  13.31  12.91  42.3  17.18  89.0  81.3  31.3  100.0  4430  4430  100.0 

34  77.0  74.0  45.0  84.3  12.90  12.11  59.2  14.34  84.8  79.4  79.4  100.0  4430  4430  10O.0 

33  77.0  72.0  43.0  34.7  12.34  11.41  34.0  15.58  34.7  77.9  77.9  100.0  4430  4430  1OO.0 

34  77.0  70.0  43.0  82.4  11.32  10.79  53.0  14.39  82.4  74.1  74.1  100.0  4430  4430  100.0 

37  77.0  48.0  45.0  81.2  11.32  10.27  50.4  14.24  81.2  74.4  73.0  93.4  4430  4230  91.3 

38  77.0  44.0  45.0  80.9  10.37  9.92  48.3  13.49  30.9  73.4  75.0  81.0  4430  1581  77.3 

39  77.0  44.0  43.0  30.3  10.30  9.43  44.4  13.24  30.8  72.7  73.0  71.3  4430  3111  47.2 

40  77.0  42.0  45.0  80.5  10.03  9.28  44.5  12.44  80.5  71.8  75.0  43.5  4430  2713  58.4 

41  77.0  4O.0  43.0  30.2   9.52  3.89  42.1  12.04  80.2  70.9  73.0  54.4  4430  2232  51.4 

42  77.0  58.0  45.0  80.1   9.24  8.70  41.0  11.73  30.1  70.2  73.0  51.3  4430  2150  44.4 

43  74.0  74.0  45.0  84.7  12.34  12.08  59.0  14.30  84.7  79.4  79.4  100.0  4430  4430  100.0 

44  72.0  72.0  45.0  84.0  11.42  10.73  52.5  14.49  84.0  77.4  77.4  100.0  4430  4430  100.0 

43  72.0  70.0  43.0  81.9  11.09  10.14  49.5  13.98  31.9  73.7  73.7  100.0  4430  4430  100.0 

44  72.0  48.0  45.0  30.8  10.57   9.49  44.9  13.32  80.3  74.4  73.0  '0.4  4430  4034  38.2 

47  72.0  44.0  45.0  30.5  10.12  9.35  44.9  12.78  80.5  73.4  73.0  78.0  4430  3422  73.9 

48  72.0  44.0  43.0  30.4   9.71   9.04  43.0  12.28  80.4  72.5  75.0  48.5  4430  2933  43.3 

49  72.0  42.0  43.0  30.2   9.33  8.74  41.4  11.84  80.2  71.7  73.0  40.7  4430  2333  33.3 

50  72.0  4O.0  45.0  79.9   3.91   8.44  39.4  11.31  79.9  70.8  73.0  53.9  4430  2272  49.1 

51  72.0  5B.0  45.0  79.7   8.43  8.08  37.2  10.72  79.7  70.0  73.0  48.0  4430  2004  43.3 

32  72.0  54.0  45.0  79.5   8.18   7.91  34.1  10.43  79.3  49.2  75.0  43.3  4430  317  17.4 

33  72.0  54.0  45.0  79.4  7.97  7.75  33.1  10.17  79.4  48.5  73.0  40.3  4430  1448  34.0 
54  72.0  52.0  45.0  79.3   7.81   7.44  34.4   9.99  79.3  47.8  73.0  37.5  4430  1347  33.4 

53  70.0  70.0  45.0  82.3  11.40  10.59  52.0  14.41  32.3  74.0  74.0  1OO.0  4430  4430  100.0 

54  48.0  48.0  45.0  80.5  10.01   9.24  44.4  12.44  30.5  74.3  75.0  68.1  4430  3939  35.3 


(double  bin  data  run) 
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SYSTEH 

CHILLER KBBKOIl— CHUM  TDNER 

cnilled  Niter  supply  teig.  held  it  45.0  :eq.7 


SYSTBI    STEADY    STATE    OPERATION 


CASE 

.aaoient. 

tub    ma 

Ten 

Tel 

Ot 

Oc 

Ten 

Tci 

Tc.in       I 

p 

PL8 

Cfll  tnr 

cm 

: 

1 
57 

67.0 

50.  0 

dtq.f.... 
43.0 

80.2 

tans 

9.38 

OS 

8.79 

I 
41.5 

tons 
11.88 

80.2      73.3 

IDturl 

.r.cfi 

... 

CFlttw! 

73.0 

74.7 

4430 

3239 

70.4 

sa 

67.0 

64.0 

45.0 

79.9 

8.94 

3.44 

39.5 

11.35 

79.9 

72.3 

75.0 

65.1 

4430 

2791 

40.3 

59 

67.0 

62.0 

45.0 

79.7 

8.56 

8.13 

37.3 

10.38 

79.7 

71.5 

75.0 

57.3 

4430 

2425 

52.4 

M 

67.0 

60.0 

45.0 

79.5 

8.18 

7.90 

34.1 

10.43 

79.5 

70.4 

73.0 

30.9 

4430 

2132 

44.1 

u 

67.0 

58.0 

45.0 

79.3 

7.77 

7.41 

34.2 

9.93 

79.3 

49.8 

75.0 

45.3 

4630 

1384 

40.7 

62 

67.0 

56.0 

45.0 

79.1 

7.33 

7.30 

32.3 

9.41 

79.1 

49.0 

75.0 

40.4 

4630 

1473 

34.1 

63 

67.0 

54.0 

45.0 

78.9 

7.07 

7.11 

St.] 

9.09 

7B.9 

48.3 

75.0 

34.3 

4430 

1513 

32.3 

64 

67.0 

52.0 

43.0 

78.3 

6.87 

4.97 

30.2 

3.35 

73.8 

47.6 

75.0 

34.0 

4630 

1398 

30.2 

65 

47.0 

50.0 

45.0 

78.8 

4.72 

4.87 

29.5 

3.58 

78.8 

66.9 

75.0 

31.7 

4630 

838 

18.1 

66 

67.0 

48.0 

45.0 

7B.7 

4.40 

4.78 

7.9.0 

8.53 

78.7 

44.3 

75.0 

29.3 

4430 

7B4 

14.9 

67 

64.0 

64.0 

45.0 

79.3 

8.74 

8.31 

73.4 

11.10 

79.8 

72.3 

75.0 

64.2 

4430 

2747 

59.3 

H 

62.0 

62.0 

45.0 

79.3 

7.80 

7.43 

34.4 

9.97 

79.3 

71.3 

75.0 

53.9 

4430 

226S 

49.0 

69 

62.0 

60.0 

45.0 

■  79.1 

7.43 

7.37 

32.7 

J.32 

79.1 

70.4 

75.0 

47.6 

4430 

1986 

42.9 

70 

62.0 

58.0 

45.0 

78.9 

7.07 

7.11 

31.1 

9.09 

78.9 

49.4 

75.0 

42.3 

4630 

1734 

37.9 

71 

62.0 

56.0 

45.0 

78.7 

6.48 

4.83 

29.3 

8.62 

73.7 

68.3 

75.0 

37.7 

4420 

1535 

33.6 

72 

62.0 

54.0 

45.0 

78.5 

6.22 

4.52 

27.3 

3.08 

78.5 

68.0 

75.0 

33.5 

4430 

337 

19.2 

73 

.62.0 

52.0 

45.0 

78.4 

5.99 

4.34 

24.2 

7.79 

78.4 

47.3 

75.0 

30.6 

4630 

807 

17.4 

74 

62.0 

50.0 

45.0 

78.3 

5.79 

4.22 

23.3 

7.55 

78.3 

64.7 

75.0 

28.2 

4630 

1132 

24.9 

75 

62.0 

48.0 

45.0 

78.2 

5.64 

4.11 

24.8 

7.40 

78.2 

64.0 

75.0 

24.3 

4430 

1074 

76 

62.0 

44.0 

45.0 

78.1 

5.53 

5.97 

34.2 

7.23 

78.1 

4S.4 

75.0 

24.6 

4630 

1004 

2l!7 

77 

60.0 

60.0 

45.0 

79.0 

7.19 

7.19 

J1.4 

9.23 

79.0 

70.4 

73.0 

44.5 

4630 

1937 

41.8 

78 

SB.O 

58.0 

45.0 

78.4 

6.45 

4.47 

28.3 

8.33 

78.6 

49.5 

73.0 

39.4 

4430 

1435 

33.3 

79 

57.0 

56.0 

45.0 

78.4 

5.99 

4.34 

24.2 

7.30 

73.4 

43.7 

73.0 

34.3 

4630 

1429 

30.9 

80 

57.0 

54.0 

45.0 

78.2 

5.58 

4.02 

24.4 

7.29 

78.2 

47.9 

75.0 

30.7 

4630 

311 

17.5 

31 

57.0 

52.0 

45.0 

77.9 

5.15 

5.:4 

22.5 

4.72 

77.9 

47.1 

75.0 

27.0 

4430 

1101 

23.3 

82 

57.0 

50.0 

45.0 

77.3 

4.92 

5.29 

21.5 

4.42 

77.8 

44.4 

75.0 

24.3 

4430 

997 

21.3 

83 

57.0 

48.0 

45.0 

77.7 

4.71 

5.04 

20.4 

6.15 

77.7 

45.3 

75.0' 

22.4 

4430 

911 

19.7 

84 

57.0 

46.0 

45.0 

77.4 

4.59 

4.93 

20.0 

5.99 

77.4 

65.2 

73.0 

20.9 

4430 

348 

18.3 

85 

57.0 

44.0 

45.0 

77.3 

4.46 

4.79 

19.4 

5.82 

77.5 

64.4 

75.0 

19.5 

4630 

739 

17.0 

36 

57.0 

42.0 

45.0 

77.5 

4.33 

4.63 

18.9 

5.66 

77.3 

64.0 

75.0 

18.2 

4630 

738 

13.9 

87 

54.0 

54.0 

45.0 

77.9 

5.19 

5.59 

22.7 

6.78 

77.9 

47.8 

73.0 

23.9 

4630 

1130 

25.5 

88 

52.0 

52.0 

45.0 

77.4 

4.52 

4. 84 

19.7 

3.91 

77.4 

46.9 

73.0 

24.1 

4430 

980 

21.2 

89 
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(8 

48 

0.0 

37.5 

39.5 

0.0 

29.5 

0.0 

Z.  5 

60 

2373 

.24 

212 

.27 

97 

n 

46 

0.0 

35.0 

35.0 

0.0 

35.0 

3.0 

101 

3340 

.36 

315 

.40 

98 

47 

14 

0.0 

32.4 

32.4 

0.0 

32.4 

0.0 

:.'° 
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4078 

.41 

343 

.46 

99 

47 

42 

0.0 

30.3 

30.3 

0.0 

30.3 

0.0 

2.7 

103 

3121 

.31 

277 

.35 

100 

47 

40 

0.0 

23.8 

28.3 

0.0 

23.3 

0.0 

2.6 

34 

980 

.10 

37 

.11 

101 

47 

38 

0.0 

27.5 

27.5 

0.0 

37.3 

0.0 

2.4 

24 

715 

.07 

43 

.08 

102 

47 

56 

0.0 

24.0 

26.0 

0.0 

26.0 

0,3 

2.3 

5 

130 

.01 

12 

.01 

103 

44 

44 

0.0 

27.3 

27.5 

0.0 

27.5 

0.0 

2.4 

24 

660 

.07 

=Q 

.07 

10* 

42 

42 

0.0 

20.7 

20.7 

0.0 

20.7 

0.0 

1.8 

75 

1551 

.16 

137 

105 

42 

40 

0.0 

13.4 

13.4 

0.0 

18.4 

0.0 

1.6 

54 

992 

.10 

33 

J! 

106 

42 

53 

0.0 

14.9 

14.9 

3.0 

16.9 

0.3 

1.5 

82 

1387 

.14 

123 

.  13 

107 

42 

56 

0.0 

15.5 

15.5 

0.0 

15.5 

0.0 

1.4 

54 

334 

.08 

74 

yp 

108 

42 

54 

0.0 

3.8 

3.3 

0.0 

3.3 

0.0 

.3 

11 

96 

.01 

9 

.01 

109 

40 

40 

0.0 

13.7 

13.7 

0.0 

13.7 

0.0 

1.2 

96 

.01 

9 

.01 

lit 

a 

a 

0.0 

3.9 

3.9 

0.0 

8.9 

3.0 

.3 

34 

302 

.03 

27 

.03 

111 

" 

54 

0.0 

4.4 

4.4 

0.0 

4.0 

0.0 

.4 

15 

208 

.02 

18 

Annual  pre-coDling 

coil  load 

• 

33.66  fflBtu 

Annual  after-coolinq 

coil  load 

58.72  HIlBtu 

Annual  sensible  cool 

nq  load 

14.76  fflBtu 

Annual  latent 

cool: 

nq  load 

* 

77.62  MIBtu 

Annual  cooling  load  =      :a?.;3  HIlBtu 

Annual  hours  a-f  operation  *         8537  Hours 

Annual  Chiller  Pcner  consuaption  =       79oS9  Ka-hrs 


(double  bin  data  run) 
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APPENDIX   G 


COMPUTER  PROGRAM  FOR  DETERMINING  THE 
WEBB-VILLACRES  COOLING  TOWER  CONSANTS 


10!  125 

20  !  'TOMER  FILL' 

30 

40  !  THIS  PRQSRAd  IS  BASED  ON  THE  MODEL  DEVELOPED  BY  R.L.MEBB  WO  A.VILLACRES 

50  !   AMD  IS  OSED  TO  CALCULATE  THE  FILL  CHARACTERISTICS  OF  A  C00LIN6  TOWER 

60 

70 

80  DIK  TOWER* [301 

90  ! 

100  DISP  ■  Please  input  the  nane  of  the  cooling  toner  sanuf acture. ■ 

110  DISP 

120  INPUT  TOWER* 

130  ! 

140  DISP  '  Please  input  the  «odel  nuiber  of  the  cooling  toner.' 

150  DISP 

160  input  moat 

170  ! 

180  DISP  "  Please  input  the  intonation  requested  for  design  conditions" 

190  DISP 

200  ! 

240  DISP  '  Please  input  the  design  cooling  toner  AIR  FLOW  RATE  in  CFH  ■ 

250  DISP 

260  INPUT  CFH  toner 

262  ! 

270  DISP  "  Please  input  the  design  aibient  MET-BUL3  teap. .RANEE, APPROACH. I  SATES  FLQK  RATE 

272  DISP  ■  (deg.F)     tdeg.FI  (deg.F)'      (gpi)" 

280  DISP 

290  INPUT  T«B,R,Al,SPm 

300  ! 

310  DISP  ■  Please  input  another  APPROACH  TEHPERATURE  and  MATER  FLOW  RATE  for' 

320  DISP  '   the  COOLINS  TOWER  to  achieve  that  RANGE, WET-BULB — (den . F  I  6PK)  * 

330  INPUT  A2,6PH2 

340  ! 

350  SOSUB  TOWERID  ! 

360  ! 

370  Rlg=6PHl*9.33/(CFN  toner*. 07121 

380  A=A1 

390  Rlq!ll=Rlg 

400  SOSUB  COUNT 

410  Ntur(l)=CT 

420  Rlg=6PH2*8.33/(CFH  toner*.0712) 

430  Rlg(2)=Rlg 

440  A=« 

450  SOSUB  COUNT 

460  Ntur(2)=CT 

470  XN=-(LST  (Ntur tl) /Htur (2) ) /LST  (Rlg(l)/Rlg(2D) 

480  C=Ntur(l)*Rlg(l)AXN 

490  ! 

500  PRINT 

510  PRINT  '  COOLING  TOWER  FILL  CHARACTERISTICS" 

520  PRINT  '         ■ 

530  PRINT 

540  PRINT  USIN6  "13I,AAA,DD.DDDDD,!51,AA,DD.DDDDD"  ;  "XN=',XN ,"C=",C 

550  PRINT 

360  PRINT 

570  END 

580  ! 

590  COUNT:  !       SUBROUTINE  TO  FIND  KaV/L  FOR  C.T. 

600  ! 

610  CT=0 

620  Yd )=.  1 

630  Y(2)=.3 

640  Y(3)=.2 

650  Y(4)=.3 

660  Teip-TWB 

670  SOSUB  HVTA 

680  HA=H 

690  T»=TWB+A 

700  FOR  1=1  TO  4 

710   TW=TWtY(I)*R 

720   HA=HA+Rlg*Y(I)*R 

730   Te»p=TW 

740   60SUB  HVTA 

750   HT=H 


SUBROUTINE  TO  FIND  THE  ENTHALPY  OF  SIS 


SUBROUTINE  TO  IDENTIFY  THE  COOLING  TQKER 


; TOMER*; * 


BODE  ";«ODEL* 


760      DH=HT-HA 

HO      CT=CTtl/DH 

730  NEJT  I 

790  CT=CT»R/4 

300  RETURN 

810  ! 

820  ! 

930  HVTA:  ! 

S40  PB=101325 

850  SDSUB  PSflT 

840  HS=.622*(PS/(PB-PS)) 

870  H=. 24*Te«p+HS* I 1041+. 444*Te«p ) 

880  RETURN 

890  ! 

900  TOKERID: 

910  ! 

920  PRINT  ' 

?30  PRINT 

940  PRINT  ' 

950  PRINT  ■ 

960  PRINT 

970  PRINT  " 

981  PRINT  ■ 

982  PRINT  " 
990  PRINT  ■ 
1020  PRINT  ■ 
1030  PRINT  ■ 
1040  PRINT  • 
1050  PRINT 
1060  PRINT 
1070  RETURN 
1080  ! 
1090  ! 

1100  PSAT:  ! 
1110  ! 
1120  All, 11=6.41533947 
1130  A(t,2)=4.41542599 
1140  A (1,31  =6. 44302266 
1130  A(l.4)=6. 53947833 
1160  A(2,l)=.0821033478 
1170  A(2, 21=. 0724023398 
1130  A(2,3)=.06988O3796 
1190  A(2,4I=,0656725584 
1200  A(3, 11=-. 00034006899! 
1210  A(3, 21=-. 000264361973 
1220  A(3,31=-2.0606744496E-4 
1230  A(3, 41=-. 000160854949 
1240  ! 

1250  TT=(Te«n-321*5/9 

1260  IF  TT>22.0322  THEN  GOTO  1270  ELSE  GOTO  1310 

1270   IF  TK46.0304  THEN  J=3  ELSE  SOTO  1290 

1280   GOTO  1370 

1290   IF  TT<72.0476  THEN  J=4  ELSE  GOTO  1350 

1300   GOTO  1370 

1310  IF  TT). 00895  THEN  J=2  ELSE  GOTO  1330 

1320   GOTO  1370 

1330   IF  TT>-19.94  THEN  J=l  ELSE  GOTO  1350 

1340   GOTO  1370 

1350  PRINT  '  THE  TEMPERATURE  FALLS  OUTSIDE  OF  OUR  RANGE 

1360  PAUSE 

1370  PS=EXP  (A(l,Jltm(A(2,JltTT»A(3,JD) 

1380  RETURN 

1390  ! 
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COOLINB  TONER  SPECIFICATIONS' 


DESI6N  SET-BULB  TEMPERATURE 
DESIGN  RANGE  OF  THE  C.T. 
DESIGN  AIR  FLOW  RATE 
FIRST  DEGIGN  WATER  FLOK  RATE 

FIRST  DESIGN  APPROACH  OF  THE  C.T. 

SECOND  DESIGN  APPROACH  OF  THE  C.T. 

SECOND  DESIGN  HATER  aQH  RATE  OF  THE  C.T. 


SUBROUTINE  TO  FIND  SATURATION  PRESSURE 


tdeg.F) 

=";THB 

(deg.F) 
(c?«) 

=-;R 

=";CFN  taner 

Iqpil 
Ideg.F) 

=';6P?ir 

=';A1 

(deg.F) 

=";A2 

(gps) 

=';SP!12 
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ABSTRACT 

Refrigeration/dehumidification  systems  operate  the  majority  of  the 
time  at  conditions  less  severe  than  those  of  design.  A  generalized  model 
describing  reciprocating  vapor-compression  ref rigeraton/dehuraidification 
system  performance  was  developed  with  cylinder  unloading  utilized  as  the 
basis  for  system  capacity  control  at  partial  load  operation.  Normalized 
parameters  were  used  where  possible  to  reduce  the  number  of  parameters 
required  to  accurately  represent  equipment  size  and  operational 
characteristics.  The  procedures  outlined  allow  the  model  to  reproduce 
catalog  data  (+3-1/2%)  for  a  wide  range  of  equipment  sizes.  Because  of  the 
generalized  nature  of  the  model,  most  of  the  reciprocating 
refrigeration/dehumidification  systems  available  are  covered  by  the  model. 


